A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Fibroblast growth factor receptor 4 (FGFR4) expression in newborn murine calvaria and primary osteoblast cultures

Fibroblast growth factor receptor (FGFR) signalling is important in the initiation and regulation of osteogenesis. Although mutations in FGFR1, 2 and 3 genes are known to cause skeletal deformities, the expression of FGFR4 in bony tissue remains unclear. We have investigated the expression pattern of FGFR4 in the neonatal mouse calvaria and compared it to the expression pattern in cultures of primary osteoblasts. Immunohistochemistry demonstrated that FGFR4 was highly expressed in rudimentary membranous bone and strictly localised to the cellular components (osteoblasts) between the periosteal and endosteal layers. Cells in close proximity to the newly formed osteoid (preosteoblasts) also expressed FGFR4 on both the endosteal and periosteal surfaces. Immunocytochemical analysis of primary osteoblast cultures taken from the same cranial region also revealed high levels of FGFR4 expression, suggesting a similar pattern of cellular expression in vivo and in vitro. RT-PCR and Western blotting for FGFR4 confirmed its presence in primary osteoblast cultures. These results suggest that FGFR4 may be an important regulator of osteogenesis with involvement in preosteoblast proliferation and differentiation as well as osteoblast functioning during intramembranous ossification. The consistent expression of FGFR4 in vivo and in vitro supports the use of primary osteoblast cultures for elucidating the role of FGFR4 during osteogenesis.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Dinoflagellate symbiosis: strategies and adaptations for the aquisition and fixation of inorganic carbon

Dinoflagellates exist in symbiosis with a number of marine invertebrates including giant clams, which are the largest of these symbiotic organisms. The dinoflagellates (Symbiodinium sp.) live intercellularly within tubules in the mantle of the host clam. The transport of inorganic carbon (Ci) from seawater to Symbiodinium (=zooxanthellae) is an essential function of hosts that derive the majority of their respiratory energy from the photosynthate exported by the zooxanthellae. Immunolocalisation studies show that the host has adapted its physiology to acquire, rather than remove CO2, from the haemolymph and clam tissues. Two carbonic anhydrase (CA) isoforms (32 and 70 kDa) play an essential part in this process. These have been localised to the mantle and gill tissues where they catalyse the interconversion of HCO3- to CO2, which then diffuses into the host tissues. The zooxanthellae exhibit a number of strategies to maximise Ci acquisition and utilisation. This is necessary as they express a form II Rubisco that has poor discrimination between CO2 and O-2. Evidence is presented for a carbon concentrating mechanism (CCM) to overcome. this disadvantage. The CCM incorporates the presence of a light-activated CA activity, a capacity to take up both HCO3- and CO2, an ability to accumulate an elevated concentration of Ci within the algal cell, and localisation of Rubisco to the pyrenoid. These algae also express both external and intracellular CAs, with the intracellular isoforms being localised to the thylakoid lumen and pyrenoid. These results have been incorporated into a model that explains the transport of Ci from seawater through the clam to the zooxanthellae.

Modelling activated sludge flocculation using population balances

Activated sludge flocculation was modelled using population balances. The model followed the dynamics of activated sludge flocculation providing a good approximation of the change in mean floe size with time. Increasing the average velocity gradient decreased the final floe size. The breakage rate coefficient and collision efficiency also varied with the average velocity gradient. A power law relationship was found for the increase in breakage rate coefficient with increasing average velocity gradient. Further investigation will be conducted to determine the relationship between the collision efficiency and particle size to provide a better approximation of dynamic changes in the floe size distribution during flocculation. (C) 2002 Published by Elsevier Science B.V.

SAM (dependent) I AM: the S-adenosylmethionine-dependent methyltransferase fold

The S-adenosylmethionine-dependent methyltransferase enzymes share little sequence identity, but incorporate a highly conserved structural fold. Surprisingly, residues that bind the common cofactor are poorly conserved, although the binding site is localised to the same region of the fold. The substrate-binding region of the fold varies enormously. Over the past two years, there has been a significant increase in the number of structures that are known to incorporate this fold, including several uncharacterised proteins and two proteins that lack methyltransferase activity.

A new method for analysing the equilibrium and time-dependent behaviour of Markovian models

Many large-scale stochastic systems, such as telecommunications networks, can be modelled using a continuous-time Markov chain. However, it is frequently the case that a satisfactory analysis of their time-dependent, or even equilibrium, behaviour is impossible. In this paper, we propose a new method of analyzing Markovian models, whereby the existing transition structure is replaced by a more amenable one. Using rates of transition given by the equilibrium expected rates of the corresponding transitions of the original chain, we are able to approximate its behaviour. We present two formulations of the idea of expected rates. The first provides a method for analysing time-dependent behaviour, while the second provides a highly accurate means of analysing equilibrium behaviour. We shall illustrate our approach with reference to a variety of models, giving particular attention to queueing and loss networks. (C) 2003 Elsevier Ltd. All rights reserved.

Existence of multiple solutions for finite difference approximations to second-order boundary value problems

Error condition detected We consider discrete two-point boundary value problems of the form D-2 y(k+1) = f (kh, y(k), D y(k)), for k = 1,...,n - 1, (0,0) = G((y(0),y(n));(Dy-1,Dy-n)), where Dy-k = (y(k) - Yk-I)/h and h = 1/n. This arises as a finite difference approximation to y" = f(x,y,y'), x is an element of [0,1], (0,0) = G((y(0),y(1));(y'(0),y'(1))). We assume that f and G = (g(0), g(1)) are continuous and fully nonlinear, that there exist pairs of strict lower and strict upper solutions for the continuous problem, and that f and G satisfy additional assumptions that are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. Under these assumptions we show that there are at least three distinct solutions of the discrete approximation which approximate solutions to the continuous problem as the grid size, h, goes to 0. (C) 2003 Elsevier Science Ltd. All rights reserved.

Product form approximations for highly linear loss networks with trunk reservation

Admission controls, such as trunk reservation, are often used in loss networks to optimise their performance. Since the numerical evaluation of performance measures is complex, much attention has been given to finding approximation methods. The Erlang Fixed-Point (EFP) approximation, which is based on an independent blocking assumption, has been used for networks both with and without controls. Several more elaborate approximation methods which account for dependencies in blocking behaviour have been developed for the uncontrolled setting. This paper is an exploratory investigation of extensions and synthesis of these methods to systems with controls, in particular, trunk reservation. In order to isolate the dependency factor, we restrict our attention to a highly linear network. We will compare the performance of the resulting approximations against the benchmark of the EFP approximation extended to the trunk reservation setting. By doing this, we seek to gain insight into the critical factors in constructing an effective approximation. (C) 2003 Elsevier Ltd. All rights reserved.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

Biaxial failure model for fiber reinforced concrete

Steel fiber reinforced concrete (SFRC) is widely applied in the construction industry. Numerical elastoplastic analysis of the macroscopic behavior is complex. This typically involves a piecewise linear failure curve including corner singularities. This paper presents a single smooth biaxial failure curve for SFRC based on a semianalytical approximation. Convexity of the proposed model is guaranteed so that numerical problems are avoided. The model has sufficient flexibility to closely match experimental results. The failure curve is also suitable for modeling plain concrete under biaxial loading. Since this model is capable of simulating the failure states in all stress regimes with a single envelope, the elastoplastic formulation is very concise and simple. The finite element implementation is developed to demonstrate the conciseness and the effectiveness of the model. The computed results display good agreement with published experimental data.

Clinical aspects of gestational trophoblastic disease: A review based partly on 25-year experience of a statewide registry

Background: Gestational trophoblastic disease is a fascinating group of pregnancy disorders characterised by abnormal proliferation of trophoblast, ranging from benign to malignant. Because the disease is uncommon, there is a need to formulate management with the assistance of collective information. Methodology: A review of available information from English written literature was undertaken especially data reported by registries around the world (Charing Cross Hospital in England, the North-western University and the New England area in the USA as well as our own experience in Queensland, Australia). Where possible, collated data from relevant studies were analysed to answer some of the questions posed in clinical practice, with reference to metastatic disease to liver and brain, twinning of molar gestation and coexisting fetus, and placental-site tumour. Results: We found that molar gestation can be classified according to its clinical presentation which influences the time taken to reach human chorionic gonadotropin (HCG) 'negativity' and the risk of persisting disease. Categorisation of risk is the basis for choice of chemotherapy to achieve good outcomes. Metastases to liver and brain remain problems in management; the development of 'new' metastases during chemotherapy is a very poor prognostic factor. In the variant of twinning with molar gestation and coexisting fetus, it is important to elucidate the fetal karyotype in planning management: a 69XXX fetus is not salvageable but a normal 46XX or 46XY fetus faces the prospect of early preterm delivery. The placental-site tumour is very rare; localised disease is curable by surgery; chemotherapy is less effective in disseminated disease. From collated worldwide data, the recurrence rate after one mole is 1.3% and after two or more is 20%. Reproductive outcome in subsequent pregnancies, even after multidrug chemotherapy, is not different from the general population. Because of the increased risk long-term of second tumours after multidrug chemotherapy a closer surveillance of these patients is necessary Conclusion: In general, the disease in its persisting or malignant form is 'a cancer model par excellence' because of an identifiable precursor condition, a reliable HCG marker, and sensitivity of the disease to cytotoxic drugs. With current management, retention of fertility is possible and normal reproductive outcome assured.

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

Photon-added detection

The production of conditional quantum states and quantum operations based on the result of measurement is now seen as a key tool in quantum information and metrology. We propose a different type of photon number detector. It functions nondeterministically, but when successful, it has high fidelity. The detector, which makes use of an n-photon auxiliary Fock state and high efficiency homodyne detection, allows a tunable trade-off between fidelity and probability. By sacrificing probability of operation, an excellent approximation to a photon-number detector is achieved.

Resolution enhancement by subtraction of confocal signals taken at different pinhole sizes

Subtractive imaging in confocal fluorescence light microscopy is based on the subtraction of a suitably weighted widefield image from a confocal image. An approximation to a widefield image can be obtained by detection with an opened confocal pinhole. The subtraction of images enhances the resolution in-plane as well as along the optic axis. Due to the linearity of the approach, the effect of subtractive imaging in Fourier-space corresponds to a reduction of low spatial frequency contributions leading to a relative enhancement of the high frequencies. Along the direction of the optic axis this also results in an improved sectioning. Image processing can achieve a similar effect. However, a 3D volume dataset must be acquired and processed, yielding a result essentially identical to subtractive imaging but superior in signal-to-noise ratio. The latter can be increased further with the technique of weighted averaging in Fourier-space. A comparison of 2D and 3D experimental data analysed with subtractive imaging, the equivalent Fourier-space processing of the confocal data only, and Fourier-space weighted averaging is presented. (C) 2003 Elsevier Ltd. All rights reserved.