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.

Optical precursors and Beer's law violations; non-exponential propagation losses in water

We develop a model for exponential decay of broadband pulses, and examine its implications for experiments on optical precursors. One of the signature features of Brillouin precursors is attenuation with a less rapid decay than that predicted by Beer's Law. Depending on the pulse parameters and the model that is adopted for the dielectric properties of the medium, the limiting z-dependence of the loss has been described as z(-1/2), z(-1/3), exponential, or, in more detailed descriptions, some combination of the above. Experimental results in the search for precursors are examined in light of the different models, and a stringent test for sub-exponential decay is applied to data on propagation of 500 femtosecond pulses through 1-5 meters of water. (C) 2005 Optical Society of America.

A scalarization technique for computing the power and exponential moments of gaussian random matrices

We consider the problems of computing the power and exponential moments EXs and EetX of square Gaussian random matrices X=A+BWC for positive integer s and real t, where W is a standard normal random vector and A, B, C are appropriately dimensioned constant matrices. We solve the problems by a matrix product scalarization technique and interpret the solutions in system-theoretic terms. The results of the paper are applicable to Bayesian prediction in multivariate autoregressive time series and mean-reverting diffusion processes.

A worldwide correlation for exponential bed particle size variation in subaerial aqueous flows

The particle size of the bed sediments in or on many natural streams, alluvial fans, laboratory flumes, irrigation canals and mine waste deltas varies exponentially with distance along the stream. A plot of the available worldwide exponential bed particle size diminution coefficient data against stream length is presented which shows that all the data lie within a single narrow band extending over virtually the whole range of stream lengths and bed sediment particle sizes found on Earth. This correlation applies to both natural and artificial flows with both sand and gravel beds, irrespective of either the solids concentration or whether normal or reverse sorting occurs. This strongly suggests that there are common mechanisms underlying the exponential diminution of bed particles in subaerial aqueous flows of all kinds. Thus existing models of sorting and abrasion applicable to some such flows may be applicable to others. A comparison of exponential laboratory abrasion and field diminution coefficients suggests that abrasion is unlikely to be significant in gravel and sand bed streams shorter than about 10 km to 100 km, and about 500 km, respectively. Copyright (C) 1999 John Wiley & Sons, Ltd.

Spin-glass dynamics

Spin glasses are magnetic systems with conflicting and random interactions between the individual spins. The dynamics of spin glasses, as of structural glasses, reflect their complexity. Both in experimental and numerical work the relaxation below the freezing temperature depends strongly on the annealing conditions (aging) and, above the freezing point, relaxation in equilibrium is slow and non-exponential, In this Forum, observed characteristics of the dynamics were summarized and the physical models proposed to explain them were outlined. (C) 1998 Elsevier Science B.V. All rights reserved.

Anthocyanin synthesis, growth and nutrient uptake in suspension cultures of strawberry cells

Strawberry (Fragaria ananassa cv. Shikinari) cell suspension cultures carried out in shake flasks for 18 d were closely examined for cell growth, anthocyanin synthesis and the development of pigmented cells in relation to the uptake of carbohydrate, extracellular PO4, NO3, NH4, and calcium. Cell viability, extracellular anthocyanin content, pH and electrical conductivity of the broth were also monitored. The specific growth rate of strawberry cells at exponential phase was 0.27 and 0.28 d(-1) based on fresh and dry weight, respectively. Anthocyanin synthesis was observed to increase continuously to a maximum value of 0.86 mg/g fresh cell weight (FCW) at day 6, and was partially growth-associated. Anthocyanin synthesis was linearly related to the increase in pigmented cell ratio, which increased with time and reached a maximum value of ca. 70% at day 6 due to reduction in cell viability and depletion of substrate. Total carbohydrate uptake was closely associated with increase in cell growth, and glucose was utilized in preference to fructose. Nitrate and ammonia were consumed until 9 d of culture, but phosphate was completely absorbed within 4 d. Calcium was assimilated throughout the growth cycle. After 9 d, cell lysis was observed which resulted in the leakage of intracellular substances and a concomitant pH rise. Anthocyanin was never detected in the broth although the broth became darkly pigmented during the lysis period. This suggests that anthocyanin was synthesized only by viable pigmented cells, and degraded rapidly upon cell death and lysis. Based on the results of kinetic analysis, a model was developed by incorporating governing equations for the ratio of pigmented cells into a Bailey and Nicholson's model. This was verified by comparison with the experimental data. The results suggest Bat the model satisfactorily describes the strawberry cell culture process, and may thus be used for process optimization.

Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Control of fed-batch fermentations

Fed-batch fermentation is used to prevent or reduce substrate-associated growth inhibition by controlling nutrient supply. Here we review the advances in control of fed-batch fermentations. Simple exponential feeding and inferential methods are examined, as are newer methods based on fuzzy control and neural networks. Considerable interest has developed in these more advanced methods that hold promise for optimizing fed-batch techniques for complex fermentation systems. (C) 1999 Elsevier Science Inc. All rights reserved.

Rectification of the olfactory cyclic nucleotide-gated channel by intracellular polyamines

Polyamine-induced inward rectification of cyclic nucleotide-gated channels was studied in inside-out patches from rat olfactory neurons. The polyamines, spermine, spermidine and putrescine, induced an 'instantaneous' voltage-dependent inhibition with K-d values at 0 mV of 39, 121 mu M and 2.7 mM, respectively. Hill coefficients for inhibition were significantly < 1, suggesting an allosteric inhibitory mechanism. The Woodhull model for voltage-dependent block predicted that all 3 polyamines bound to a site 1/3 of the electrical distance through the membrane from the internal side. Instantaneous inhibition was relieved at positive potentials, implying significant polyamine permeation. Spermine also induced exponential current relaxations to a 'steady-state' impermeant level. This inhibition was also mediated by a binding site 1/3 of the electrical distance through the pore, but with a K-d of 2.6 mM. Spermine inhibition was explained by postulating two spermine binding sites at a similar depth. Occupation of the first site occurs rapidly and with high affinity, but once a spermine molecule has bound, it inhibits spermine occupation of the second binding site via electrostatic repulsion. This repulsion is overcome at higher membrane potentials, but results in a lower apparent binding affinity for the second spermine molecule. The on-rate constant for the second spermine binding saturated at a low rate (similar to 200 sec(-1) at +120 mV), providing further evidence for an allosteric mechanism. Polyamine-induced inward rectification was significant at physiological concentrations.