Hydraulic Design of Stepped Spillways and Downstream Energy Dissipators for Embankment Dams

In recent years, the design flows of many dams were re-evaluated, often resulting in discharges larger than the original design. In many cases, the occurrence of the revised flows could result in dam overtopping because of insufficient storage and spillway capacity. An experimental study was conducted herein to gain a better understanding of the flow properties in stepped chutes with slopes typical of embankment dams. The work was based upon a Froude similitude in large-size experimental facilities. A total of 10 configurations were tested including smooth steps, steps equipped with devices to enhance energy dissipation and rough steps. The present results yield a new design procedure. The design method includes some key issues not foreseen in prior studies : e.g., gradually varied flow, type of flow regime, flow resistance. It is believed that the outcomes are valid for a wide range of chute geometry and flow conditions typical of embankment chutes.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

New calibration technique for multiple-component stress wave force balances

Force measurement in hypervelocity expansion tubes is not possible using conventional techniques. The stress wave force balance technique can be applied in expansion tubes to measure forces despite the short test times involved. This paper presents a new calibration technique for multiple-component stress wave force balances where an impulse response created using a load distribution is required and no orthogonal surfaces on the model exist.. This new technique relies on the tensorial superposition of single-component impulse responses analogous to the vectorial superposition of the calibration loads. The example presented here is that of a scale model of the Mars Pathfinder, but the technique is applicable to any geometry and may be useful for cases where orthogonal loads cannot be applied.

Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Solution structure and proposed binding mechanism of a novel potassium channel toxin kappa-conotoxin PVIIA

Background: kappa-PVIIA is a 27-residue polypeptide isolated from the venom of Conus purpurascens and is the first member of a new class of conotoxins that block potassium channels. By comparison to other ion channels of eukaryotic cell membranes, voltage-sensitive potassium channels are relatively simple and methodology has been developed for mapping their interactions with small-peptide toxins, PVIIA, therefore, is a valuable new probe of potassium channel structure. This study of the solution structure and mode of channel binding of PVIIA forms the basis for mapping the interacting residues at the conotoxin-ion channel interface. Results: The three-dimensional structure of PVIIA resembles the triple-stranded beta sheet/cystine-knot motif formed by a number of toxic and inhibitory peptides. Subtle structural differences, predominantly in loops 2 and 4, are observed between PVIIA and other conotoxins with similar structural frameworks, however. Electrophysiological binding data suggest that PVIIA blocks channel currents by binding in a voltage-sensitive manner to the external vestibule and occluding the pore, Comparison of the electrostatic surface of PVIIA with that of the well-characterised potassium channel blocker charybdotoxin suggests a likely binding orientation for PVIIA, Conclusions: Although the structure of PVIIA is considerably different to that of the alpha K scorpion toxins, it has a similar mechanism of channel blockade. On the basis of a comparison of the structures of PVIIA and charybdotoxin, we suggest that Lys19 of PVIIA is the residue which is responsible for physically occluding the pore of the potassium channel.

A novel, open access, elliptical cross-section magnet for paediatric MRI

Almost all clinical magnetic resonance imaging systems are based on circular cross-section magnets. Recent advances in elliptical cross-section RF probe and gradient coil hardware raise the question of the possibility of using elliptical cross-section magnet systems, This paper presents a methodology for calculating rapidly the magnetic fields generated by a multi-turn coil of elliptical cross-section and incorporates this in a stochastic optimization method for magnet design, An open magnet system of elliptical cross-section is designed that both reduces the claustrophobia for the patients and allows ready access by attending physicians, The magnet system is optimized for paediatric use, The coil geometry produced by the optimization method has several novel features.

Two dimensional computational fluid dynamic models for waste stabilisation ponds

Traditional waste stabilisation pond (WSP) models encounter problems predicting pond performance because they cannot account for the influence of pond features, such as inlet structure or pond geometry, on fluid hydrodynamics. In this study, two dimensional (2-D) computational fluid dynamics (CFD) models were compared to experimental residence time distributions (RTD) from literature. In one of the-three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, flow patterns in the other two geometries were not well described due to the difficulty of representing the three dimensional (3-D) experimental inlet in the 2-D CFD model, and the sensitivity of the model results to the assumptions used to characterise the inlet. Neither a velocity similarity nor geometric similarity approach to inlet representation in 2-D gave results correlating with experimental data. However. it was shown that 2-D CFD models were not affected by changes in values of model parameters which are difficult to predict, particularly the turbulent inlet conditions. This work suggests that 2-D CFD models cannot be used a priori to give an adequate description of the hydrodynamic patterns in WSP. (C) 1998 Elsevier Science Ltd. All rights reserved.

Modification of MCM-41 by surface silylation with trimethylchlorosilane and adsorption study

Siliceous MCM-41 samples were modified by silylation using trimethylchlorosilane (TMCS). The surface coverage of functional groups was studied systematically in this work. The role of surface silanol groups during modification was evaluated using techniques of FTIR and Si-29 CP/MAS NMR. Adsorption of water and benzene on samples of various hydrophobicities was measured and compared. It was found that the maximum degree of surface attachments of trimethylsilyl (TMS) groups was about 85%, corresponding to the density of TMS groups of 1.9 per nm(2). The degree of silylation is found to linearly increase with increasing pre-outgassing temperature prior to silylation. A few types of silanol groups exist on MCM-41 surfaces, among which both free and geminal ones are responsible for active silylation. Results of water adsorption show that aluminosilicate MCM-41 materials are more or less hydrophilic, giving a type IV isotherm, similar to that of nitrogen adsorption, whereas siliceous MCM-41 are hydrophobic, exhibiting a type V adsorption isotherm. The fully silylated Si-MCM-41 samples are more hydrophobic giving a type III adsorption isotherm. Benzene adsorption on all MCM-41 samples shows type IV isotherms regardless of the surface chemistry. Capillary condensation occurs at a higher relative pressure for the silylated MCM-41 than that for the unsilylated sample, though the pore diameter was found reduced markedly by silylation. This is thought attributed to the diffusion constriction posed by the attached TMS groups. The results show that the surface chemistry plays an important role in water adsorption, whereas benzene adsorption is predominantly determined by the pore geometry of MCM-41.

Copper(II) complexes of mono- and di-nucleating hexaamines

The potentially hexadentate polyamines N,N,N',N'-tetrakis(2-aminoethyl)ethane-1,2-diamine (L-1) and the octamethylated analogue N,N,N',N'-tetrakis(2-dimethylaminoethyl)ethane 1,2-diamine (L-2) have been complexed with copper(II) and the crystal structures of their complexes determined. A trigonal-bipyramidal co-ordination geometry for [Cu(HL1)][ClO4](3) was found where one aminoethyl arm is not co-ordinated. By contrast, a dinuclear structure of formula [(H2O)Cu(L-2)Cu(OH)](3+) was determined for the N-methylated analogue, where the hexaamine acts as a bridging ligand between the two square-pyramidal metal centres. Electronic and EPR spectroscopy are both consistent with these structures being maintained in solution.

Thermochromism and structure of piperazinium tetrachlorocuprate(II) complexes

The tetrachlorocuprate(II) ion can crystallize in two different structures with the piperazinium dication (pipzH(2)). Both structures contain discrete CuCl42- species. A yellow compound (pipzH(2))[CuCl4]. 2H(2)O (1) is monoclinic (C2/c, Z = 4, a = 10.538(3) Angstrom, b = 7.4312(5) Angstrom, c = 17.281(4) Angstrom, beta = 111.900(10)degrees) and contains the CuCl42- ion as a distorted tetrahedron. A green compound (pipzH(2))(2)[CuCl4]. Cl-2. 3H(2)O (2) is triclinic (P (1) over bar, Z = 2, a = 9.264(3) Angstrom, b = 10.447(2) Angstrom, c = 11.366(2) Angstrom, alpha = 68.38 degrees, beta = 82.86(2)degrees, gamma = 83.05(2)degrees) and contains the CuCl42- ion with a square planar geometry. This latter compound shows thermo/photochromism, changing from green to yellow upon heating or laser irradiation.

Spectroscopic identification of a dinuclear metal centre in manganese(II)-activated aminopeptidase P from Escherichia coli: implications for human prolidase

Electron paramagnetic resonance (EPR) spectra and X-ray absorption (EXAFS and XANES) data have been recorded for the manganese enzyme aminopeptidase P (AMPP, PepP protein) from Escherichia coli. The biological function of the protein, a tetramer of 50-kDa subunits, is the hydrolysis of N-terminal Xaa-Pro peptide bonds. Activity assays confirm that the enzyme is activated by treatment with Mn2+. The EPR spectrum of Mn2+-activated AMPP at liquid-He temperature is characteristic of an exchange-coupled dinuclear Mn(II) site, the Mn-Mn separation calculated from the zero-field splitting D of the quintet state being 3.5 (+/- 0.1) Angstrom. In the X-ray absorption spectrum of Mn2+-activated AMPP at the Mn K edge, the near-edge features are consistent with octahedrally coordinated Mn atoms in oxidation state +2. EXAFS data, limited to k less than or equal to 12 Angstrom(-1) by traces of Fe in the protein, are consistent with a single coordination shell occupied predominantly by O donor atoms at an average Mn-ligand distance of 2.15 Angstrom, but the possibility of a mixture of O and N donor atoms is not excluded. The Mn-Mn interaction at 3.5 Angstrom, is not detected in the EXAFS, probably due to destructive interference from light outer-shell atoms. The biological function, amino acid sequence and metal-ion dependence of E. coli AMPP are closely related to those of human prolidase, an enzyme that specifically cleaves Xaa-Pro dipeptides. Mutations that lead to human prolidase deficiency and clinical symptoms have been identified. Several known inhibitors of prolidase also inhibit AMPP. When these inhibitors are added to Mn2+-activated AMPP, the EPR spectrum and EXAFS remain unchanged. It can be inferred that the inhibitors either do not bind directly to the Mn centres, or substitute for existing Mn ligands without a significant change in donor atoms or coordination geometry. The conclusions from the spectroscopic measurements on AMPP have been verified by, and complement, a recent crystal structure analysis.

Three-dimensional solution structure of alpha-conotoxin MII by NMR spectroscopy: Effects of solution environment on helicity

alpha-Conotoxin MII, a 16-residue polypeptide from the venom of the piscivorous cone snail Conus magus, is a potent and highly specific blocker of mammalian neuronal nicotinic acetylcholine receptors composed of alpha 3 beta 2 subunits. The role of this receptor type in the modulation of neurotransmitter release and its relevance to the problems of addiction and psychosis emphasize the importance of a structural understanding of the mode of interaction of MII with the alpha 3 beta 2 interface. Here we describe the three-dimensional solution structure of MIT determined using 2D H-1 NMR spectroscopy. Structural restraints consisting of 376 interproton distances inferred from NOEs and 12 dihedral restraints derived from spin-spin coupling constants were used as input for simulated annealing calculations and energy minimization in the program X-PLOR. The final set of 20 structures is exceptionally well-defined with mean pairwise rms differences over the whole molecule of 0.07 Angstrom for the backbone atoms and 0.34 Angstrom for all heavy atoms. MII adopts a compact structure incorporating a central segment of alpha-helix and beta-turns at the N- and C-termini. The molecule is stabilized by two disulfide bonds, which provide cross-links between the N-terminus and both the middle and C-terminus of the structure. The susceptibility of the structure to conformational change was examined using several different solvent conditions. While the global fold of MII remains the same, the structure is stabilized in a more hydrophobic environment provided by the addition of acetonitrile or trifluoroethanol to the aqueous solution. The distribution of amino acid side chains in MII creates distinct hydrophobic and polar patches on its surface that may be important for the specific interaction with the alpha 3 beta 2 neuronal nAChR. A comparison of the structure of MII with other neuronal-specific alpha-conotoxins provides insights into their mode of interaction with these receptors.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.