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.

Potential errors in the application of mixture theory to multifrequency bioelectrical impedance analysis

Potential errors in the application of mixture theory to the analysis of multiple-frequency bioelectrical impedance data for the determination of body fluid volumes are assessed. Potential sources of error include: conductive length; tissue fluid resistivity; body density; weight and technical errors of measurement. Inclusion of inaccurate estimates of body density and weight introduce errors of typically < +/-3% but incorrect assumptions regarding conductive length or fluid resistivities may each incur errors of up to 20%.

When oceans meet: a teleost shows secondary intergradation at an Indian-Pacific interface

The Indo-West Pacific is characterized by extraordinary marine species diversity. The evolutionary mechanisms responsible for generating this diversity remain puzzling, but are often linked to Pleistocene sea level fluctuations. The impact of these sea level changes on the population genetic architecture of the estuarine fish Lates calcarifer are investigated via a natural experiment in a region of the Indo-West Pacific known to have undergone considerable change during the Pleistocene. L. calcarifer, a coastline-restricted catadromous teleost, provides an excellent model for studying the effects of sea level change as its habitat requirements potentially make it sensitive to the region's physical history. Evidence was found for a large phylogenetic break (4% mtDNA control region; 0.47% ATPase 6 and 8) either side of the Torres Strait, which separates the Western Pacific and Indian Oceans, although some mixing of the clades was evident. This suggests clinal secondary introgression of the clades via contemporary gene flow. Further, populations on Australia's east coast appear to have passed through a bottleneck. This was linked to the historical drying of the Great Barrier Reef coastal lagoon, which resulted in a significant loss of habitat and forced retreat into isolated refugia. These results suggest that historical eustatic changes have left a significant imprint on the molecular diversity within marine species as well as among them in the Indo-West Pacific.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Response (re-)programming in aging: A kinematic analysis

Background. Age-related motor slowing may reflect either motor programming deficits, poorer movement execution, or mere strategic preferences for online guidance of movement. We controlled such preferences, limiting the extent to which movements could be programmed. Methods. Twenty-four young and 24 older adults performed a line drawing task that allowed movements to he prepared in advance in one case (i.e., cue initially available indicating target location) and not in another (i.e., no cue initially available as to target location). Participants connected large or small targets illuminated by light-emitting diodes upon a graphics tablet that sampled pen tip position at 200 Hz. Results. Older adults had a disproportionate difficulty initiating movement when prevented from programming in advance. Older adults produced slower, less efficient movements, particularly when prevented from programming under greater precision requirements. Conclusions. The slower movements of older adults do not simply reflect a preference for online control, as older adults have less efficient movements when forced to reprogram their movements. Age-related motor slowing kinematically resembles that seen in patients with cerebellar dysfunction.

Epidermal iontophoresis: II. Application of the ionic mobility-pore model to the transport of local anesthetics

Purpose, An in vitro study was carried out to determine the iontophoretic permeability of local anesthetics through human epidermis. The relationship between physicochemical structure and the permeability of these solutes was then examined using an ionic mobility-pore model developed to define quantitative relationships. Methods. The iontophoretic permeability of both ester-type anesthetics (procaine, butacaine, tetracaine) and amide-type anesthetics (prilocaine, mepivacaine, lidocaine, bupivacaine, etidocaine, cinchocaine) were determined through excised human epidermis over 2 hrs using a constant d.c. current and Ag/AgCl electrodes. Individual ion mobilities were determined from conductivity measurements in aqueous solutions. Multiple stepwise regression was applied to interrelate the iontophoretic permeability of the solutes with their physical properties to examine the appropriateness of the ionic mobility-pore model and to determine the best predictor of iontophoretic permeability of the local anesthetics. Results. The logarithm of the iontophoretic permeability coefficient (log PCj,iont) for local anesthetics was directly related to the log ionic mobility and MW for the free volume form of the model when other conditions are held constant. Multiple linear regressions confirmed that log PCj,iont was best defined by ionic mobility (and its determinants: conductivity, pK(a) and MW) and MW. Conclusions. Our results suggest that of the properties studied, the best predictors of iontophoretic transport of local anesthetics are ionic mobility (or pK(a)) and molecular size. These predictions are consistent with the ionic mobility pore model determined by the mobility of ions in the aqueous solution, the total current, epidermal permselectivity and other factors as defined by the model.

A state-contingent production approach to principal-agent problems with an application to point-source pollution control

Two basic representations of principal-agent relationships, the 'state-space' and 'parameterized distribution' formulations, have emerged. Although the state-space formulation appears more natural, analytical studies using this formulation have had limited success. This paper develops a state-space formulation of the moral-hazard problem using a general representation of production under uncertainty. A closed-form solution for the agency-cost problem is derived. Comparative-static results are deduced. Next we solve the principal's problem of selecting the optimal output given the agency-cost function. The analysis is applied to the problem of point-source pollution control. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

A relative gradient theory for layered materials

It is possible to remedy certain difficulties with the description of short wave length phenomena and interfacial slip in standard models of a laminated material by considering the bending stiffness of the layers. If the couple or moment stresses are assumed to be proportional to the relative deformation gradient, then the bending effect disappears for vanishing interface slip, and the model correctly reduces to an isotropic standard continuum. In earlier Cosserat-type models this was not the case. Laminated materials of the kind considered here occur naturally as layered rock, or at a different scale, in synthetic layered materials and composites. Similarities to the situation in regular dislocation structures with couple stresses, also make these ideas relevant to single slip in crystalline materials. Application of the theory to a one-dimensional model for layered beams demonstrates agreement with exact results at the extremes of zero and infinite interface stiffness. Moreover, comparison with finite element calculations confirm the accuracy of the prediction for intermediate interfacial stiffness.

Observation of a cubic phase in mesoporous silicate films grown at the air/water interface

X-Ray diffraction is reported from mesoporous silicate films grown at the air/water interface. The films were studied both as powdered films, and oriented on silicon or mica sheets. At early stages of growth we observe Bragg diffraction from a highly ordered cubic phase, with both long and short d-spacing peaks. We have assigned this as a discontinuous micellar Pm3n phase in which the silica is partly ordered. Later films retain only the known hexagonal p6m peaks and have lost any order both at short d-spacings and the longer d-spacing Bragg peaks characteristic of the cubic structure. The silica framework is considerably expanded from that in bulk amorphous silica, average Si Si distances are some 30% greater. Incorporation of glycerol or polyethylene glycol preserves the earlier cubic structure. To be consistent with earlier, in situ, X-ray and neutron reflectivity data we infer that both structures are produced after a phase transition from a less-ordered him structure late in the induction phase. The structural relations between the film Pm3n and p6m phase(s) and the known bulk SBA-1 and MCM-41 phases are briefly discussed.

Determination of chemical shielding tensor of an indole carbon and application of tryptophan orientation of a membrane peptide

Using tryptophan C-13-enriched at the C-4 (C epsilon(3)) of the indole, the orientation of the C epsilon(3) chemical shift tensor relative to the C epsilon(3)-H dipolar axis was determined from the C-13 chemical shift/C-13-H-1 dipolar 2D NMR powder pattern. The principal values obtained were 208, 137 and 15 ppm with sigma(33) perpendicular to the indole plane, and sigma(11) (least shielded direction) 5 degrees off the C epsilon(3)-H bond toward C xi(3). The side off the C epsilon(3)-H bond was determined by comparing the reduced chemical shift anisotropies obtained by solid-state NMR and from molecular dynamics calculations of [4-C-13] tryptophans in gramicidin A aligned in phospholipid membranes. (C) 1999 Elsevier Science B.V. All rights reserved.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Targeting local tissues by transdermal application: Understanding drug physicochemical properties that best exploit protein binding and blood flow effects

The targeting of topically applied drug molecules into tissues below a site of application requires an understanding of the complex interrelationships between the drug, its formulation, the barrier properties of the skin, and the physiological processes occurring below the skin that are responsible for drug clearance from the site, tissue, and/or systemic distribution and eventual elimination. There is still a certain amount of controversy over the ability of topically applied drugs to penetrate into deeper tissues by diffusion or whether this occurs by redistribution in the systemic circulation. The major focus of our work in this area has been in determining how changes in drug structure and physicochemical properties, such as protein binding and lipophilicity, affect drug clearance into the local dermal microcirculation and lymphatics, as well as subsequent distribution into deeper tissues below an application site. The present study outlines our recent thinking on the drug molecule optimal physical attributes, in terms of plasma and tissue partitioning behaviour, that offer the greatest potential for deep tissue targeting. Drug Dev. Res. 46:309-315, 1999. (C) 1999 Wiley-Liss, Inc.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Incorporating phase retardation effects into radiofrequency resonator models: application to magnetic resonance microscopy

A method is presented for including path propagation effects into models of radiofrequency resonators for use in magnetic resonance imaging. The method is based on the use of Helmholtz retarded potentials and extends our previous work on current density models of resonators based on novel inverse finite Hilbert transform solutions to the requisite integral equations. Radiofrequency phase retardation effects are most pronounced at high field strengths (frequencies) as are static field perturbations due to the magnetic materials in the resonators themselves. Both of these effects are investigated and a novel resonator structure presented for use in magnetic resonance microscopy.