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.

Characterization of calcium-activated chloride channels in patches excised from the dendritic knob of mammalian olfactory receptor neurons

We investigated the properties of calcium-activated chloride channels in inside-out membrane patches from the dendritic knobs of acutely dissociated rat olfactory receptor neurons. Patches typically contained large calcium-activated currents, with total conductances in the range 30-75 nS. The dose response curve for calcium exhibited an EC50 of about 26 mu M. In symmetrical NaCl solutions, the current-voltage relationship reversed at 0 mV and was linear between -80 and +70 mV. When the intracellular NaCl concentration was progressively reduced from 150 to 25 mM, the reversal potential changed in a manner consistent with a chloride-selective conductance. Indeed, modeling these data with the Goldman-Hodgkin-Katz equation revealed a P-Na/P-Cl of 0.034. The halide permeability sequence was P-Cl > P-F > P-I > P-Br indicating that permeation through the channel was dominated by ion binding sites with a high field strength. The channels were also permeable to the large organic anions, SCN-, acetate(-), and gluconate(-), with the permeability sequence P-Cl > P-SCN > gluconaie. Significant permeation to gluconate ions suggested that the channel pore had a minimum diameter of at least 5.8 Angstrom.

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.

Overflow characteristics of circular weirs: Effects of inflow conditions

The most common types of weirs are the broad-crested weir, the sharp-crested weir, the circular-crested weir, and nowadays, the ogee crest weir, Advantages of the cylindrical weir shape include the stable overflow pattern, the ease to pass floating debris, the simplicity of design compared to ogee crest design, and the associated lower costs. in this study, the writers describe new experiments of circular weir overflows, with eight cylinder sizes, for several weir heights and for five types of inflow conditions: partially developed inflow, fully developed inflow, upstream ramp, upstream undular hydraulic jump, and upstream (breaking) hydraulic jump. Within the range of the experiments, the cylinder size, the weir height DIR and the presence of an upstream ramp had no effect on the discharge coefficient, flow depth at crest, and energy dissipation. But the inflow conditions had substantial effects on the discharge characteristics and flow properties at the crest. Practically, the results indicate that discharge measurements with circular weirs are significantly affected by the upstream flow conditions.

Modeling ex vivo hematopoiesis using chemical engineering metaphors

Ex vivo hematopoiesis is increasingly used for clinical applications. Models of ex vivo hematopoiesis are required to better understand the complex dynamics and to optimize hematopoietic culture processes. A general mathematical modeling framework is developed which uses traditional chemical engineering metaphors to describe the complex hematopoietic dynamics. Tanks and tubular reactors are used to describe the (pseudo-) stochastic and deterministic elements of hematopoiesis, respectively. Cells at any point in the differentiation process can belong to either an immobilized, inert phase (quiescent cells) or a mobile, active phase (cycling cells). The model describes five processes: (1) flow (differentiation), (2) autocatalytic formation (growth),(3) degradation (death), (4) phase transition from immobilized to mobile phase (quiescent to cycling transition), and (5) phase transition from mobile to immobilized phase (cycling to quiescent transition). The modeling framework is illustrated with an example concerning the effect of TGF-beta 1 on erythropoiesis. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

NMR structure of C-terminally tagged gramicidin channels

A biotin group was covalently attached to the C terminus of gramicidin A (gA) through a linker arm comprising a glycine residue with either one (gAXB) or two caproyl groups (gAXXB). High-resolution two-dimensional NMR spectroscopy was used to determine the structure of these modified gA analogues and [Lys(16)]gramicidin A (gA-Lys) in sodium dodecyl-d(25) sulphate micelles. Gated gA ion channels based on linking a receptor group to these gA analogues have been used recently as a component in a sensing device. The conformations of the gA backbones and amino acid side chains of lysinated gA and biotinylated gA in detergent micelles were found to be almost identical to that of native gA, i.e. that of an N-terminal to N-terminal (head to head) dimer formed by two right-handed, single-stranded beta(6.3) helices. The biotin tail of the gAXB and gAXXB and the lysine extremity of gA-Lys appeared to lie outside the micelle. Thus it appears that the covalent attachment of functional groups to the C terminus of gA does not disrupt the peptide's helical configuration. Further, single channel measurements of all three gA analogues showed that functioning ion channels were preserved within a membrane environment. (C) 1999 Elsevier Science B.V. All rights reserved.

A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

A wavelet visible difference predictor

In this paper, we describe a model of the human visual system (HVS) based on the wavelet transform. This model is largely based on a previously proposed model, but has a number of modifications that make it more amenable to potential integration into a wavelet based image compression scheme. These modifications include the use of a separable wavelet transform instead of the cortex transform, the application of a wavelet contrast sensitivity function (CSP), and a simplified definition of subband contrast that allows us to predict noise visibility directly from wavelet coefficients. Initially, we outline the luminance, frequency, and masking sensitivities of the HVS and discuss how these can be incorporated into the wavelet transform. We then outline a number of limitations of the wavelet transform as a model of the HVS, namely the lack of translational invariance and poor orientation sensitivity. In order to investigate the efficacy of this wavelet based model, a wavelet visible difference predictor (WVDP) is described. The WVDP is then used to predict visible differences between an original and compressed (or noisy) image. Results are presented to emphasize the limitations of commonly used measures of image quality and to demonstrate the performance of the WVDP, The paper concludes with suggestions on bow the WVDP can be used to determine a visually optimal quantization strategy for wavelet coefficients and produce a quantitative measure of image quality.

Photolytic manipulation of [Ca2+](i) reveals slow kinetics of potassium channels underlying the afterhyperpolarization in hipppocampal pyramidal neurons

The identity of the potassium channel underlying the slow, apamin-insensitive component of the afterhyperpolarization current (sl(AHP)) remains unknown. We studied sl(AHP) in CA1 pyramidal neurons using simultaneous whole-cell recording, calcium fluorescence imaging, and flash photolysis of caged compounds. Intracellular calcium concentration ([Ca2+](i)) peaked earlier and decayed more rapidly than sl(AHP). Loading cells with low concentrations of the calcium chelator EGTA slowed the activation and decay of sl(AHP). In the presence of EGTA, intracellular calcium decayed with two time constants. When [Ca2+](i) was increased rapidly after photolysis of DM-Nitrophen, both apamin-sensitive and apamin-insensitive outward currents were activated. The apamin-sensitive current activated rapidly (<20 msec), whereas the apamin-insensitive current activated more slowly (180 msec). The apamin-insensitive current was reduced by application of serotonin and carbachol, confirming that it was caused by sl(AHP) channels. When [Ca2+](i) was decreased rapidly via photolysis of diazo-2, the decay of sl(AHP) was similar to control (1.7 sec). All results could be reproduced by a model potassium channel gated by calcium, suggesting that the channels underlying sl(AHP) have intrinsically slow kinetics because of their high affinity for calcium.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.