The critical power function is dependent on the duration of the predictive exercise tests chosen

The linear relationship between work accomplished (W-lim) and time to exhaustion (t(lim)) can be described by the equation: W-lim = a + CP.t(lim). Critical power (CP) is the slope of this line and is thought to represent a maximum rate of ATP synthesis without exhaustion, presumably an inherent characteristic of the aerobic energy system. The present investigation determined whether the choice of predictive tests would elicit significant differences in the estimated CP. Ten female physical education students completed, in random order and on consecutive days, five art-out predictive tests at preselected constant-power outputs. Predictive tests were performed on an electrically-braked cycle ergometer and power loadings were individually chosen so as to induce fatigue within approximately 1-10 mins. CP was derived by fitting the linear W-lim-t(lim) regression and calculated three ways: 1) using the first, third and fifth W-lim-t(lim) coordinates (I-135), 2) using coordinates from the three highest power outputs (I-123; mean t(lim) = 68-193 s) and 3) using coordinates from the lowest power outputs (I-345; mean t(lim) = 193-485 s). Repeated measures ANOVA revealed that CPI123 (201.0 +/- 37.9W) > CPI135 (176.1 +/- 27.6W) > CPI345 (164.0 +/- 22.8W) (P < 0.05). When the three sets of data were used to fit the hyperbolic Power-t(lim) regression, statistically significant differences between each CP were also found (P < 0.05). The shorter the predictive trials, the greater the slope of the W-lim-t(lim) regression; possibly because of the greater influence of 'aerobic inertia' on these trials. This may explain why CP has failed to represent a maximal, sustainable work rate. The present findings suggest that if CP is to represent the highest power output that an individual can maintain for a very long time without fatigue then CP should be calculated over a range of predictive tests in which the influence of aerobic inertia is minimised.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Decoherence in ion traps due to laser intensity and phase fluctuations

We consider one source of decoherence for a single trapped ion due to intensity and phase fluctuations in the exciting laser pulses. For simplicity we assume that the stochastic processes involved are white noise processes, which enables us to give a simple master equation description of this source of decoherence. This master equation is averaged over the noise, and is sufficient to describe the results of experiments that probe the oscillations in the electronic populations as energy is exchanged between the internal and electronic motion. Our results are in good qualitative agreement with recent experiments and predict that the decoherence rate will depend on vibrational quantum number in different ways depending on which vibrational excitation sideband is used.

Integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons

A general graded reflection equation algebra is proposed and the corresponding boundary quantum inverse scattering method is formulated. The formalism is applicable to all boundary lattice systems where an invertible R-matrix exists. As an application, the integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons are investigated. The diagonal boundary K-matrices are found and a class of integrable boundary terms are determined. The boundary system is solved by means of the coordinate space Bethe ansatz technique and the Bethe ansatz equations are derived. As a sideline, it is shown that all R-matrices associated with a quantum affine superalgebra enjoy the crossing-unitarity property. (C) 1998 Elsevier Science B.V.

Predicting sprint running times from isokinetic and squat lift tests: A regression analysis

This study examined the relationship between isokinetic hip extensor/hip flexor strength, 1-RM squat strength, and sprint running performance for both a sprint-trained and non-sprint-trained group. Eleven male sprinters and 8 male controls volunteered for the study. On the same day subjects ran 20-m sprints from both a stationary start and with a 50-m acceleration distance, completed isokinetic hip extension/flexion exercises at 1.05, 4.74, and 8.42 rad.s(-1), and had their squat strength estimated. Stepwise multiple regression analysis showed that equations for predicting both 20-m maximum velocity nm time and 20-m acceleration time may be calculated with an error of less than 0.05 sec using only isokinetic and squat strength data. However, a single regression equation for predicting both 20-m acceleration and maximum velocity run times from isokinetic or squat tests was not found. The regression analysis indicated that hip flexor strength at all test velocities was a better predictor of sprint running performance than hip extensor strength.

On the graded quantum Yang-Baxter and reflection equations

We clarify the extra signs appearing in the graded quantum Yang-Baxter reflection equations, when they are written in a matrix form. We find the boundary K-matrix for the Perk-Schultz six-vertex model, thus give a general solution to the graded reflection equation associated with it.

Coalbed methane sorption related to coal composition

Gas sorption by coal is closely related to its physical and chemical properties, which are, in turn, governed by coal type and rank. The role of coal type (sensu maceral composition) is not fully established but it is clear that coal type may affect both adsorption capacity and desorption rate. Adsorption capacity is closely related to micropore (pores <2 nm) development, which is rank and maceral dependent. Adsorption isotherms indicate that in most cases bright (vitrinite-rich) coals have a greater adsorption capacity than their dull (often inertinite-rich) equivalents. However, no differences, or even the opposing trend, may be observed in relation to coal type. Desorption rate investigations have been performed using selected bright and dull coal samples in a high pressure microbalance. Interpretation of results using unipore spherical and bidisperse pore models indicate the importance of the pore structure. Bright, vitrinite-rich coals usually have the slowest desorption rates which is associated with their highly microporous structure. However, rapid desorption in bright coals may be related to development of extensive, unmineralised fracture systems. Both macro-and micro-pore systems are implicated in the more rapidly desorbing dull coals. Some dull, inertinite-rich coals may rapidly desorb due to a predominance of large, open cell lumina. Mineral matter is essentially nonadsorbent to coal gases and acts as a simple diluent. However, mineral-rich coals may be associated with more rapid desorption. Coal rank and type (maceral composition) per se do not appear to be the critical factors in controlling gas sorption, but rather the influence they exert over pore structure development. (C) 1998 Elsevier Science B.V.

Probability backflow for a Dirac particle

The phenomenon of probability backflow, previously quantified for a free nonrelativistic particle, is considered for a free particle obeying Dirac's equation. It is known that probability backflow can occur in the opposite direction to the momentum; that is to say, there exist positive-energy states in which the particle certainly has a positive momentum in a given direction, but for which the component of the probability flux vector in that direction is negative. It is shown thar the maximum possible amount of probability that can flow backwards, over a given time interval of duration T, depends on the dimensionless parameter epsilon = root 4h/mc(2)T, where m is the mass of the particle and c is the speed of light. At epsilon = 0, the nonrelativistic value of approximately 0.039 for this maximum is recovered. Numerical studies suggest that the maximum decreases monotonically as epsilon increases from 0, and show that it depends on the size of m, h, and T, unlike the nonrelativistic case.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

Calculating current densities and fields due to shielded bi-planar radio-frequency coils

A method for the accurate computation of the current densities produced in a wide-runged bi-planar radio-frequency coil is presented. The device has applications in magnetic resonance imaging. There is a set of opposing primary rungs, symmetrically placed on parallel planes and a similar arrangement of rungs on two parallel planes surrounding the primary serves as a shield. Current densities induced in these primary and shielding rungs are calculated to a high degree of accuracy using an integral-equation approach, combined with the inverse finite Hilbert transform. Once these densities are known, accurate electrical and magnetic fields are then computed without difficulty. Some test results are shown. The method is so rapid that it can be incorporated into optimization software. Some preliminary fields produced from optimized coils are presented.

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

Thickness and stability of adsorbed film in cylindrical mesopores

The adsorbed film in small cylindrical mesopores is studied by using MCM-41 samples of uniform cylindrical channels as model systems. It is found that at a given relative pressure, the smaller the pore radius, the thicker the adsorbed film is, as postulated by Broekhoff and De Beer. Thermodynamics analysis established that the stability of the adsorbed film is determined by interface curvature and the potential of interaction between adsorbate and adsorbent. A semiempirical equation is proposed to describe the state of stable adsorbed films in cylindrical mesopores. It is also shown to be useful in calculations of pore size distributions of mesoporous solids.

Epidermal iontophoresis: I. Development of the ionic mobility-pore model

Purpose, An integrated ionic mobility-pore model for epidermal iontophoresis is developed from theoretical considerations using both the free volume and pore restriction forms of the model for a range of solute radii (r(j)) approaching the pore radii (r(p)) as well as approximation of the pore restriction form for r(j)/r(p) < 0.4. In this model, we defined the determinants for iontophoresis as solute size (defined by MV, MW or radius), solute mobility, solute shape, solute charge, the Debye layer thickness, total current applied, solute concentration, fraction ionized, presence of extraneous ions (defined by solvent conductivity), epidermal permselectivity, partitioning rates to account for interaction of unionized and ionized lipophilic solutes with the wall of the pore and electroosmosis. Methods, The ionic mobility-pore model was developed from theoretical considerations to include each of the determinants of iontophoretic transport. The model was then used to reexamine iontophoretic flux conductivity and iontophoretic flux-fraction ionized literature data on the determinants of iontophoretic flux. Results. The ionic mobility-pore model was found to be consistent with existing experimental data and determinants defining iontophoretic transport. However, the predicted effects of solute size on iontophoresis are more consistent with the pore-restriction than free volume form of the model. A reanalysis of iontophoretic flux-conductivity data confirmed the model's prediction that, in the absence of significant electroosmosis, the reciprocal of flux is linearly related to either donor or receptor solution conductivity. Significant interaction with the pore walls, as described by the model, accounted for the reported pH dependence of the iontophoretic transport for a range of ionizable solutes. Conclusions. The ionic mobility-pore iontophoretic model developed enables a range of determinants of iontophoresis to be described in a single unifying equation which recognises a range of determinants of iontophoretic flux.

Theory of multidimensional parametric band-gap simultons

Multidimensional spatiotemporal parametric simultons (simultaneous solitary waves) are possible in a nonlinear chi((2)) medium with a Bragg grating structure, where large effective dispersion occurs near two resonant band gaps for the carrier and second-harmonic field, respectively. The enhanced dispersion allows much reduced interaction lengths, as compared to bulk medium parametric simultons. The nonlinear parametric band-gap medium permits higher-dimensional stationary waves to form. In addition, solitons can occur with lower input powers than conventional nonlinear Schrodinger equation gap solitons. In this paper, the equations for electromagnetic propagation in a grating structure with a parametric nonlinearity are derived from Maxwell's equation using a coupled mode Hamiltonian analysis in one, two, and three spatial dimensions. Simultaneous solitary wave solutions are proved to exist by reducing the equations to the coupled equations describing a nonlinear parametric waveguide, using the effective-mass approximation (EMA). Exact one-dimensional numerical solutions in agreement with the EMA solutions are also given. Direct numerical simulations show that the solutions have similar types of stability properties to the bulk case, providing the carrier waves are tuned to the two Bragg resonances, and the pulses have a width in frequency space less than the band gap. In summary, these equations describe a physically accessible localized nonlinear wave that is stable in up to 3 + 1 dimensions. Possible applications include photonic logic and switching devices. [S1063-651X(98)06109-1].