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.

From sabotage to camouflage: viral evasion of cytotoxic T lymphocyte and natural killer cell-mediated immunity

The outcome of a virus infection is strongly influenced by interactions between host immune defences and virus 'anti-defence' mechanisms. For many viruses, their continued survival depends on, the speed of their attach: their capacity to replicate and transmit to uninfected hosts prior to their elimination by an effective immune response. In contrast, the success of persistent viruses lies in their capacity for immunological subterfuge: the evasion of host defence mechanisms by either mutation (covered elsewhere in this issue, by Gould and Bangham, pp. 321-328) or interference with the action of host cellular proteins that are important components of the immune response. This review will focus on the strategies employed by persistent viruses against two formidable host defences against virus infection: the CD8+ cytotoxic T lymphocyte (CTL) and natural killer (NK) cell responses.

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.

Nest temperature and respiratory gases during natural incubation in the broad-shelled river turtle, Chelodina expansa (Testudinata : Chelidae)

Temperature was monitored in three natural nests, and oxygen and carbon dioxide partial pressure monitored in one natural nest of the broad-shelled river turtle, Chelodina expansa, throughout incubation. Nest temperature decreased after nest construction in autumn, remained low during winter and gradually increased in spring to a maximum in summer. In a nest where temperature was recorded every hour, temperature typically fluctuated through a 2 degrees C cycle on a daily basis throughout the entire incubation period, and the nest always heated faster than it cooled. Oxygen and carbon dioxide partial pressures in this nest were similar to soil oxygen and carbon dioxide partial pressures for the first 5 months of incubation, but nest respiratory gas tensions deviated from the surrounding soil over the last three months of incubation. Nest respiratory gas tensions were not greatly different from those in the atmosphere above the ground except after periods of rain. After heavy rain during the last 3 months of incubation the nest became moderately hypoxic (P-O2 similar to 100 Torr) and hypercapnic (P-CO2 similar to 50 Torr) for several successive days. These short periods of hypoxia and hypercapnia were not lethal.

Nitrogen relations of natural and disturbed plant communities in tropical Australia

Nitrogen relations of natural and disturbed tropical plant communities in northern Australia (Kakadu National Park) were studied. Plant and soil N characteristics suggested that differences in N source utilisation occur at community and species level. Leaf and xylem sap N concentrations of plants in different communities were correlated with the availability of inorganic soil N (NH4+ and NO3-). In general, rates of leaf NO3- assimilation were low. Even in communities with a higher N status, including deciduous monsoon forest, disturbed wetland, and a revegetated mine waste rock dump, levels of leaf nitrate reductase, xylem and leaf NO3 levels were considerably lower than those that have been reported for eutrophic communities. Although NO3- assimilation in escarpment and eucalypt woodlands, and wetland, was generally low, within these communities there was a suite of species that exhibited a greater capacity for NO3- assimilation. These high-NO3- species were mainly annuals, resprouting herbs or deciduous trees that had leaves with high N contents. Ficus, a high-NO3- species, was associated with soil exhibiting higher rates of net mineralisation and net nitrification. Low-NO3- species were evergreen perennials with low leaf N concentrations. A third group of plants, which assimilated NO3- (albeit at lower rates than the high-NO3- species), and had high-N leaves, were leguminous species. Acacia species, common in woodlands, had the highest leaf N contents of all woody species. Acacia species appeared to have the greatest potential to utilise the entire spectrum of available N sources. This versatility in N source utilisation may be important in relation to their high tissue N status and comparatively short life cycle. Differences in N utilisation are discussed in the context of species life strategies and mycorrhizal associations.

Subantarctic Macquarie Island - a model ecosystem for studying animal-derived nitrogen sources using N-15 natural abundance

Plants collected from diverse sites on subantarctic Macquarie Island varied by up to 30 parts per thousand in their leaf delta(15)N values. N-15 natural abundance of plants, soils, animal excrement and atmospheric ammonia suggest that the majority of nitrogen utilised by plants growing in the vicinity of animal colonies or burrows is animal-derived. Plants growing near scavengers and animal higher in the food chain had highly enriched delta(15)N values (mean = 12.9 parts per thousand), reflecting the highly enriched signature of these animals' excrement, while plants growing near nesting penguins and albatross, which have an intermediate food chain position, had less enriched delta(15)N values (> 6 parts per thousand). Vegetation in areas affected by rabbits had lower delta(15)N values (mean = 1.2 parts per thousand), while the highly depleted delta(15)N values (below -5 parts per thousand) of plants at upland plateau sites inland of penguin colonies, suggested that a portion of their nitrogen is derived from ammonia (mean N-15 = -10 parts per thousand) lost during the degradation of penguin guano. Vegetation in a remote area had delta(15)N values near -2 parts per thousand. These results contrast with arctic and subarctic studies that attribute large variations in plant N-15 values to nitrogen partitioning in nitrogen-limited environments. Here, plant N-15 reflects the N-15 Of the likely nitrogen sources utilised by plants.

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.

The chemistry of copper patination

The chemistry of copper patination was investigated by two series of experiments. The chemistry of an aqueous copper-sulphate solution was studied at concentrations and pH values near those predicted in an electrolyte on copper exposed to the atmosphere. The electrochemical reactions in an electrolyte in contact with cuprite were investigated in a reaction vessel which used cuprite powder in artificial rainwater to study the electrochemistry of the atmospheric corrosion and patination of copper. Typical sulphate concentrations in rainwater are sufficient to precipitate posnjakite (Cu4SO4(OH)(6)2H(2)O)), a possible precursor to brochantite, within an hour of wetting a cuprite surface. Brochantite (Cu4SO4(OH)(6)), the most commonly found copper salt in natural patinas is responsible for their green appearance. Precipitation of brochantite from the electrolyte resulted from an increase in pH due to the cathodic reduction of oxygen and an increase in cupric ion concentrations by cuprite oxidation. (C) 1998 Elsevier Science Ltd. All rights reserved.

Bandwidth considerations for multilevel converters

Multilevel converters can achieve an overall effective switch frequency multiplication and consequent ripple reduction through the cancellation of the lowest order switch frequency terms. This paper investigates the harmonic content and the frequency response of these multimodulator converters. It is shown that the transfer function of uniformly sampled modulators is a bessel function associated with the inherent sampling process. Naturally sampled modulators have a flat transfer function, but multiple switchings per switch cycle will occur unless the input is slew-rate limited. Lower sideband harmonics of the effective carrier frequency and, in uniform converters, harmonics of the input signal also limit the useful bandwidth. Observations about the effect of the number of converters, their type (naturally or uniformly sampled), and the ratio of modulating frequency and switch frequency are made.

The N-15 natural abundance (delta N-15) of ecosystem samples reflects measures of water availability

We assembled a globally-derived data set for site-averaged foliar delta(15)N, the delta(15)N of whole surface mineral soil and corresponding site factors (mean annual rainfall and temperature, latitude, altitude and soil pH). The delta(15)N of whole soil was related to all of the site variables (including foliar delta(15)N) except altitude and, when regressed on latitude and rainfall, provided the best model of these data, accounting for 49% of the variation in whole soil delta(15)N. As single linear regressions, site-averaged foliar delta(15)N was more strongly related to rainfall than was whole soil delta(15)N. A smaller data set showed similar, negative correlations between whole soil delta(15)N, site-averaged foliar delta(15)N and soil moisture variations during a single growing season. The negative correlation between water availability (measured here by rainfall and temperature) and soil or plant delta(15)N fails at the landscape scale, where wet spots are delta(15)N-enriched relative to their drier surroundings. Here we present global and seasonal data, postulate a proximate mechanism for the overall relationship between water availability and ecosystem delta(15)N and, newly, a mechanism accounting for the highly delta(15)N-depleted values found in the foliage and soils of many wet/cold ecosystems. These hypotheses are complemented by documentation of the present gaps in knowledge, suggesting lines of research which will provide new insights into terrestrial N-cycling. Our conclusions are consistent with those of Austin and Vitousek (1998) that foliar (and soil) delta(15)N appear to be related to the residence time of whole ecosystem N.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

Cytomegalovirus evasion of natural killer cell responses

Natural killer (NK) cells are an important component of the innate cellular immune system. They are particularly important during the early immune responses following virus infection, prior to the induction of cytotoxic T cells (CTL). Unlike CTL, which recognize specific peptides displayed on the surface of cells by class I MHC, NK cells respond to aberrant expression of cell surface molecules, in particular class I MHC, in a non-specific manner. Thus, cells expressing low levels of surface class I MHC are susceptible to recognition by NK cells, with concomitant triggering of cytolytic and cytokine-mediated responses. Many viruses, including the cytomegaloviruses, downregulate cell surface MHC class I: this is likely to provide protection against CTL-mediated clearance of infected cells, but may also render infected cells sensitive to NK-cell attack. This review focuses upon cytomegalovirus-encoded proteins that are believed to promote evasion of NK-cell-mediated immunity. The class I MHC homologues, encoded by all cytomegaloviruses characterised to date, have been implicated as molecular 'decoys', which may mimic the ability of cellular MHC class I to inhibit NK-cell functions. Results from studies in vitro are not uniform, but in general they support the proposal that the class I homologues engage inhibitory receptors from NK cells and other cell types that normally interact with cellular class I. Consistent with this, in vivo studies of murine cytomegalovirus indicate that the class I homologue is required for efficient evasion of NK-cell-mediated clearance. Recently a second murine cytomegalovirus protein, a C-C chemokine homologue, has been implicated as promoting evasion of NK and T-cell-mediated clearance in vivo.

A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.