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.

Bioelectrical impedance analysis for the estimation of body composition in rats

Bioelectrical impedance analysis (BIA) was used to assess body composition in rats fed on either standard laboratory diet or on a high-fat diet designed to induce obesity. Bioelectrical impedance analysis predictions of total body water and thus fat-free mass (FFM) for the group mean values were generally within 5% of the measured values by tritiated water ((H2O)-H-3) dilution. The limits of agreement for the procedure were, however, large, approximately +/-25%, limiting the applicability of the technique for measurement of body composition in individual animals.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - II. Testing simulation capabilities of wheat growth

To simulate cropping systems, crop models must not only give reliable predictions of yield across a wide range of environmental conditions, they must also quantify water and nutrient use well, so that the status of the soil at maturity is a good representation of the starting conditions for the next cropping sequence. To assess the suitability for this task a range of crop models, currently used in Australia, were tested. The models differed in their design objectives, complexity and structure and were (i) tested on diverse, independent data sets from a wide range of environments and (ii) model components were further evaluated with one detailed data set from a semi-arid environment. All models were coded into the cropping systems shell APSIM, which provides a common soil water and nitrogen balance. Crop development was input, thus differences between simulations were caused entirely by difference in simulating crop growth. Under nitrogen non-limiting conditions between 73 and 85% of the observed kernel yield variation across environments was explained by the models. This ranged from 51 to 77% under varying nitrogen supply. Water and nitrogen effects on leaf area index were predicted poorly by all models resulting in erroneous predictions of dry matter accumulation and water use. When measured light interception was used as input, most models improved in their prediction of dry matter and yield. This test highlighted a range of compensating errors in all modelling approaches. Time course and final amount of water extraction was simulated well by two models, while others left up to 25% of potentially available soil water in the profile. Kernel nitrogen percentage was predicted poorly by all models due to its sensitivity to small dry matter changes. Yield and dry matter could be estimated adequately for a range of environmental conditions using the general concepts of radiation use efficiency and transpiration efficiency. However, leaf area and kernel nitrogen dynamics need to be improved to achieve better estimates of water and nitrogen use if such models are to be use to evaluate cropping systems. (C) 1998 Elsevier Science B.V.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.

Two-body and three-body abrasion: A critical discussion

It is argued that the common classification of abrasive wear into 'two-body abrasion' and 'three-body abrasion' is seriously flawed. No definitions have been agreed upon for these terms, and indeed there are two quite different interpretations, the implications of which are mutually inconsistent. In the dominant interpretation, the primary thrust of the two-body/three-body concept is to describe whether the abrasive particles are constrained (two-body) or free to roll (three-body). In this view, two-body abrasion is generally much more severe than three-body. The alternative interpretation emphasises the presence (three-body) or absence (two-body) of a rigid counterface backing the abrasive. In this view, three-body abrasion is equated to high-stress (or grinding) abrasion and is generally more severe than two-body (low-stress) abrasion. This paper recommends that the 'two-body/three-body' terminology be abandoned, to be replaced by an alternative classification scheme based directly upon the manifest severity of wear. (C) 1998 Elsevier Science S.A.

Bioimpedance spectrometry in the determination of body water compartments: Accuracy and clinical significance

Bioelectrical impedance analysis (BIA) offers the potential for a simple, portable and relatively inexpensive technique for the in vivo measurement of total body water (TBW). The potential of BIA as a technique of body composition analysis is even greater when one considers that body water can be used as a surrogate measure of lean body mass. However, BIA has not found universal acceptance even with the introduction of multi-frequency BIA (MFBIA) which, potentially, may improve the predictive accuracy of the measurement. There are a number of reasons for this lack of acceptance, although perhaps the major reason is that no single algorithm has been developed which can be applied to all subject groups. This may be due, in part, to the commonly used wrist-to-ankle protocol which is not indicated by the basic theory of bioimpedance, where the body is considered as five interconnecting cylinders. Several workers have suggested the use of segmental BIA measurements to provide a protocol more in keeping with basic theory. However, there are other difficulties associated with the application of BIA, such as effects of hydration and ion status, posture and fluid distribution. A further putative advantage of MFBIA is the independent assessment not only of TBW but also of the extracellular fluid volume (ECW), hence heralding the possibility of,being able to assess the fluid distribution between these compartments. Results of studies in this area have been, to date, mixed. Whereas strong relationships of impedance values at low frequencies with ECW, and at high frequencies with TBW, have been reported, changes in impedance are not always well correlated with changes in the size of the fluid compartments (assessed by alternative and more direct means) in pathological conditions. Furthermore, the theoretical advantages of Cole-Cole modelling over selected frequency prediction have not always been apparent. This review will consider the principles, methodology and applications of BIA. The principles and methodology will,be considered in relation to the basic theory of BIA and difficulties experienced in its application. The relative merits of single and multiple frequency BIA will be addressed, with particular attention to the latter's role in the assessment of compartmental fluid volumes. (C) 1998 Elsevier Science Ltd. All rights reserved.

Velocimetry and thermometry of supersonic flow around a cylindrical body

The supersonic flow around a cylindrical body has been studied using two optical techniques. For both sets of measurements, the cylinder was mounted from the side of the tunnel, allowing investigation of the bow shock region as well as in the wake. A new technique, laser-enhanced ionization flow tagging, was used for streamwise velocity determinations behind the body. From these measurements, it was found that the downstream velocity outside the wake was (1.90 +/- 0.06) km/s, whereas inside the wake the velocity was about 0-500 m/s in the upstream direction. Planar laser induced fluorescence of nitric oxide was employed for temperature determinations. It was established that the freestream temperature was (2120 +/- 100) K, decreasing to around (1550 +/- 400) K in the wake.

Erosion/productivity modelling of maize farming in the Philippine uplands part II: Simulation of alternative farming methods

A version of the Agricultural Production Systems Simulator (APSIM) capable of simulating the key agronomic aspects of intercropping maize between legume shrub hedgerows was described and parameterised in the first paper of this series (Nelson et al., this issue). In this paper, APSIM is used to simulate maize yields and soil erosion from traditional open-field farming and hedgerow intercropping in the Philippine uplands. Two variants of open-field farming were simulated using APSIM, continuous and fallow, for comparison with intercropping maize between leguminous shrub hedgerows. Continuous open-field maize farming was predicted to be unsustainable in the long term, while fallow open-field farming was predicted to slow productivity decline by spreading the effect of erosion over a larger cropping area. Hedgerow intercropping was predicted to reduce erosion by maintaining soil surface cover during periods of intense rainfall, contributing to sustainable production of maize in the long term. In the third paper in this series, Nelson et al. (this issue) use cost-benefit analysis to compare the economic viability of hedgerow intercropping relative to traditional open-field farming of maize in relatively inaccessible upland areas. (C) 1998 Elsevier Science Ltd. All rights reserved.

Flush air data system calibration using numerical simulation

The use of computational fluid dynamics simulations for calibrating a flush air data system is described, In particular, the flush air data system of the HYFLEX hypersonic vehicle is used as a case study. The HYFLEX air data system consists of nine pressure ports located flush with the vehicle nose surface, connected to onboard pressure transducers, After appropriate processing, surface pressure measurements can he converted into useful air data parameters. The processing algorithm requires an accurate pressure model, which relates air data parameters to the measured pressures. In the past, such pressure models have been calibrated using combinations of flight data, ground-based experimental results, and numerical simulation. We perform a calibration of the HYFLEX flush air data system using computational fluid dynamics simulations exclusively, The simulations are used to build an empirical pressure model that accurately describes the HYFLEX nose pressure distribution ol cr a range of flight conditions. We believe that computational fluid dynamics provides a quick and inexpensive way to calibrate the air data system and is applicable to a broad range of flight conditions, When tested with HYFLEX flight data, the calibrated system is found to work well. It predicts vehicle angle of attack and angle of sideslip to accuracy levels that generally satisfy flight control requirements. Dynamic pressure is predicted to within the resolution of the onboard inertial measurement unit. We find that wind-tunnel experiments and flight data are not necessary to accurately calibrate the HYFLEX flush air data system for hypersonic flight.

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.

Growth hormone binding protein correlates strongly with leptin and percentage body fat in GH-deficient adults, is increased by GH replacement but does not predict IGF-I response

GH-binding protein (GHBP) corresponds to the extracellular domain of the GH receptor (GHR) and has been shown to be closely related to body fat. This study aimed to examine the inter-relationship between GHBP, leptin and body fat, and to test the hypothesis that GHBP is modified by GH replacement in GH-deficient adults and predicts IGF-I response. Twenty adults, mean age 47 years (range 20-69) with proven GH deficiency were randomly allocated to either GH (up to 0.25 U/kg/week in daily doses) or placebo for 3 months before cross-over to the opposite treatment. Plasma GHBP and leptin were measured at baseline and 2, 4, 8 and 12 weeks after each treatment. Whole body composition was measured at baseline by dual-energy X-ray absorptiometry (DEXA). There was a strong correlation between baseline leptin and GHBP (r = 0.88, P < 0.0001) and between baseline GHBP and percentage body fat, (r = 0.83, P < 0.0001). Mean GHBP levels were higher on GH compared with placebo, 1.53 +/- 0.28 vs 1.41 +/- 0.25 nM, P = 0.049. There was no correlation between baseline IGF-I and GHBP (r = -0.049, P = 0.84), and GHBP did not predict IGF-I response to GH replacement. The close inter-relationship between GHBP, leptin and body fat suggests a possible role for GHBP in the regulation of body composition. GHBP is increased by GH replacement in GH-deficient adults, but does not predict biochemical response to GH replacement. (C) 1999 Churchill Livingstone.

New techniques in nutritional assessment: body composition methods

New techniques in air-displacement plethysmography seem to have overcome many of the previous problems of poor reproducibility and validity. These have made body-density measurements available to a larger range of individuals, including children, elderly and sick patients who often have difficulties in being submerged underwater in hydrodensitometry systems. The BOD POD air-displacement system (BOD POD body composition system; Life Measurement Instruments, Concord, CA, USA) is more precise than hydrodensitometry, is simple and rapid to operate (approximately 1 min measurements) and the results agree closely with those of hydrodensitometry (e.g. +/-3.4% for estimation of body fat). Body line scanners employing the principles of three-dimensional photography are potentially able to measure the surface area and volume of the body and its segments even more rapidly (approximately 10 s), but the validity of the measurements needs to be established. Advances in i.r. spectroscopy and mathematical modelling for calculating the area under the curve have improved precision for measuring enrichment of (H2O)-H-2 in studies of water dilution (CV 0.1-0.9% within the range of 400-1000 mu l/l) in saliva, plasma and urine. The technique is rapid and compares closely with mass spectrometry (bias 1 (SD 2) %). Advances in bedside bioelectrical-impedance techniques are making possible potential measurements of skinfold thicknesses and limb muscle mass electronically. Preliminary results suggest that the electronic method is more reproducible (intra-and inter-individual reproducibility for measuring skinfold thicknesses) and associated with less bias (+ 12%), than anthropometry (+ 40%). In addition to these selected examples, the 'mobility' or transfer of reference methods between centres has made the distinction between reference and bedside or field techniques less distinct than in the past.