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.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. I. The del-operator MEs in a two-shell composite Gelfand-Paldus basis

This is the first in a series of three articles which aimed to derive the matrix elements of the U(2n) generators in a multishell spin-orbit basis. This is a basis appropriate to many-electron systems which have a natural partitioning of the orbital space and where also spin-dependent terms are included in the Hamiltonian. The method is based on a new spin-dependent unitary group approach to the many-electron correlation problem due to Gould and Paldus [M. D. Gould and J. Paldus, J. Chem. Phys. 92, 7394, (1990)]. In this approach, the matrix elements of the U(2n) generators in the U(n) x U(2)-adapted electronic Gelfand basis are determined by the matrix elements of a single Ll(n) adjoint tensor operator called the del-operator, denoted by Delta(j)(i) (1 less than or equal to i, j less than or equal to n). Delta or del is a polynomial of degree two in the U(n) matrix E = [E-j(i)]. The approach of Gould and Paldus is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. Hence, to generalize this approach, we need to obtain formulas for the complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The nonzero shift coefficients are uniquely determined and may he evaluated by the methods of Gould et al. [see the above reference]. In this article, we define zero-shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis which are appropriate to the many-electron problem. By definition, these are proportional to the corresponding two-shell del-operator matrix elements, and it is shown that the Racah factorization lemma applies. Formulas for these coefficients are then obtained by application of the Racah factorization lemma. The zero-shift adjoint reduced Wigner coefficients required for this procedure are evaluated first. All these coefficients are needed later for the multishell case, which leads directly to the two-shell del-operator matrix elements. Finally, we discuss an application to charge and spin densities in a two-shell molecular system. (C) 1998 John Wiley & Sons.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

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.

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.

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.

Administration of growth hormone or IGF-I to pregnant rats on a reduced diet throughout pregnancy does not prevent fetal intrauterine growth retardation and elevated blood pressure in adult offspring

Increasing evidence from human epidemiological studies suggests that poor growth before birth is associated with postnatal growth retardation and the development of cardiovascular disease in adulthood. We have shown previously that nutritional deprivation in the pregnant rat leads to intrauterine growth retardation (IUGR), postnatal growth failure, changes in the endocrine parameters of the somatotrophic axis, and to increased blood pressure in later life. In the present study, we investigated whether administration of insulin-like growth factor-I (IGF-I) or bovine growth hormone (GH) during pregnancy could prevent IUGR and/or alter long-term outcome. Dams h-om day 1 of pregnancy throughout gestation received a diet of nd libitum available food or a restricted dietary intake of 30% of ad libitum fed dams. From day 10 of gestation, dams were treated for 10 days with three times daily subcutaneous injections of saline (100 mu l), IGF-I (2 mu g/g body weight) or GH (2 mu g/g body weight). Maternal weight gain was significantly increased (P

Can tissue drug concentrations be monitored by microdialysis during or after isolated limb perfusion for melanoma treatment?

Isolated limb perfusion (ILP) with melphalan is used to treat recurrent melanoma. This study aimed to develop a microdialysis technique for melphalan tissue concentration measurement during ILP. The effects of melphalan concentration (50-600 mu g/ml), microdialysis flow rate (0.55-17.5 mu l/min), probe length (5-50 mm) and temperature (25-41.5 degrees C) on in vitro recovery were studied. In addition, in vivo recovery was measured in rat hindlimbs perfused with melphalan using 50 mm microdialysis probes implanted subcutaneously and into muscle. Both dialysate and tissue sample melphalan concentrations were determined by high performance liquid chromatography. The in vitro recovery of melphalan was not affected by melphalan concentration or temperature, but increased with probe length and decreased with flow rate. The melphalan concentrations in subcutaneous and muscle dialysates were not significantly different. A linear relationship was found between tissue dialysate concentrations and actual tissue concentrations of melphalan (r(2) = 0.97). Microdialysis is a potential method for tissue drug monitoring which may assist in the efficacious use of cytotoxics in human ILP. (C) 2000 Lippincott Williams & Wilkins.

Glucocorticoid receptors in human preadipocytes: regional and gender differences

Glucocorticoid excess causes visceral obesity and its accompanying insulin resistance, dyslipidemia and hypertension. Glucocorticoids enhance preadipocyte (PA) differentiation and increase their aromatase activity (oestrogen production) and there is regional variability in these PA processes. Therefore, we studied human PAs for the presence of, and any regional or gender differences in, glucocorticoid receptors (GRs). Confluent subcultured human subcutaneous (Sc) and visceral (Vis) PAs from both genders contained GRs as assessed by GR gene expression and specific glucocorticoid (dexamethasone) binding. The dissociation constant was similar to that of other human cells and there was no difference between Sc and Vis sites or between males and females. There was significantly less GR mRNA in Vis PAs compared with Sc PAs in females (P=0.008) but not in males. There was less glucocorticoid binding in Vis compared with Sc PAs in females, measured by maximal binding capacity (P=0.035) or single saturating dose glucocorticoid binding (Bssd) (P=0.019). There was no regional difference in specific glucocorticoid binding in males. There was a gender difference with fewer GRs in Vis PAs in females compared with males measured by Bssd (P=0.006). In summary, GRs are present in human PAs. There is a lower GR density in Vis compared with Sc PAs in females, and females have fewer GRs in Vis PAs compared with males. These differences are likely to affect regional aromatase activity and to contribute to the smaller visceral fat mass in females compared with males.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Residues of zeta-cypermethrin in bovine tissues and milk following pour-on and spray application

The depletion of zeta-cypermethrin residues in bovine tissues and milk was studied. Beef cattle were treated three times at 3-week intervals with 1 ml 10 kg(-1) body weight of a 25 g litre(-1) or 50 g litre(-1) pour-on formulation (2.5 and 5.0 mg zeta-cypermethrin kg(-1) body weight) or 100 mg kg(-1) spray to simulate a likely worst-case treatment regime. Friesian and Jersey dairy cows were treated once with 2.5 mg zeta-cypermethrin kg(-1) in a pour-on formulation. Muscle, liver and kidney residue concentrations were generally less than the limit of detection (LOD = 0.01 mg kg(-1)). Residues in renal-fat and back-fat samples from animals treated with 2.5 mg kg(-1) all exceeded the limit of quantitation (LOQ = 0.05 mg kg(-1)), peaking at 10 days after treatment. Only two of five kidney fat samples were above the LOQ after 34 days, but none of the back-fat samples exceeded the LOQ at 28 days after treatment. Following spray treatments, fat residues were detectable in some animals but were below the LOQ at all sampling intervals. Zeta-cypermethrin was quantifiable (LOQ = 0.01 mg kg(-1)) in only one whole-milk sample from the Friesian cows (0.015 mg kg(-1), 2 days after treatment). In whole milk from Jersey cows, the mean concentration of zeta-cypermethrin peaked 1 day after treatment, at 0.015 mg kg(-1), and the highest individual sample concentration was 0.025 mg kg(-1) at 3 days after treatment. Residues in milk were not quantifiable beginning 4 days after treatment. The mean concentrations of zeta-cypermethrin in milk fat from Friesian and Jersey cows peaked two days after treatment at 0.197 mg kg(-1) and 0.377 mg kg(-1), respectively, and the highest individual sample concentrations were 2 days after treatment at 0.47 mg kg(-1) and 0.98 mg kg(-1), respectively. (C) 2001 Society of Chemical Industry.

A method for the design of MRI radiofrequency coils based on triangular and pulse basis functions

This paper presents a numerical technique for the design of an RF coil for asymmetric magnetic resonance imaging (MRI) systems. The formulation is based on an inverse approach where the cylindrical surface currents are expressed in terms of a combination of sub-domain basis functions: triangular and pulse functions. With the homogeneous transverse magnetic field specified in a spherical region, a functional method is applied to obtain the unknown current coefficients. The current distribution is then transformed to a conductor pattern by use of a stream function technique. Preliminary MR images acquired using a prototype RF coil are presented and validate the design method. (C) 2002 Elsevier Science B.V. All rights reserved.

Relationship of interleukin-6 and tumor necrosis factor-alpha with muscle mass and muscle strength in elderly men and women: The health ABC study

Background. A decline in muscle mass and muscle strength characterizes normal aging. As clinical and animal studies show it relationship between higher cytokine levels and low muscle mass, the aim of this study was to investigate whether markers, of inflammation are associated with muscle mass and strength in well-functioning elderly persons. Methods. We Used baseline data (1997-1998) of the Health, Aging, and Body Composition (Health ABC) Study on 3075 black and white men and women aged 70-79 years. Midthigh muscle cross-sectional area (computed tomography), appendicular muscle mass (dual-energy x-ray ab absorptiometry). isokinetic knee extensor strength (KinCom). and isometric inip strength were measured. plasma levels of interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-alpha) were assessed by enzyme-linked immunosorbent assay (ELISA). Results. Higher cytokine levels were generally associated with lower muscle mass and lower muscle strength. The most consistent relationship across the gender and race groups was observed for IL-6 and grip strength: per SD increase in IL-6, grip strength was 1.1 to 2.4 kg lower (p < .05) after adjustment for age, clinic Site. health status, medications, physical activity. smoking. height. and body fat. Ail overall measure of elevated cytokine level was created by combining the levels of IL-6 and TNF-alpha. With the exception of white men, elderly persons having high levels of IL-6 (> 1.80 pg/ml) as well as high levels of TNF-alpha (> 3.20 pg/ml) had a smaller muscle area, less appendicular mass. a lower knee extensor strength. and a lower grip strength compared to those with low levels of both cytokines. Conclusions. Higher plasma concentrations of IL-6 and TNF-alpha are associated with lower muscle mass and lower muscle strength in well-functioning older men and women. Higher cytokine levels. as often observed in healthy older persons. may contribute to the loss Of muscle mass and strength that accompanies aging.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.