Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

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.

Advice on exercise from a family physician can help sedentary patients to become active

Objectives: To test the effectiveness, in the setting of primary health care, of verbal advice on exercise from a family physician (FP) combined with supporting written information. Design: A controlled trial with subjects allocated to a control group or one of two intervention groups using a balanced design based on day of the week. Setting: Ten general practices in Perth, Western Australia. Subjects: All sedentary patients consulting an FP. Intervention: Verbal advice on exercise from the FP and a pamphlet on exercise mailed to the patient's home address within 2 days of his/her visit to the doctor. Main outcome measure: Level of physical activity at followup. Results: 6,351 adult patients attending an FP practice completed a screening questionnaire, and 763 sedentary adults were recruited to the project. The response to follow-up, via a postal survey at 1, 6, and 12 months after the index consultation was 70%, 60%, and 57%, respectively. At 1 month a subsample of the control and intervention subjects were contacted for a telephone interview to verify self-reported levels of activity (n = 136). Treating all nonresponders as sedentary, at 1 month significantly more subjects in the combined intervention groups reported doing some physical activity (40%) compared with the control group (31%). Similarly, at 6 months, 30% of the control group and 38% of the combined intervention groups were now active. There was very little change at followup at 12 months (31% control and 36% intervention groups, respectively). Conclusion: A simple intervention aimed at the promotion of physical activity to sedentary patients in general practice can help reduce inactivity.

Two-level atom in a squeezed vacuum with finite bandwidth

The resonance fluorescence of a two-level atom driven by a coherent laser field and damped by a finite bandwidth squeezed vacuum is analysed. We extend the Yeoman and Barnett technique to a non-zero detuning of the driving field from the atomic resonance and discuss the role of squeezing bandwidth and the detuning in the level shifts, widths and intensities of the spectral lines. The approach is valid for arbitrary values of the Rabi frequency and detuning but for the squeezing bandwidths larger than the natural linewidth in order to satisfy the Markoff approximation. The narrowing of the spectral lines is interpreted in terms of the quadrature-noise spectrum. We find that, depending on the Rabi frequency, detuning and the squeezing phase, different factors contribute to the line narrowing. For a strong resonant driving field there is no squeezing in the emitted field and the fluorescence spectrum exactly reveals the noise spectrum. In this case the narrowing of the spectral lines arises from the noise reduction in the input squeezed vacuum. For a weak or detuned driving field the fluorescence exhibits a large squeezing and, as a consequence, the spectral lines have narrowed linewidths. Moreover, the fluorescence spectrum can be asymmetric about the central frequency despite the symmetrical distribution of the noise. The asymmetry arises from the absorption of photons by the squeezed vacuum which reduces the spontaneous emission. For an appropriate choice of the detuning some of the spectral lines can vanish despite that there is no population trapping. Again this process can be interpreted as arising from the absorption of photons by the squeezed vacuum. When the absorption is large it may compensate the spontaneous emission resulting in the vanishing of the fluorescence lines.

Spectral linewidth narrowing by a weak squeezed field of an arbitrary bandwidth

We study the spectral and noise properties of the fluorescence field emitted from a two-level atom driven by a beam of squeezed light. For a weak driving field we derive simple analytical formulae for the fluorescence and quadrature-noise spectra which are valid for an arbitrary bandwidth of the squeezed field. We analyse the spectra in the regime where the squeezing bandwidth is smaller or comparable to the atomic linewidth, the area where non-Markovian effects are important. We emphasize that there is a noticable difference between the fluorescence spectra for the thermal and squeezed field excitations. In both cases the spectrum can be narrower than any bandwidth involved in the process. However, as we point out for the squeezed driving field the linewidth narrowing, being much larger than in the thermal-field case, can be attributed to the squeezing of the fluctuations in the driving held. We also calculate the quadrature-noise spectrum of the emitted fluorescence, and find that for a detuned squeezed field the fluorescence spectrum does not reveal the quadrature-noise spectrum. In contrast to the fluorescence spectrum having two peaks, the quadrature-noise spectrum exhibits three peaks. We explain this difference as arising from the competiting three-photon scattering processes. (C) 1998 Elsevier Science B.V. All rights reserved.

Generalized optimal lattice covering of finite-dimensional Euclidean space

A generalization of the classical problem of optimal lattice covering of R-n is considered. Solutions to this generalized problem are found in two specific classes of lattices. The global optimal solution of the generalization is found for R-2. (C) 1998 Elsevier Science Inc. All rights reserved.

Determination of pindolol enantiomers in human plasma and urine by simple liquid-liquid extraction and high-performance liquid chromatography

A simple method for the measurement of pindolol enantiomers by HPLC is presented. Alkalinized serum or urine is extracted with ethyl acetate and the residue remaining after evaporation of the organic layer is then derivatised with (S)-(-)-alpha-methylbenzyl isocyanate. The diastereoisomers of derivatised pindolol and metoprolol (internal standard) are separated by high-performance liquid chromatography (HPLC) using a C-18 silica column and detected using fluorescence (excitation lambda: 215 nm, emission lambda: 320 nm). The assay displays reproducible linearity for pindolol enantiomers with a correlation coefficient of r(2) greater than or equal to 0.998 over the concentration range 8-100 ng ml(-1) for plasma and 0.1-2.5 mu g ml(-1) for urine. The coefficient of variation for accuracy and precision of the quality control samples for both plasma and urine are consistently

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Behavior and interactions of children in cooperative groups in lower and middle elementary grades

The study investigated the behaviors and interactions of children in structured and unstructured groups as they worked together on a 6-week social studies activity each term for 3 school terms. Two hundred and twelve children in Grade 1 and 184 children in Grade 3 participated in the study. Stratified random assignment occurred so that each gender-balanced group consisted of 1 high-, 2 medium-, and 1 low-ability student. The results show that the children in the structured groups were consistently more cooperative and they provided more elaborated and nonelaborated help than did their peers in the unstructured groups. The children in the structured groups in Grade 3 obtained higher reading and learning outcome scores than their peers in the unstructured groups.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Maintenance of cooperative and helping behaviors in reconstituted groups

The authors conducted a 1-year investigation of whether children, who had been trained in the previous school gear to cooperate, were able to use the skills they had been taught in reconstituted groups without additional training. Sixty-four 4th graders, who had participated in training in cooperative group behaviors in the previous gear, were assigned to the trained condition; 84 4th graders, who had not received any training, were assigned to the untrained condition. The children worked in mixed-ability (high, medium, low) gender-balanced groups (2 boys, 2 girls). Results indicate that the children in the trained groups were consistently more cooperative and helpful than their peers in the untrained groups, although they had not received refresher training in cooperative group behaviors.

A simple analysis of the US emission permits auctions

We examine a stylized version of EPA auctions when agents know the list of values of sellers and buyers. There are inefficient equilibria where no goods are traded and efficient equilibria where all exchange occurs at a uniform price. We also provide examples under incomplete information when the uniform price equilibrium holds and when it does not hold. (C) 1999 Elsevier Science S.A. All rights reserved. JEL classification: D44; Q29.