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.

Delayed postural contraction of transversus abdominis in low back pain associated with movement of the lower limb

The temporal parameters of the response of the trunk muscles associated with movement of the lower limb were investigated in people with and without low back pain (LBP). The weight shift component of the task was completed voluntarily prior to a stimulus to move to allow investigation of the movement component of the response. In the control subjects the onset of electromyographic (EMG) activity of all trunk muscles preceded that of the muscle responsible for limb movement, thus contributing to the feed forward postural response. The EMG onset of transversus abdominis was delayed in the LBP subjects with movement in each direction, while the EMG onsets of rectus abdominis, erector spinae, and oblique abdominal muscles were delayed with specific movement directions. This result provides evidence of a change in the postural control of the trunk in people with LBP.

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.

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.

Preparatory trunk motion accompanies rapid upper limb movement

Evaluation of trunk movements, trunk muscle activation, intra-abdominal pressure and displacement of centres of pressure and mass was undertaken to determine whether trunk orientation is a controlled variable prior to and during rapid bilateral movement of the upper limbs. Standing subjects performed rapid bilateral symmetrical upper limb movements in three directions (flexion, abduction and extension). The results indicated a small (0.4-3.3 degrees) but consistent initial angular displacement between the segments of the trunk in a direction opposite to that produced by the reactive moments resulting from limb movement. Phasic activation of superficial trunk muscles was consistent with this pattern of preparatory motion and with the direction of motion of the centre of mass. In contrast, activation of the deep abdominal muscles was independent of the direction of limb motion, suggesting a non-direction specific contribution to spinal stability. The results support the opinion that feedforward postural responses result in trunk movements, and that orientation of the trunk and centre of mass are both controlled variables in relation to rapid limb movements.

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.

Relationships between movement initiation times and movement-related cortical potentials in Parkinson's disease

Movement-related cortical potentials recorded from the scalp reveal increasing cortical activity occurring prior to voluntary movement. Studies of set-related cortical activity recorded from single neurones within premotor and supplementary motor areas in monkeys suggest that such premovement activity may act to prime activity of appropriate motor units in readiness to move, thereby facilitating the movement response. Such a role of early stage premovement activity in movement-related cortical potentials was investigated by examining the relationship between premovement cortical activity and movement initiation or reaction times. Parkinson's disease and control subjects performed a simple button-pressing reaction time task and individual movement-related potentials were averaged for responses with short compared with long reaction times. For Parkinson's disease subjects but not for the control subjects, early stage premovement cortical activity was significantly increased in amplitude for faster reaction times, indicating that there is indeed a relationship between premovement cortical activity amplitude and movement initiation or reaction times. In support of studies of set-related cortical activity in monkeys, it is therefore suggested that early stage premovement activity reflects the priming of appropriate motor units of primary motor cortex, thereby reducing movement initiation or reaction times. (C) 1999 Elsevier Science B.V. All rights reserved.

The moderating role of task characteristics in determining responses to a stressful work simulation

Two experimental studies were conducted to examine whether the stress-buffering effects of behavioral control on work task responses varied as a function of procedural information. Study 1 manipulated low and high levels of task demands, behavioral control, and procedural information for 128 introductory psychology students completing an in-basket activity. ANOVA procedures revealed a significant three-way interaction among these variables in the prediction of subjective task performance and task satisfaction. It was found that procedural information buffered the negative effects of task demands on ratings of performance and satisfaction only under conditions of low behavioral control. This pattern of results suggests that procedural information may have a compensatory effect when the work environment is characterized by a combination of high task demands and low behavioral control. Study 2 (N = 256) utilized simple and complex versions of the in-basket activity to examine the extent to which the interactive relationship among task demands, behavioral control, and procedural information varied as a function of task complexity. There was further support for the stress-buffering role of procedural information on work task responses under conditions of low behavioral control. This effect was, however, only present when the in-basket activity was characterized by high task complexity, suggesting that the interactive relationship among these variables may depend on the type of tasks performed at work. Copyright (C) 1999 John Wiley & Sons, Ltd.

Movement-related potentials in Parkinson's disease: External cues and attentional strategies

Hypokinetic movement can be greatly improved in Parkinson's disease patients by the provision of external cues to guide movement. It has recently been reported, however, that movement performance in parkinsonian patients can be similarly improved in the absence of external cues by using attentional strategies, whereby patients are instructed to consciously attend to particular aspects of the movement which would normally be controlled automatically. To study the neurophysiological basis of such improvements in performance associated with the use of attentional strategies, movement-related cortical potentials were examined in Parkinson's disease and control subjects using a reaction time paradigm. One group of subjects were explicitly instructed to concentrate on internally timed responses to anticipate the presentation of a predictably timed go signal. Other subjects were given no such instruction regarding attentional strategies. Early-stage premovement activity of movement-related potentials was significantly increased in amplitude and reaction times were significantly faster for Parkinson's disease subjects when instructed to direct their attention toward internally generating responses rather than relying on external cues. It is therefore suggested that the use of attentional strategies may allow movement to be mediated by less automatic and more conscious attentional motor control processes which may be less impaired by basal ganglia dysfunction, and thereby improve movement performance in Parkinson's disease.

Season-specific directional movement in migratory Australian butterflies

Large numbers of adults of certain species of butterfly flying in an apparently 'purposeful' manner are often noted by entomologists and the general public. Occasionally, these are recorded in the literature. Using these records we summarise information regarding the direction of movement in Australian butterflies and test whether there are consistent patterns that could account for known seasonal shifts in geographical range. The data were analysed using contingency tables and directionality statistics. Vanessa itea, Vanessa kershawi, Danaus plexippus, Danaus chrysippus and Badamia exclamationis flew predominately south in the spring-summer and north in the autumn-winter. Tirumala hamata has a strong southern component to its flight in spring but, as in Euploea core, appears non-directional in the autumn. For many supposedly known migratory species, the number of literature records are few, particularly in one season (mainly autumn). Thus, for Appias paulina, four of seven records were south in the spring-summer, as were six of nine records for Catopsilia pomona, and three of five for Zizina labradus. For Belenois java, flight records were only available for the spring and these showed geographical differences; predominantly north-west in northern Australia (Queensland) and south-west in southern Australia (Victoria, New South Wales). There were too few records for Papilio demoleus in the literature (four only) to draw any conclusions. Major exceptions to the seasonal trend of south in the spring and north in the autumn were Junonia villida, which showed a predominant north-westward direction in both seasons, and Eurema smilax, with a predominant southern or western flight in both seasons. We discuss these species specific trends in migration direction in relation to seasonal shifts in suitable habitat conditions, possible cues used in orientation and in timing changes in direction.