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.

Screening for leprosy in immigrants: A decision analysis model

Almost all leprosy cases reported in industrialized countries occur amongst immigrants or refugees from developing countries where leprosy continues to be an important health issue. Screening for leprosy is an important question for governments in countries with immigration and refugee programmes. A decision analysis framework is used to evaluate leprosy screening. The analysis uses a set of criteria and parameters regarding leprosy screening, and available data to estimate the number of cases which would be detected by a leprosy screening programme of immigrants from countries with different leprosy prevalences, compared with a policy of waiting for immigrants who develop symptomatic clinical diseases to present for health care. In a cohort of 100,000 immigrants from high leprosy prevalence regions (3.6/10,000), screening would detect 32 of the 42 cases which would arise in the destination country over the 14 years after migration; from medium prevalence areas (0.7/10,000) 6.3 of the total 8.1 cases would be detected, and from low prevalence regions (0.2/10,600) 1.8 of 2.3 cases. Using Australian data, the migrant mix would produce 74 leprosy cases from 10 years intake; screening would detect 54, and 19 would be diagnosed subsequently after migration. Screening would only produce significant case-yield amongst immigrants from regions or social groups with high leprosy prevalence. Since the number of immigrants to Australia from countries of higher endemnicity is not large routine leprosy screening would have a small impact on case incidence.

Incorporating demand-side aspects into regional policy: variations in the importance of private investment decision factors across regions

Most regional programs focus on the supply side of regions, emphasizing the attraction conditions offered, such as infrastructure, labor skills, tax incentives, etc. This study analyzes one aspect of the demand side, that is, how investment decisions of private firms are made by asking the question: ""Do corporations decide the same way on investments in different parts of the territory?"" The paper analyzes the investments of 373 large Brazilian firms during 1996-2004. Based on the investment decisions of these firms, the role of sales, cash-flow, external financing, and working capital is investigated through regression analysis. The regional influence is captured by explanatory variables representing regional and firm characteristics, and by interaction dummies between the region and the main investment determinants. The results indicate significant differences across regions in the importance of investment determinants. This information is important for regional development policy, because different mechanisms should be used in different regions to foster private investments.

A spectral clustering algorithm for manufacturing cell formation

A graph clustering algorithm constructs groups of closely related parts and machines separately. After they are matched for the least intercell moves, a refining process runs on the initial cell formation to decrease the number of intercell moves. A simple modification of this main approach can deal with some practical constraints, such as the popular constraint of bounding the maximum number of machines in a cell. Our approach makes a big improvement in the computational time. More importantly, improvement is seen in the number of intercell moves when the computational results were compared with best known solutions from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Decision rules for egg recognition are related to functional roles and chemical cues in the queenless ant Dinoponera quadriceps

The capacity to distinguish colony members from strangers is a key component in social life. In social insects, this extends to the brood and involves discrimination of queen eggs. Chemical substances communicate colony affiliation for both adults and brood; thus, in theory, all colony members should be able to recognize fellow nestmates. In this study, we investigate the ability of Dinoponera quadriceps workers to discriminate nestmate and non-nestmate eggs based on cuticular hydrocarbon composition. We analyzed whether cuticular hydrocarbons present on the eggs provide cues of discrimination. The results show that egg recognition in D. quadriceps is related to both age and the functional role of workers. Brood care workers were able to distinguish nestmate from non-nestmate eggs, while callow and forager workers were unable to do so.

Australia-wide predictions of soil properties using decision trees

This paper describes the construction of Australia-wide soil property predictions from a compiled national soils point database. Those properties considered include pH, organic carbon, total phosphorus, total nitrogen, thickness. texture, and clay content. Many of these soil properties are used directly in environmental process modelling including global climate change models. Models are constructed at the 250-m resolution using decision trees. These relate the soil property to the environment through a suite of environmental predictors at the locations where measurements are observed. These models are then used to extend predictions to the continental extent by applying the rules derived to the exhaustively available environmental predictors. The methodology and performance is described in detail for pH and summarized for other properties. Environmental variables are found to be important predictors, even at the 250-m resolution at which they are available here as they can describe the broad changes in soil property.

The effect of defendant gender on juror decision-making

Research investigating the role of stereotypes in jury decision-making has typically considered stereotypes as acting as peripheral cues in determin ing the credibility of experts or likelihood of guilt of defendants — with counter-stereotypic courtroom participants faring less well. The present study investigated the possibility that the extent to which courtroom participants are stereotypic can alter the mode of information processing. Students (N = 78) read a transcript of a case in which either a male or female allegedly committed an armed robbery. As predicted, the female counter-stereotypic defendant was distracting and effortful processing only occurred when the defendant was male. The male was seen as more guilty and the prosecution's case was more convincing when the prosecution had a strong, but not weak, case. There were no effects of case strength for the female defendant. Results are discussed in terms of the role of stereotypes in the jury decision-making.

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.