Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones.

The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented three-dimensional fault zone has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. (C) 2004 Published by Elsevier B.V.

A novel method to prepare mesoporous nano-zirconia

A novel method to prepare mesoporous zirconia was developed. The synthesis was carried out in the presence of PEO surfactants via solid-state reaction. The materials exhibit strong diffraction peak at low 2-theta angle and their nitrogen adsorption/desorption isotherms are typical of IV type with H3 hysteresis loops. The pore structure examined by TEM can be described as wormhole domains. The tetragonal zirconia nanocrystals are uniform in size (around 1.5nm) and their pores center at around 4.6nm. The zirconia nanocrystal growth is mainly via an aggregation mechanism. This study also reveals that the PEO surfactants can interact with the Zr-O-Zr framework to reinforce the thermal stability of zirconia. The ratio of NaOH to ZrOCl2, crystallization and calcination temperature play an important role in the synthesis of mesoporous zirconia.

Quantifying the effect of stress relaxation on excavation

Stress relaxation is relevant to the design of both civil and mining excavations. While many authors refer to the adverse effect of stress relaxation on excavation stability, some present compelling empirical evidence indicating that stress relaxation does not have a significant effect. Establishing clear definitions of stress relaxation was critical to understanding and quantifying stress relaxation of the various types that have been referred to in the literature. This paper defines three types of stress relaxation – partial relaxation, full relaxation and tangential relaxation. Once clear definitions were determined, it became clear that the theoretical arguments and empirical evidence presented by various authors to support their respective cases are not contradictory; rather, the different conclusions can be attributed to different types of stress relaxation. In particular, when the minor principal stress is negative the intermediate principal stress has been identified as significantly affecting jointed rock mass behaviour. The aim of the study was to review and evaluate existing methods of quantifying the effect of stress relaxation around underground excavations and, if necessary, propose a new set of recommendations. An empirical stope stability model, that has been termed the Extended Mathews stability chart, was considered to be the most appropriate method of quantifying the effects of stress relaxation. A new set of guidelines to account for the effect of stress relaxation on excavation stability in the Extended Mathews stability chart has been proposed from a back-analysis of 55 case histories of stress relaxation.

Distortional buckling of I-beams by finite element method

Distortional buckling, unlike the usual lateral-torsional buckling in which the cross-section remains rigid in its own plane, involves distortion of web in the cross-section. This type of buckling typically occurs in beams with slender web and stocky flanges. Most of the published studies assume the web to deform with a cubic shape function. As this assumption may limit the accuracy of the results, a fifth order polynomial is chosen here for the web displacements. The general line-type finite element model used here has two nodes and a maximum of twelve degrees of freedom per node. The model not only can predict the correct coupled mode but also is capable of handling the local buckling of the web.

The evaluation of the spatial and temporal stability of sugarcane farm performance based on yield and commercial cane sugar

In broader catchment scale investigations, there is a need to understand and ultimately exploit the spatial variation of agricultural crops for an improved economic return. In many instances, this spatial variation is temporally unstable and may be different for various crop attributes and crop species. In the Australian sugar industry, the opportunity arose to evaluate the performance of 231 farms in the Tully Mill area in far north Queensland using production information on cane yield (t/ha) and CCS ( a fresh weight measure of sucrose content in the cane) accumulated over a 12-year period. Such an arrangement of data can be expressed as a 3-way array where a farm x attribute x year matrix can be evaluated and interactions considered. Two multivariate techniques, the 3-way mixture method of clustering and the 3-mode principal component analysis, were employed to identify meaningful relationships between farms that performed similarly for both cane yield and CCS. In this context, farm has a spatial component and the aim of this analysis was to determine if systematic patterns in farm performance expressed by cane yield and CCS persisted over time. There was no spatial relationship between cane yield and CCS. However, the analysis revealed that the relationship between farms was remarkably stable from one year to the next for both attributes and there was some spatial aggregation of farm performance in parts of the mill area. This finding is important, since temporally consistent spatial variation may be exploited to improve regional production. Alternatively, the putative causes of the spatial variation may be explored to enhance the understanding of sugarcane production in the wet tropics of Australia.

GCMC-surface area of carbonaceous materials with N-2 and Ar adsorption as an alternative to the classical BET method

In this paper we apply a new method for the determination of surface area of carbonaceous materials, using the local surface excess isotherms obtained from the Grand Canonical Monte Carlo simulation and a concept of area distribution in terms of energy well-depth of solid–fluid interaction. The range of this well-depth considered in our GCMC simulation is from 10 to 100 K, which is wide enough to cover all carbon surfaces that we dealt with (for comparison, the well-depth for perfect graphite surface is about 58 K). Having the set of local surface excess isotherms and the differential area distribution, the overall adsorption isotherm can be obtained in an integral form. Thus, given the experimental data of nitrogen or argon adsorption on a carbon material, the differential area distribution can be obtained from the inversion process, using the regularization method. The total surface area is then obtained as the area of this distribution. We test this approach with a number of data in the literature, and compare our GCMC-surface area with that obtained from the classical BET method. In general, we find that the difference between these two surface areas is about 10%, indicating the need to reliably determine the surface area with a very consistent method. We, therefore, suggest the approach of this paper as an alternative to the BET method because of the long-recognized unrealistic assumptions used in the BET theory. Beside the surface area obtained by this method, it also provides information about the differential area distribution versus the well-depth. This information could be used as a microscopic finger-print of the carbon surface. It is expected that samples prepared from different precursors and different activation conditions will have distinct finger-prints. We illustrate this with Cabot BP120, 280 and 460 samples, and the differential area distributions obtained from the adsorption of argon at 77 K and nitrogen also at 77 K have exactly the same patterns, suggesting the characteristics of this carbon.

Is the method of signal analysis and test selection important for measuring standing balance in subjects with persistent whiplash?

Dizziness and or unsteadiness, associated with episodes of loss of balance, are frequent complaints in those suffering from persistent problems following a whiplash injury. Research has been inconclusive with respect to possible aetiology, discriminative tests and analyses used. The aim of this pilot research was to identify the test conditions and the most appropriate method for the analysis of sway that may differentiate subjects with persistent whiplash associated disorders (WAD) from healthy controls. The six conditions of the Clinical Test for Sensory Interaction in Balance was performed in both comfortable and tandem stance in 20 subjects with persistent WAD compared to 20 control subjects. The analyses were carried out using a traditional method of measurement, total sway distance, to results obtained from the use of wavelet analysis. Subjects with WAD were significantly less able to complete the tandem stance tests on a firm surface than controls. In comfortable stance, using wavelet analysis, significant differences between subjects with WAD and the control group were evident in total energy of the trace for all test conditions apart from eyes open on the firm surface. In contrast, the results of the analysis using total sway distance revealed no significant differences between groups across all six conditions. Wavelet analysis may be more appropriate for detecting disturbances in balance in whiplash subjects because the technique allows separation of the noise from the underlying systematic effect of sway. These findings will be used to direct future studies on the aeitiology of balance disturbances in WAD. (c) 2004 Elsevier B.V. All rights reserved.

Changes in joint stability with muscle contraction measured from transmission of mechanical vibration

A non-invasive in vivo technique was developed to evaluate changes in wrist joint stability properties induced by increased co-activation of the forearm muscles in a gripping task. Mechanical vibration at 45, 50 and 55 Hz was applied to the radial head in ten healthy volunteers. Vibrations of the styloid process of the radius and the distal end of the metacarpal bone of the index finger were measured with triaxial accelerometers. Joint stability properties were quantified by the transfer function gain between accelerations on either side of the wrist-joint. Gain was calculated with the muscles at rest and at five force levels ranging from 5% to 25% of maximum grip force (%MF). During contraction the gain was significantly greater than in control trial (0%MF) for all contractions levels at 45 and 50 Hz and a trend for 15%MF and higher at 55 Hz. Group means of contraction force and gain were significantly correlated at 45 (R-2 = 0.98) and 50 Hz (R-2 = 0.72), but not at 55 Hz (R-2 = 0.10). In conclusion, vibration transmission gain may provide a method to evaluate changes in joint stability properties. (c) 2005 Published by Elsevier Ltd.

A note on the Balanced method

Recently the Balanced method was introduced as a class of quasi-implicit methods for solving stiff stochastic differential equations. We examine asymptotic and mean-square stability for several implementations of the Balanced method and give a generalized result for the mean-square stability region of any Balanced method. We also investigate the optimal implementation of the Balanced method with respect to strong convergence.

The stimulus-response curve and motor unit variability in normal subjects and subjects with amyotrophic lateral sclerosis

The behavior and stability of motor units (MUs) in response to electrical stimulation of different intensities can be assessed with the stimulus-response curve, which is a graphical representation of the size of the compound muscle action potential (CMAP) in relation to stimulus intensity. To examine MU characteristics across the whole stimulus range, the variability of CMAP responses to electrical stimulation, and the differences that occur between normal and disease states, the curve was studied in 11 normal subjects and 16 subjects with amyotrophic lateral sclerosis (ALS). In normal subjects, the curve showed a gradual increase in CMAP size with increasing stimulus intensity, although one or two discrete steps were sometimes observed in the upper half of the curve, indicating the activation of large MUs at higher intensities. In ALS subjects, large discrete steps, due to loss of MUs and collateral sprouting, were frequently present. Variability of the CMAP responses was greater than baseline variability, indicating variability of MU responses, and at certain levels this variability was up to 100 mu Vms. The stimulus-response curve shows differences between normal and ALS subjects and provides information on MU activation and variability throughout the curve.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Measurement and stability of FTY720 in human whole blood by high-performance liquid chromatography-atmospheric pressure chemical ionization-tandem mass spectrometry

We report here a validated method for the quantification of a new immunosuppressant drug FTY720, using HPLC-tandem mass spectrometry. Whole blood samples (500 mu l) were subjected to liquid-liquid extraction, in the presence of an internal standard (Y-32919). Mass spectrometric detection was by selected reaction monitoring with an atmospheric pressure chemical ionization source in positive ionization mode (FTY720: m/z 308.3 -> 255.3). The assay was linear from 0.2 to 25 mu g/l (r(2) > 0.997, n = 5). The inter- and intra-day analytical recovery and imprecision for quality control samples (0.5, 7 and 15 mu g/l) were 95.8-103.2 and < 5.5%, respectively. At the lower limit of quantification (0.2 mu g/l) the interand intra-day analytical recovery was 99.0-102.8% with imprecision of < 7.6% (n = 5). The assay had a mean relative recovery of 100.5 +/- 5.8% (n = 15). Extracted samples were stable for 16 h. IFTY720 quality control samples were stable at room temperature for 16 h at 4 degrees C for at least 8 days and when taken through at least three freeze-thaw cycles. In conclusion, the method described displays analytical performance characteristics that are suitable for pharmacokinetic studies in humans. (c) 2006 Elsevier B.V. All rights reserved.

Eigenvalue sensitivity analysis for probabilistic small signal stability assessment

This paper presents a method to analyze the first order eigenvalue sensitivity with respect to the operating parameters of a power system. The method is based on explicitly expressing the system state matrix into sub-matrices. The eigenvalue sensitivity is calculated based on the explicitly formed system state matrix. The 4th order generator model and 4th order exciter system model are used to form the system state matrix. A case study using New England 10-machine 39-bus system is provided to demonstrate the effectiveness of the proposed method. This method can be applied into large scale power system eigenvalue sensitivity with respect to operating parameters.

Optimum switching point assessment of parallel braking resistor and serial current limiter for power system stability enhancement

A new control algorithm using parallel braking resistor (BR) and serial fault current limiter (FCL) for power system transient stability enhancement is presented in this paper. The proposed control algorithm can prevent transient instability during first swing by immediately taking away the transient energy gained in faulted period. It can also reduce generator oscillation time and efficiently make system back to the post-fault equilibrium. The algorithm is based on a new system energy function based method to choose optimal switching point. The parallel BR and serial FCL resistor can be switched at the calculated optimal point to get the best control result. This method allows optimum dissipation of the transient energy caused by disturbance so to make system back to equilibrium in minimum time. Case studies are given to verify the efficiency and effectiveness of this new control algorithm.

Monitoring the condition of savanna riparian zones in Northern Australia using high spatial resolution imagery and a field based rapid assessment method

Government agencies responsible for riparian environments are assessing the combined utility of field survey and remote sensing for mapping and monitoring indicators of riparian zone condition. The objective of this work was to compare the Tropical Rapid Appraisal of Riparian Condition (TRARC) method to a satellite image based approach. TRARC was developed for rapid assessment of the environmental condition of savanna riparian zones. The comparison assessed mapping accuracy, representativeness of TRARC assessment, cost-effectiveness, and suitability for multi-temporal analysis. Two multi-spectral QuickBird images captured in 2004 and 2005 and coincident field data covering sections of the Daly River in the Northern Territory, Australia were used in this work. Both field and image data were processed to map riparian health indicators (RHIs) including percentage canopy cover, organic litter, canopy continuity, stream bank stability, and extent of tree clearing. Spectral vegetation indices, image segmentation and supervised classification were used to produce RHI maps. QuickBird image data were used to examine if the spatial distribution of TRARC transects provided a representative sample of ground based RHI measurements. Results showed that TRARC transects were required to cover at least 3% of the study area to obtain a representative sample. The mapping accuracy and costs of the image based approach were compared to those of the ground based TRARC approach. Results proved that TRARC was more cost-effective at smaller scales (1-100km), while image based assessment becomes more feasible at regional scales (100-1000km). Finally, the ability to use both the image and field based approaches for multi-temporal analysis of RHIs was assessed. Change detection analysis demonstrated that image data can provide detailed information on gradual change, while the TRARC method was only able to identify more gross scale changes. In conclusion, results from both methods were considered to complement each other if used at appropriate spatial scales.