Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project

To foster ongoing international cooperation beyond ACES (APEC Cooperation for Earthquake Simulation) on the simulation of solid earth phenomena, agreement was reached to work towards establishment of a frontier international research institute for simulating the solid earth: iSERVO = International Solid Earth Research Virtual Observatory institute (http://www.iservo.edu.au). This paper outlines a key Australian contribution towards the iSERVO institute seed project, this is the construction of: (1) a typical intraplate fault system model using practical fault system data of South Australia (i.e., SA interacting fault model), which includes data management and editing, geometrical modeling and mesh generation; and (2) a finite-element based software tool, which is built on our long-term and ongoing effort to develop the R-minimum strategy based finite-element computational algorithm and software tool for modelling three-dimensional nonlinear frictional contact behavior between multiple deformable bodies with the arbitrarily-shaped contact element strategy. A numerical simulation of the SA fault system is carried out using this software tool to demonstrate its capability and our efforts towards seeding the iSERVO Institute.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Bifurcation analysis of a dynamic power system with SVC devices

Bifurcation analysis is a very useful tool for power system stability assessment. In this paper, detailed investigation of power system bifurcation behaviour is presented. One and two parameter bifurcation analysis are conducted on a 3-bus power system. We also examined the impact of FACTS devices on power system stability through Hopf bifurcation analysis by taking static Var compensator (SVC) as an example. A simplified first-order model of the SVC device is included in the 3-bus sample system. Real and reactive powers are used as bifurcation parameter in the analysis to compare the system oscillatory properties with and without SVC. The simulation results indicate that the linearized system model with SVC enlarge the voltage stability boundary by moving Hopf bifurcation point to higher level of loading conditions. The installation of SVC increases the dynamic stability range of the system, however complicates the Hopf bifurcation behavior of the system

Tyre-force estimation by Kalman inverse filtering: applications to off-highway mining trucks

Kalman inverse filtering is used to develop a methodology for real-time estimation of forces acting at the interface between tyre and road on large off-highway mining trucks. The system model formulated is capable of estimating the three components of tyre-force at each wheel of the truck using a practical set of measurements and inputs. Good tracking is obtained by the estimated tyre-forces when compared with those simulated by an ADAMS virtual-truck model. A sensitivity analysis determines the susceptibility of the tyre-force estimates to uncertainties in the truck's parameters.

In vivo demonstration of load-induced fluid flow in the rat tibia and its potential implications for processes associated with functional adaptation

Load-induced extravascular fluid flow has been postulated to play a role in mechanotransduction of physiological loads at the cellular level. Furthermore, the displaced fluid serves as a carrier for metabolites, nutrients, mineral precursors and osteotropic agents important for cellular activity. We hypothesise that load-induced fluid flow enhances the transport of these key substances, thus helping to regulate cellular activity associated with processes of functional adaptation and remodelling. To test this hypothesis, molecular tracer methods developed previously by our group were applied in vivo to observe and quantify the effects of load-induced fluid flow under four-point-bending loads. Preterminal tracer transport studies were carried out on 24 skeletally mature Sprague Dawley rats. Mechanical loading enhanced the transport of both small- and larger-molecular-mass tracers within the bony tissue of the tibial mid-diaphysis. Mechanical loading showed a highly significant effect on the number of periosteocytic spaces exhibiting tracer within the cross section of each bone. For all loading rates studied, the concentration of Procion Red tracer was consistently higher in the tibia subjected to pure bending loads than in the unloaded, contralateral tibia, Furthermore, the enhancement of transport was highly site-specific. In bones subjected to pure bending loads, a greater number of periosteocytic spaces exhibited the presence of tracer in the tension band of the cross section than in the compression band; this may reflect the higher strains induced in the tension band compared with the compression band within the mid-diaphysis of the rat tibia. Regardless of loading mode, the mean difference between the loaded side and the unloaded contralateral control side decreased with increasing loading frequency. Whether this reflects the length of exposure to the tracer or specific frequency effects cannot be determined by this set of experiments. These in vivo experimental results corroborate those of previous ex vivo and in vitro studies, Strain-related differences in tracer distribution provide support for the hypothesis that load-induced fluid flow plays a regulatory role in processes associated with functional adaptation.

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.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

A quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-SiO2 system

A structurally based viscosity model for fully liquid silicate slags has been proposed and applied to the Al2O3-CaO-'FeO'-SiO2 system at metallic iron saturation. The model links the slag viscosity to the internal structure of melts through the concentrations of various anion/cation structural units (SUs). The concentrations of structural units are equivalent to the second nearest neighbor bond concentrations calculated by the quasi-chemical thermodynamic model. This viscosity model describes experimental data over the entire temperature and composition range within the Al2O3-CaO-'FeO'-SiO2 system at metallic iron saturation and can be extended to other industrial slag systems.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 1. Description of the model and its application to the MgO, MgO-SiO2, Al2O3-MgO and CaO-MgO sub-systems

A structurally-based quasi-chemical viscosity model for fully liquid slags in the Al2O3 CaO-'FeO'-MgO-SiO2 system has been developed. The model links the slag viscosities to the internal structures of the melts through the concentrations of various Si0.5O, Me2/nn+O and Me1/nn+Si0.25O viscous flow structural units. The concentrations of these structural units are derived from a quasi-chemical thermodynamic model of the system. The model described in this series of papers enables the viscosities of liquid slags to be predicted within experimental uncertainties over the whole range of temperatures and compositions in the Al2O3 CaOMgO-SiO2 system.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 2. A review of the experimental data and the model predictions for the Al2O3-CaO-MgO, CaO-MgO-SiO2 and Al2O3-MgO-SiO2 systems

A structurally-based quasi-chemical viscosity model has been developed for the Al2O3 CaO-'FeO'-MgO-SiO2 system. The model links the slag viscosity to the internal structure of melts through the concentrations of various anion/cation Si0.5O, Me2/nn+O and Me1/nn+Si0.25O viscous flow structural units. The concentrations of structural units are derived from the quasi-chemical thermodynamic model. The focus of the work described in the present paper is the analysis of experimental data and the viscosity models for fully liquid slags in the Al2O3-CaO-MgO, Al2O3 MgO-SiO2 and CaO-MgO-SiO2 systems.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 3. Summary of the model predictions for the Al2O3-CaO-MgO-SiO2 system and its sub-systems

A structurally-based quasi-chemical viscosity model for fully liquid slags in the Al2O3 CaO-'FeO'-MgOSiO2 system has been developed. The focus of the work described in the present paper is the analysis of the experimental data and viscosity models in the quaternary system Al2O3 CaO-MgO-SiO2 and its subsystems. A review of the experimental data, viscometry methods used and viscosity models available in the Al2O3 CaO-MgO-SiO2 and its sub-systems is reported. The quasi-chemical viscosity model is shown to provide good agreement between experimental data and predictions over the whole compositional range.

Prediction of vapor-moisture equilibriums for a wood-moisture system using a modified UNIQUAC model

A modified UNIQUAC model has been extended to describe and predict the equilibrium relative humidity and moisture content for wood. The method is validated over a range of moisture content from oven-dried state to fiber saturation point, and over a temperature range of 20-70 degrees C. Adjustable parameters and binary interaction parameters of the UNIQUAC model were estimated from experimental data for Caribbean pine and Hoop pine as well as data available in the literature. The two group-interaction parameters for the wood-moisture system were consistent with using function group contributions for H2O, -OH and -CHO. The result reconfirms that the main contributors to water adsorption in cell walls are the hydroxyl groups of the carbohydrates in cellulose and hemicelluloses. This provides some physical insight into the intermolecular force and energy between bound water and the wood material. (c) 2006 Elsevier Ltd. All rights reserved.

Polyethism and comparability of termite choice assays in a model system using Microcerotermes turneri (Termitidae : Termitinae): Implications for standardised testing techniques

Traditional measures of termite food preference assess consequences of foraging behavior such as wood consumption, aggregation and/or termite survivorship. Although studies have been done to investigate the specifics of foraging behavior this is not generally integrated into choice assay experiments. Here choice assays were conducted with small isolated (orphaned) groups of workers and compared with choice assays involving foragers from whole nests (non-orphaned) in the laboratory. Aggregation to two different wood types was used as a measure of preference. Specific worker caste and instars participating in initial exploration were compared between assay methods, with samples of termites taken from nest carton material and sites where termites were feeding. Aggregation results differ between choice assay techniques. Castes and instars responsible for initial exploration, as determined in whole nest trials, were not commonly found exploring in isolated group trials, nor were they numerous in termites taken from active feeding sites. Consequently the use of small groups of M. turneri worker termites extracted from active feeding sites may not be appropriate for use in choice assays.