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.

A variable stepsize implementation for stochastic differential equations

Stochastic differential equations (SDEs) arise from physical systems where the parameters describing the system can only be estimated or are subject to noise. Much work has been done recently on developing higher order Runge-Kutta methods for solving SDEs numerically. Fixed stepsize implementations of numerical methods have limitations when, for example, the SDE being solved is stiff as this forces the stepsize to be very small. This paper presents a completely general variable stepsize implementation of an embedded Runge Kutta pair for solving SDEs numerically; in this implementation, there is no restriction on the value used for the stepsize, and it is demonstrated that the integration remains on the correct Brownian path.

Spatio-temporal evaluation of neck muscle activation during postural perturbations in healthy subjects

The purpose of this study was to examine the spatio-temporal activation of the sternocleidomastoid (SCM) and cervical extensor (CE) muscles with respect to the deltoid muscle onset during rapid voluntary upper limb movement in healthy volunteers. The repeatability and reliability of the spatio-temporal aspects of the myoelectric signals were also examined. Ten subjects performed bilateral and unilateral rapid upper limb flexion, abduction and extension in response to a visual stimulus. EMG onsets and normalised root mean square (nRMS) values were calculated for the SCM and CE muscles. Subjects attended three testing sessions over non-consecutive days allowing the repeatability and reliability of these measures to be assessed. The SCM and CE muscles demonstrated feed-forward activation (activation within 50 ms of deltoid onset) during rapid arm movements in all directions. The sequence and magnitude of neck muscle activation displayed directional specificity, however, the neck flexor and extensor muscles displayed co-activation during all perturbations. EMG onsets demonstrated high repeatability in terms of repeated measure precision (nSEM in the range 1.9-5.7%). This was less evident for the repeatability of nRMS values. The results of this study provide a greater understanding of cervical neuromotor control strategies. During bilateral and unilateral upper limb perturbations, the SCM and CE muscles demonstrate feed-forward co-activation. It seems apparent that feed-forward activation of neck muscles is a mechanism necessary to achieve stability for the visual and vestibular systems, whilst ensuring stabilisation and protection of the cervical spine. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical simulation of stochastic ordinary differential equations in biomathematical modelling

In this work we discuss the effects of white and coloured noise perturbations on the parameters of a mathematical model of bacteriophage infection introduced by Beretta and Kuang in [Math. Biosc. 149 (1998) 57]. We numerically simulate the strong solutions of the resulting systems of stochastic ordinary differential equations (SDEs), with respect to the global error, by means of numerical methods of both Euler-Taylor expansion and stochastic Runge-Kutta type. (C) 2003 IMACS. Published by Elsevier B.V. All rights reserved.

Quantification of lean bodyweight

Background: Lean bodyweight (LBW) has been recommended for scaling drug doses. However, the current methods for predicting LBW are inconsistent at extremes of size and could be misleading with respect to interpreting weight-based regimens. Objective: The objective of the present study was to develop a semi-mechanistic model to predict fat-free mass (FFM) from subject characteristics in a population that includes extremes of size. FFM is considered to closely approximate LBW. There are several reference methods for assessing FFM, whereas there are no reference standards for LBW. Patients and methods: A total of 373 patients (168 male, 205 female) were included in the study. These data arose from two populations. Population A (index dataset) contained anthropometric characteristics, FFM estimated by dual-energy x-ray absorptiometry (DXA - a reference method) and bioelectrical impedance analysis (BIA) data. Population B (test dataset) contained the same anthropometric measures and FFM data as population A, but excluded BIA data. The patients in population A had a wide range of age (18-82 years), bodyweight (40.7-216.5kg) and BMI values (17.1-69.9 kg/m(2)). Patients in population B had BMI values of 18.7-38.4 kg/m(2). A two-stage semi-mechanistic model to predict FFM was developed from the demographics from population A. For stage 1 a model was developed to predict impedance and for stage 2 a model that incorporated predicted impedance was used to predict FFM. These two models were combined to provide an overall model to predict FFM from patient characteristics. The developed model for FFM was externally evaluated by predicting into population B. Results: The semi-mechanistic model to predict impedance incorporated sex, height and bodyweight. The developed model provides a good predictor of impedance for both males and females (r(2) = 0.78, mean error [ME] = 2.30 x 10(-3), root mean square error [RMSE] = 51.56 [approximately 10% of mean]). The final model for FFM incorporated sex, height and bodyweight. The developed model for FFM provided good predictive performance for both males and females (r(2) = 0.93, ME = -0.77, RMSE = 3.33 [approximately 6% of mean]). In addition, the model accurately predicted the FFM of subjects in population B (r(2) = 0.85, ME -0.04, RMSE = 4.39 [approximately 7% of mean]). Conclusions: A semi-mechanistic model has been developed to predict FFM (and therefore LBW) from easily accessible patient characteristics. This model has been prospectively evaluated and shown to have good predictive performance.

Predicting residue-wise contact orders in proteins by support vector regression

Background: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

The HI Parkes All Sky Survey: southern observations, calibration and robust imaging

The acquisition of HI Parkes All Shy Survey (HIPASS) southern sky data commenced at the Australia Telescope National Facility's Parkes 64-m telescope in 1997 February, and was completed in 2000 March. HIPASS is the deepest HI survey yet of the sky south of declination +2 degrees, and is sensitive to emission out to 170 h(75)(-1) Mpc. The characteristic root mean square noise in the survey images is 13.3 mJy. This paper describes the survey observations, which comprise 23 020 eight-degree scans of 9-min duration, and details the techniques used to calibrate and image the data. The processing algorithms are successfully designed to be statistically robust to the presence of interference signals, and are particular to imaging point (or nearly point) sources. Specifically, a major improvement in image quality is obtained by designing a median-gridding algorithm which uses the median estimator in place of the mean estimator.

Quadriceps concentric and eccentric exercise 1: Changes in contractile and electrical activity following eccentric and concentric exercise

The purpose of this study was to determine whether or not losses of strength or endurance following eccentric and concentric exercise are associated with reduced excitation. The effects of eccentric and concentric work on maximal voluntary isometric contraction (MVC) and surface electromyogram (EMG) of the quadriceps were studied in 10 healthy male subjects following bench-stepping for 20 min with a constant leading leg. Prior to stepping and at 0, 0.25, 0.50, 0.75, 1, 3. 24 and 48 h afterwards the subjects performed a 30 s leg extension MVC with each leg during which the isometric force and the root mean square voltage of the EMG were recorded. In the eccentrically exercised muscles (ECC), MVC0-3 (force during the first 3 s of contraction) fen immediately after the bench-stepping exercise to 88 +/- 2% (mean SE) of the pre-exercise value and remained significantly lower than the concentrically exercised muscles (p < 0.05). The muscle weakness in the ECC could not be attributed to central fatigue as surface EMG amplitude at MVC0-3 increased during the recovery period. Muscle weakness after eccentric exercise appears to be due to contractile failure, which is not associated with a reduction in excitation as assessed by surface EMG. Muscular fatigue over 30 s did not change in the two muscle groups after exercise (p = 0.79), indicating that the ECC were weaker but not more fatiguable after exercise.

EMG activity normalization for trunk muscles in subjects with and without back pain

Purpose: The aims of the present study were to examine electromyographic (EMG) activity of six bilateral trunk muscles during maximal contraction in three cardinal planes, and to determine the direction of contraction that gives maximal activation for each muscle. both for healthy subjects and back-pain patients. Methods: Twenty-eight healthy subjects and 15 back-pain patients performed maximum voluntary contractions in three cardinal planes, Surface EMG signals were recorded from rectus abdominis, external oblique, internal oblique, latissimus dorsi, iliocostalis lumborum, and multifidus bilaterally. Root mean square values of the EMG data were calculated to quantify I the amplitude of EMG signals. Results: For both healthy subjects and back-pain patients. one single direction of contraction was found to give the maximum EMG signals for most muscles. Rectus abdominis demonstrated maximal activity in trunk flexion, external oblique in lateral flexion. internal oblique in axial rotation, and multifidus in extension. For the latissimus dorsi and iliocostalis lumborum. maximal activity was demonstrated in more than one cardinal plane. Conclusion: This study has implications for future research involving normalization of muscle activity to maximal levels required in many trunk EMG studies. As the latissimus dorsi and iliocostalis lumborum demonstrate individual differences in the plane that gives maximal activity, these muscles may require testing in more than one plane.

An electromyographic analysis of the deep cervical flexor muscles in performance of craniocervical flexion

Background and Purpose. This study evaluated an electromyographic technique for the measurement of muscle activity of the deep cervical flexor (DCF) muscles. Electromyographic signals were detected from the DCF, sternocleidomastoid (SCM), and anterior scalene (AS) muscles during performance of the craniocervical flexion (CCF) test, which involves performing 5 stages of increasing craniocervical flexion range of motion-the anatomical action of the DCF muscles. Subjects. Ten volunteers without known pathology or impairment participated in this study. Methods. Root-mean-square (RMS) values were calculated for the DCF, SCM, and AS muscles during performance of the CCF test. Myoelectric signals were recorded from the DCF muscles using bipolar electrodes placed over the posterior oropharyngeal wall. Reliability estimates of normalized RMS values were obtained by evaluating intraclass correlation coefficients and the normalized standard error of the mean (SEM). Results. A linear relationship was evident between the amplitude of DCF muscle activity and the incremental stages of the CCF test (F=239.04, df=36, P<.0001). Normalized SEMs in the range 6.7% to 10.3% were obtained for the normalized RMS values for the DCF muscles, providing evidence of reliability for these variables. Discussion and Conclusion. This approach for obtaining a direct measure of the DCF muscles, which differs from those previously used, may be useful for the examination of these muscles in future electromyographic applications.

Patients with chronic neck pain demonstrate altered patterns of muscle activation during performance of a functional upper limb task

Study Design. Cross-sectional study. Objective. This study compared neck muscle activation patterns during and after a repetitive upper limb task between patients with idiopathic neck pain, whiplash-associated disorders, and controls. Summary of Background Data. Previous studies have identified altered motor control of the upper trapezius during functional tasks in patients with neck pain. Whether the cervical flexor muscles demonstrate altered motor control during functional activities is unknown. Methods. Electromyographic activity was recorded from the sternocleidomastoid, anterior scalenes, and upper trapezius muscles. Root mean square electromyographic amplitude was calculated during and on completion of a functional task. Results. A general trend was evident to suggest greatest electromyograph amplitude in the sternocleidomastoid, anterior scalenes, and left upper trapezius muscles for the whiplash-associated disorders group, followed by the idiopathic group, with lowest electromyographic amplitude recorded for the control group. A reverse effect was apparent for the right upper trapezius muscle. The level of perceived disability ( Neck Disability Index score) had a significant effect on the electromyographic amplitude recorded between neck pain patients. Conclusions. Patients with neck pain demonstrated greater activation of accessory neck muscles during a repetitive upper limb task compared to asymptomatic controls. Greater activation of the cervical muscles in patients with neck pain may represent an altered pattern of motor control to compensate for reduced activation of painful muscles. Greater perceived disability among patients with neck pain accounted for the greater electromyographic amplitude of the superficial cervical muscles during performance of the functional task.

Nonlinear dynamics of a satellite with deployable solar panel arrays

The multibody dynamics of a satellite in circular orbit, modeled as a central body with two hinge-connected deployable solar panel arrays, is investigated. Typically, the solar panel arrays are deployed in orbit using preloaded torsional springs at the hinges in a near symmetrical accordion manner, to minimize the shock loads at the hinges. There are five degrees of freedom of the interconnected rigid bodies, composed of coupled attitude motions (pitch, yaw and roll) of the central body plus relative rotations of the solar panel arrays. The dynamical equations of motion of the satellite system are derived using Kane's equations. These are then used to investigate the dynamic behavior of the system during solar panel deployment via the 7-8th-order Runge-Kutta integration algorithms and results are compared with approximate analytical solutions. Chaotic attitude motions of the completely deployed satellite in circular orbit under the influence of the gravity-gradient torques are subsequently investigated analytically using Melnikov's method and confirmed via numerical integration. The Hamiltonian equations in terms of Deprit's variables are used to facilitate the analysis. (C) 2003 Published by Elsevier Ltd.

Impairment in the cervical flexors: a comparison of whiplash and insidious onset neck pain patients

There has been little investigation into whether or not differences exist in the nature of physical impairment associated with neck pain of whiplash and insidious origin. This study examined the neck flexor synergy during performance of the cranio-cervical flexion test, a test targeting the action of the deep neck flexors. Seventy-five volunteer subjects participated in this study and were equally divided between Group 1, asymptomatic control subjects, Group 2, subjects with insidious onset neck pain and Group 3, subjects with neck pain following a whiplash injury. The cranio-cervical flexion test was performed in five progressive stages of increasing cranio-cervical flexion range. Subjects' performance was guided by feedback from a pressure sensor inserted behind the neck which monitored the slight flattening of the cervical lordosis which occurs with the contraction of longus colli. Myoelectric signals (EMG) were detected from the muscles during performance of the test. The results indicated that both the insidious onset neck pain and whiplash groups had higher measures of EMG signal amplitude (normalized root mean square) in the sternocleidomastoid during each stage of the test compared to the control subjects (all P

Patients with neck pain demonstrate reduced electromyographic activity of the deep cervical flexor muscles during performance of the craniocervical flexion test

Study Design. Cross-sectional study. Objective. The present study compared activity of deep and superficial cervical flexor muscles and craniocervical flexion range of motion during a test of craniocervical flexion between 10 patients with chronic neck pain and 10 controls. Summary of Background Data. Individuals with chronic neck pain exhibit reduced performance on a test of craniocervical flexion, and training of this maneuver is effective in management of neck complaints. Although this test is hypothesized to reflect dysfunction of the deep cervical flexor muscles, this has not been tested. Methods. Deep cervical flexor electromyographic activity was recorded with custom electrodes inserted via the nose and fixed by suction to the posterior mucosa of the oropharynx. Surface electrodes were placed over the superficial neck muscles ( sternocleidomastoid and anterior scalene). Root mean square electromyographic amplitude and craniocervical flexion range of motion was measured during five incremental levels of craniocervical flexion in supine. Results. There was a strong linear relation between the electromyographic amplitude of the deep cervical flexor muscles and the incremental stages of the craniocervical flexion test for control and individuals with neck pain ( P = 0.002). However, the amplitude of deep cervical flexor electromyographic activity was less for the group with neck pain than controls, and this difference was significant for the higher increments of the task ( P < 0.05). Although not significant, there was a strong trend for greater sternocleidomastoid and anterior scalene electromyographic activity for the group with neck pain. Conclusions. These data confirm that reduced performance of the craniocervical flexion test is associated with dysfunction of the deep cervical flexor muscles and support the validity of this test for patients with neck pain.

Sampling times for monitoring tacrolimus in stable adult liver transplant recipients

The aim of this study was to determine the most informative sampling time(s) providing a precise prediction of tacrolimus area under the concentration-time curve (AUC). Fifty-four concentration-time profiles of tacrolimus from 31 adult liver transplant recipients were analyzed. Each profile contained 5 tacrolimus whole-blood concentrations (predose and 1, 2, 4, and 6 or 8 hours postdose), measured using liquid chromatography-tandem mass spectrometry. The concentration at 6 hours was interpolated for each profile, and 54 values of AUC(0-6) were calculated using the trapezoidal rule. The best sampling times were then determined using limited sampling strategies and sensitivity analysis. Linear mixed-effects modeling was performed to estimate regression coefficients of equations incorporating each concentration-time point (C0, C1, C2, C4, interpolated C5, and interpolated C6) as a predictor of AUC(0-6). Predictive performance was evaluated by assessment of the mean error (ME) and root mean square error (RMSE). Limited sampling strategy (LSS) equations with C2, C4, and C5 provided similar results for prediction of AUC(0-6) (R-2 = 0.869, 0.844, and 0.832, respectively). These 3 time points were superior to C0 in the prediction of AUC. The ME was similar for all time points; the RMSE was smallest for C2, C4, and C5. The highest sensitivity index was determined to be 4.9 hours postdose at steady state, suggesting that this time point provides the most information about the AUC(0-12). The results from limited sampling strategies and sensitivity analysis supported the use of a single blood sample at 5 hours postdose as a predictor of both AUC(0-6) and AUC(0-12). A jackknife procedure was used to evaluate the predictive performance of the model, and this demonstrated that collecting a sample at 5 hours after dosing could be considered as the optimal sampling time for predicting AUC(0-6).