An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

The effect of tidal volume and peep on CO2 and oxygenation during resuscitation of very premature lambs

Background: Over-ventilation causing low arterial carbon dioxide levels (PaCO2) has been associated with the development of neonatal chronic lung disease and adverse outcomes. This may occur very soon after birth. Aim: To investigate the effect on PaCO2 and oxygenation of very premature lambs resuscitated with different tidal volumes and PEEP. Methods: Anaesthetised lambs delivered at 126 days gestation were randomised to 15 min resuscitation with 3 regimes: (1) Laerdal resuscitation bag (B) with 100% oxygen and no PEEP, (2) fixed tidal volume (VT) of 5 mL/kg, or (3) VT of 10 mL/kg, both delivered with a Babylog 8000 ventilator in volume guarantee mode with 8 cm H2O PEEP and variable FiO2. Frequent blood gases were measured and VT, mean airway pressure (Paw), minute volume (MV), ventilation rate (VR), respiratory system compliance (Crs) and alveolar/arterial oxygen difference (AaDO2) were recorded. Results: Twenty lambs were studied. B (1) was associated with more variable VT and peak inspiratory pressures (PIP) compared to fixed tidal volumes (2 and 3). The lambs ventilated with 10 mL/kg were over-ventilated, those ventilated with 5 mL/kg were slightly under-ventilated. Those ventilated with the Laerdal bag had a mean VT of 7.5 mL/kg and were normocarbic. The different tidal volumes had little effect on oxygenation. PEEP improved oxygenation. The table shows the values at 15 minutes expressed as mean and SEM. TABLE. No caption av... TABLE. No caption av... Image Tools Conclusion: Very premature lambs can be effectively resuscitated from birth using volume guarantee ventilation. Within minutes of birth different tidal volumes had a large effect on PaCO2 and no effect on oxygenation. Studies are needed to determine the appropriate tidal volume for resuscitating very premature infants to maintain acceptable levels of PaCO2. © International Pediatrics Research Foundation, Inc. 2004. All Rights Reserved.

Modeling of columnar-to-equiaxed transition in solidified Al-Si alloys

A Cellular-Automaton Finite-Volume-Method (CAFVM) algorithm has been developed, coupling with macroscopic model for heat transfer calculation and microscopic models for nucleation and growth. The solution equations have been solved to determine the time-dependent constitutional undercooling and interface retardation during solidification. The constitutional undercooling is then coupled into the CAFVM algorithm to investigate both the effects of thermal and constitutional undercooling on columnar growth and crystal selection in the columnar zone, and formation of equiaxed crystals in the bulk liquid. The model cannot only simulate microstructures of alloys but also investigates nucleation mechanisms and growth kinetics of alloys solidified with various solute concentrations and solidification morphologies.

Effects of tidal volume and positive end-expiratory pressure during resuscitation of very premature lambs

Background: Guidelines recommend neonatal resuscitation without controlling tidal volume or positive end-expiratory pressure (PEEP). However, these may improve gas exchange, lung volume and outcome. Aim: To investigate resuscitation of very premature lambs with a Laerdal bag without PEEP versus volume guarantee ventilation with PEEP. Methods: Anaesthetized lambs (n = 20) delivered at 125 d gestation were randomized to three groups receiving 15 min resuscitation: (1) Laerdal bag and no PEEP; (2) ventilation with a tidal volume of 5 ml/kg and 8 cm H2O PEEP; (3) ventilation with 10 ml/kg and 8 cm H2O PEEP. They were then all ventilated for 2 h with tidal volumes of 5 or 10 ml/kg, and 8 cm H2O PEEP. Ventilation parameters and blood gases were recorded. Results: Different tidal volumes affected PaCO2 within minutes, with 10 ml/kg causing severe hypocarbia. PEEP had little effect on PaCO2. Oxygenation improved significantly with PEEP of 8 cm H2O, irrespective of tidal volume. Conclusion: Very premature lambs can be resuscitated effectively using volume-guarantee ventilation and PEEP. Tidal volumes affected PaCO2 within minutes but had little effect on oxygenation. PEEP halved the oxygen requirement compared with no PEEP. Resuscitating premature babies with controlled tidal volumes and PEEP might improve their outcome.

Cancer detection and mammogram volume of radiologists in a population-based screening programme

This study investigates the relationship between the number of screening mammograms read by radiologists and the screening breast cancer detection rate. Cancer detection rates for incident screens (all women aged >= 40 years) were compared by increasing categories of reader volume using Poisson regression. Data from New South Wales (NSW) for a 2 year period (2000-2001) were obtained from the BreastScreen NSW programme. Cancer detection rates increased with the number of mammograms read in the programme, reaching a plateau of approximately 40 per 10,000 after 1375 mammograms per year. No significant differences in cancer detection were evident above 875 mammograms (compared to below 875 mammograms) per year (RR = 0.79, 95% CI 0.63-0.99). (c) 2005 Elsevier Ltd. All rights reserved.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

The Qu-Prolog unification algorithm: Formalisation and correctness

Qu-Prolog is an extension of Prolog which performs meta-level computations over object languages, such as predicate calculi and lambda-calculi, which have object-level variables, and quantifier or binding symbols creating local scopes for those variables. As in Prolog, the instantiable (meta-level) variables of Qu-Prolog range over object-level terms, and in addition other Qu-Prolog syntax denotes the various components of the object-level syntax, including object-level variables. Further, the meta-level operation of substitution into object-level terms is directly represented by appropriate Qu-Prolog syntax. Again as in Prolog, the driving mechanism in Qu-Prolog computation is a form of unification, but this is substantially more complex than for Prolog because of Qu-Prolog's greater generality, and especially because substitution operations are evaluated during unification. In this paper, the Qu-Prolog unification algorithm is specified, formalised and proved correct. Further, the analysis of the algorithm is carried out in a frame-work which straightforwardly allows the 'completeness' of the algorithm to be proved: though fully explicit answers to unification problems are not always provided, no information is lost in the unification process.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Simulated-annealing-based genetic algorithm for modeling the optical constants of solids

We propose a simulated-annealing-based genetic algorithm for solving model parameter estimation problems. The algorithm incorporates advantages of both genetic algorithms and simulated annealing. Tests on computer-generated synthetic data that closely resemble optical constants of a metal were performed to compare the efficiency of plain genetic algorithms against the simulated-annealing-based genetic algorithms. These tests assess the ability of the algorithms to and the global minimum and the accuracy of values obtained for model parameters. Finally, the algorithm with the best performance is used to fit the model dielectric function to data for platinum and aluminum. (C) 1997 Optical Society of America.