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.

Coupling of a Jahn-Teller pseudorotation with a hindered internal rotation in an isolated molecule: 9-hydroxytriptycene

The irregular vibronic structure in the S-1<--S-0 resonant two-photon ionization (R2PI) spectrum of supersonically cooled triptycene is a result of a classic Exe Jahn-Teller effect [A. Furlan et al., J. Chem. Phys. 96, 7306 (1992)]. This is well characterized and can be used as an effective probe of intramolecular perturbations. Here we examine the S-1<--S-0 R2PI spectrum of 9-hydroxytriptycene and the fluorescence from various excited state vibronic levels. In this system the pseudorotation of the Jahn-Teller vibration is strongly coupled to the torsional motion of the bridgehead hydroxy group. This torsional motion results in a tunneling splitting in both the ground and excited states. The population of the upper level in the ground electronic state results in additional vibronic transitions becoming symmetry allowed in the R2PI spectrum that are forbidden in the bare triptycene molecule. The assignment of the R2PI and fluorescence spectra allows the potential energy surfaces of these vibrational modes to be accurately quantified. The full C-3v vibronic point group must be used to interpret the spectra. The time scale of the internal rotation of the-OH group and the butterfly flapping of the Jahn-Teller pseudorotation are of similar magnitude. The tunneling between the nine minima on the three dimensional potential energy surface is such that the Jahn-Teller pseudorotation occurs in concert with the-OH internal rotation. The Berry phase that is acquired during this motion is discussed. The simple physical picture emerges of the angle between two of the three benzene moieties opening in three equivalent ways in the S-1 electronic state. This geometry follows the position of the hydroxy group, which preferentially orients itself to point between these two rings. (C) 1998 American Institute of Physics. [S0021-9606(98)02348-4].

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Simulation and optimization of heated, inviscid flows in scramjet ducts

A space-marching code for the simulation and optimization of inviscid supersonic flow in three dimensions is described. The now in a scramjet module with a relatively complex three-dimensional geometry is examined and wall-pressure estimates are compared with experimental data. Given that viscous effects are not presently included, the comparison is reasonable. The thermodynamic compromise of adding heat in a diverging combustor is also examined. The code is then used to optimize the shape of a thrust surface for a simpler (box-section) scramjet module in the presence of uniform and nonuniform heat distributions. The optimum two-dimensional profiles for the thrust surface are obtained via a perturbation procedure that requires about 30-50 now solutions. It is found that the final shapes are fairly insensitive to the details of the heat distribution.