Modification of the in vivo four-point loading model for studying mechanically induced bone adaptation

We modified the noninvasive, in vivo technique for strain application in the tibiae of rats (Turner et al,, Bone 12:73-79, 1991), The original model applies four-point bending to right tibiae via an open-loop, stepper-motor-driven spring linkage, Depending on the magnitude of applied load, the model produces new bone formation at periosteal (Ps) or endocortical surfaces (Ec.S). Due to the spring linkage, however, the range of frequencies at which loads can be applied is limited. The modified system replaces this design with an electromagnetic vibrator. A load transducer in series with the loading points allows calibration, the loaders' position to be adjusted, and cyclic loading completed under load central as a closed servo-loop. Two experiments were conducted to validate the modified system: (1) a strain gauge was applied to the lateral surface of the right tibia of 5 adult female rats and strains measured at applied loads from 10 to 60 N; and (2) the bone formation response was determined in 28 adult female Sprague-Dawley rats. Loading was applied as a haversine wave with a frequency of 2 Hz for 18 sec, every second day for 10 days. Peak bending loads mere applied at 33, 40, 52, and 64 N, and a sham-loading group tr as included at 64 N, Strains in the tibiae were linear between 10 and 60 N, and the average peak strain at the Ps.S at 60 N was 2664 +/- 250 microstrain, consistent with the results of Turner's group. Lamellar bone formation was stimulated at the Ec.S by applied bending, but not by sham loading. Bending strains above a loading threshold of 40 N increased Ec Lamellar hone formation rate, bone forming surface, and mineral apposition rate with a dose response similar to that reported by Turner et al, (J Bone Miner Res 9:87-97, 1994). We conclude that the modified loading system offers precision for applied loads of between 0 and 70 N, versatility in the selection of loading rates up to 20 Hz, and a reproducible bone formation response in the rat tibia, Adjustment of the loader also enables study of mechanical usage in murine tibia, an advantage with respect to the increasing variety of transgenic strains available in bone and mineral research. (Bone 23:307-310; 1998) (C) 1998 by Elsevier Science Inc. All rights reserved.

PVT behaviour of fluids with potential models optimised for phase equilibria

NPT and NVT Monte Carlo simulations are applied to models for methane and water to predict the PVT behaviour of these fluids over a wide range of temperatures and pressures. The potential models examined in this paper have previously been presented in the literature with their specific parameters optimised to fit phase coexistence data. The exponential-6 potential for methane gives generally good prediction of PVT behaviour over the full range of temperature and pressures studied with the only significant deviation from experimental data seen at high temperatures and pressures. The NSPCE water model shows very poor prediction of PVT behaviour, particularly at dense conditions. To improve this. the charge separation in the NSPCE model is varied with density. Improvements for vapour and liquid phase PVT predictions are achieved with this variation. No improvement was found in the prediction of the oxygen-oxygen radial distribution by varying charge separation under dense phase conditions. (C) 2004 Elsevier B.V. All rights reserved.

Long-term loading inhibits ERK1/2 phosphorylation and increases FGFR3 expression in MC3T3-E1 osteoblast cells

Bone tissue homeostasis relies upon the ability of cells to detect and interpret extracellular signals that direct changes in tissue architecture. This study utilized a four-point bending model to create both fluid shear and strain forces (loading) during the time-dependent progression of MC3T3-E1 preosteoblasts along the osteogenic lineage. Loading was shown to increase cell number, alkaline phosphatase (ALP) activity, collagen synthesis, and the mRNA expression levels of Runx2, osteocalcin (OC), osteopontin, and cyclo-oxygenase-2. However, mineralization in these cultures was inhibited, despite an increase in calcium accumulation, suggesting that loading may inhibit mineralization in order to increase matrix deposition. Loading also increased fibroblast growth factor receptor-3 (FGFR3) expression coincident with an inhibition of FGFR1, FGFR4, FGF1, and extracellular signal-related kinase (ERK)1/2 phosphorylation. To examine whether these loading-induced changes in cell phenotype and FGFR expression could be attributed to the inhibition of ERK1/2 phosphorylation, cells were grown for 25 days in the presence of the MEK1/2 inhibitor, U0126. Significant increases in the expression of FGFR3, ALP, and OC were observed, as well as the inhibition of FGFR1, FGFR4, and FGF1. However, U0126 also increased matrix mineralization, demonstrating that inhibition of ERK1/2 phosphorylation cannot fully account for the changes observed in response to loading. in conclusion, this study demonstrates that preosteoblasts are mechanoresponsive, and that long-term loading, whilst increasing proliferation and differentiation of preosteoblasts, inhibits matrix mineralization. In addition, the increase in FGFR3 expression suggests that it may have a role in osteoblast differentiation.

3D statistical shape models to embed spatial relationship information

This paper presents the creation of 3D statistical shape models of the knee bones and their use to embed information into a segmentation system for MRIs of the knee. We propose utilising the strong spatial relationship between the cartilages and the bones in the knee by embedding this information into the created models. This information can then be used to automate the initialisation of segmentation algorithms for the cartilages. The approach used to automatically generate the 3D statistical shape models of the bones is based on the point distribution model optimisation framework of Davies. Our implementation of this scheme uses a parameterized surface extraction algorithm, which is used as the basis for the optimisation scheme that automatically creates the 3D statistical shape models. The current approach is illustrated by generating 3D statistical shape models of the patella, tibia and femoral bones from a segmented database of the knee. The use of these models to embed spatial relationship information to aid in the automation of segmentation algorithms for the cartilages is then illustrated.

Automatic segmentation of the knee bones using 3D active shape models

This paper presents an automated segmentation approach for MR images of the knee bones. The bones are the first stage of a segmentation system for the knee, primarily aimed at the automated segmentation of the cartilages. The segmentation is performed using 3D active shape models (ASM), which are initialized using an affine registration to an atlas. The 3D ASMs of the bones are created automatically using a point distribution model optimization scheme. The accuracy and robustness of the segmentation approach was experimentally validated using an MR database of fat suppressed spoiled gradient recall images.

Metal-insulator transition and charge ordering in the extended Hubbard model at one-quarter filling

We study, with exact diagonalization, the zero temperature properties of the quarter-filled extended Hubbard model on a square lattice. We find that increasing the ratio of the intersite Coulomb repulsion, V, to the bandwidth drives the system from a metal to a charge ordered insulator. The evolution of the optical conductivity spectrum with increasing V is in agreement with the observed optical conductivity of several layered molecular crystals with the theta and beta crystal structures.

A uniform white-noise model for fixed-point roundoff errors in digital systems

Fixed-point roundoff noise in digital implementation of linear systems arises due to overflow, quantization of coefficients and input signals, and arithmetical errors. In uniform white-noise models, the last two types of roundoff errors are regarded as uniformly distributed independent random vectors on cubes of suitable size. For input signal quantization errors, the heuristic model is justified by a quantization theorem, which cannot be directly applied to arithmetical errors due to the complicated input-dependence of errors. The complete uniform white-noise model is shown to be valid in the sense of weak convergence of probabilistic measures as the lattice step tends to zero if the matrices of realization of the system in the state space satisfy certain nonresonance conditions and the finite-dimensional distributions of the input signal are absolutely continuous.

A comparison of model laws as a starting point for the development of an enforceable international consumer protection regime

This article is concerned primarily with an examination and comparison of select aspects of the model international consumer protection laws proposed by the United Nations (UN), the European Union (EU), and the Organisation for Economic Co-operation and Development (OECD), using the Trade Practices Act 1974 (Australia) as a basis for examination and comparison. As a secondary consideration, it also broadly examines the content of, and differences between, the model laws. The motive for this article is that any future enforceable international consumer protection regime (possibly in the form of an international treaty or convention) would need to take into account the UN, EU and OECD guidelines. A cross-comparison of those model laws, and a comparison of them with the consumer protection provisions of a well established national consumer protection law, should provide a useful starting point for the development of such a regime. The 'select aspects' of the model laws in question are the various provisions of those laws which could relate to situations involving the wrong delivery or non-delivery of goods.

Quantum trajectories for the realistic measurement of a solid-state charge qubit

We present a new model for the continuous measurement of a coupled quantum dot charge qubit. We model the effects of a realistic measurement, namely adding noise to, and filtering, the current through the detector. This is achieved by embedding the detector in an equivalent circuit for measurement. Our aim is to describe the evolution of the qubit state conditioned on the macroscopic output of the external circuit. We achieve this by generalizing a recently developed quantum trajectory theory for realistic photodetectors [P. Warszawski, H. M. Wiseman, and H. Mabuchi, Phys. Rev. A 65, 023802 (2002)] to treat solid-state detectors. This yields stochastic equations whose (numerical) solutions are the realistic quantum trajectories of the conditioned qubit state. We derive our general theory in the context of a low transparency quantum point contact. Areas of application for our theory and its relation to previous work are discussed.

Influence of charge size distribution on net-power draw of tumbling mill based on DEM modelling

Modelling and optimization of the power draw of large SAG/AG mills is important due to the large power draw which modern mills require (5-10 MW). The cost of grinding is the single biggest cost within the entire process of mineral extraction. Traditionally, modelling of the mill power draw has been done using empirical models. Although these models are reliable, they cannot model mills and operating conditions which are not within the model database boundaries. Also, due to its static nature, the impact of the changing conditions within the mill on the power draw cannot be determined using such models. Despite advances in computing power, discrete element method (DEM) modelling of large mills with many thousands of particles could be a time consuming task. The speed of computation is determined principally by two parameters: number of particles involved and material properties. The computational time step is determined by the size of the smallest particle present in the model and material properties (stiffness). In the case of small particles, the computational time step will be short, whilst in the case of large particles; the computation time step will be larger. Hence, from the point of view of time required for modelling (which usually corresponds to time required for 3-4 mill revolutions), it will be advantageous that the smallest particles in the model are not unnecessarily too small. The objective of this work is to compare the net power draw of the mill whose charge is characterised by different size distributions, while preserving the constant mass of the charge and mill speed. (C) 2004 Elsevier Ltd. All rights reserved.

Superconductivity mediated by charge fluctuations in layered molecular crystals

We consider the competition between superconducting, charge ordered, and metallic phases in layered molecular crystals with the theta and beta" structures. Applying slave-boson theory to the relevant extended Hubbard model, we show that the superconductivity is mediated by charge fluctuations and the Cooper pairs have d(xy) symmetry. This is in contrast to the kappa-(BEDT-TTF)(2)X family, for which theoretical calculations give superconductivity mediated by spin fluctuations and with d(x)2(-y)2 symmetry. We predict several materials that should become superconducting under pressure.

Ladder operator for the one-dimensional Hubbard model

The one-dimensional Hubbard model is integrable in the sense that it has an infinite family of conserved currents. We explicitly construct a ladder operator which can be used to iteratively generate all of the conserved current operators. This construction is different from that used for Lorentz invariant systems such as the Heisenberg model. The Hubbard model is not Lorentz invariant, due to the separation of spin and charge excitations. The ladder operator is obtained by a very general formalism which is applicable to any model that can be derived from a solution of the Yang-Baxter equation.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.