An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Integrating spatial optimization and cellular automata for evaluating urban change

Urban growth and change presents numerous challenges for planners and policy makers. Effective and appropriate strategies for managing growth and change must address issues of social, environmental and economic sustainability. Doing so in practical terms is a difficult task given the uncertainty associated with likely growth trends not to mention the uncertainty associated with how social and environmental structures will respond to such change. An optimization based approach is developed for evaluating growth and change based upon spatial restrictions and impact thresholds. The spatial optimization model is integrated with a cellular automata growth simulation process. Application results are presented and discussed with respect to possible growth scenarios in south east Queensland, Australia.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Simulation of air traffic controllers' behaviour using the operator choice model

The Operator Choice Model (OCM) was developed to model the behaviour of operators attending to complex tasks involving interdependent concurrent activities, such as in Air Traffic Control (ATC). The purpose of the OCM is to provide a flexible framework for modelling and simulation that can be used for quantitative analyses in human reliability assessment, comparison between human computer interaction (HCI) designs, and analysis of operator workload. The OCM virtual operator is essentially a cycle of four processes: Scan Classify Decide Action Perform Action. Once a cycle is complete, the operator will return to the Scan process. It is also possible to truncate a cycle and return to Scan after each of the processes. These processes are described using Continuous Time Probabilistic Automata (CTPA). The details of the probability and timing models are specific to the domain of application, and need to be specified using domain experts. We are building an application of the OCM for use in ATC. In order to develop a realistic model we are calibrating the probability and timing models that comprise each process using experimental data from a series of experiments conducted with student subjects. These experiments have identified the factors that influence perception and decision making in simplified conflict detection and resolution tasks. This paper presents an application of the OCM approach to a simple ATC conflict detection experiment. The aim is to calibrate the OCM so that its behaviour resembles that of the experimental subjects when it is challenged with the same task. Its behaviour should also interpolate when challenged with scenarios similar to those used to calibrate it. The approach illustrated here uses logistic regression to model the classifications made by the subjects. This model is fitted to the calibration data, and provides an extrapolation to classifications in scenarios outside of the calibration data. A simple strategy is used to calibrate the timing component of the model, and the results for reaction times are compared between the OCM and the student subjects. While this approach to timing does not capture the full complexity of the reaction time distribution seen in the data from the student subjects, the mean and the tail of the distributions are similar.

Combining influence maps and cellular automata for reactive game agents

Agents make up an important part of game worlds, ranging from the characters and monsters that live in the world to the armies that the player controls. Despite their importance, agents in current games rarely display an awareness of their environment or react appropriately, which severely detracts from the believability of the game. Some games have included agents with a basic awareness of other agents, but they are still unaware of important game events or environmental conditions. This paper presents an agent design we have developed, which combines cellular automata for environmental modeling with influence maps for agent decision-making. The agents were implemented into a 3D game environment we have developed, the EmerGEnT system, and tuned through three experiments. The result is simple, flexible game agents that are able to respond to natural phenomena (e.g. rain or fire), while pursuing a goal.

The use of calcium hydroxide, antibiotics and biocides as antimicrobial medicaments in endodontics

Bacteria have been implicated in the pathogenesis and progression of pulp and periapical diseases. The primary aim of endodontic treatment is to remove as many bacteria as possible from the root canal system and then to create an environment in which any remaining organisms cannot survive. This can only be achieved through the use of a combination of aseptic treatment techniques, chemomechanical preparation of the root canal, antimicrobial irrigating solutions and intracanal medicaments. The choice of which intracanal medicament to use is dependent on having an accurate diagnosis of the condition being treated, as well as a thorough knowledge of the type of organisms likely to be involved and. their mechanisms of growth and survival. Since the disease is likely to have been caused by the presence of bacteria within the root canal, the use of an antimicrobial agent is essential. Many medicaments have been used in an attempt to achieve the above aims, but no single preparation has been found to be completely predictable or effective. Commonly used medicaments include calcium hydroxide, antibiotics; non-phenolic biocides, phenolic biocides and iodine compounds. Each has advantages and disadvantages, and further research is required to determine which is best suited for root canal infections.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Outcome with calcium channel antagonists after myocardial infarction: A community-based study

Objectives. We sought to estimate the risk of death and recurrent myocardial infarction associated with the use of calcium antagonists after myocardial infarction in a population-based cohort study. Background. Calcium antagonists are commonly prescribed after myocardial infarction, but their long-term effects are not well established. Methods. Patients 25 to 69 years old with a suspected myocardial infarction were identified and followed up through a community-based register of myocardial infarction and cardiac death (part of the World Health Organization Monitoring Trends and Determinants in Cardiovascular Disease [MONICA] Project in Newcastle, Australia). Data were collected by review of medical records, in-hospital interview and review of death certificates. Results. From 1989 to 1993, 3,982 patients with a nonfatal suspected myocardial infarction were enrolled in the study. At hospital discharge, 1,001 patients were treated with beta-adrenergic blocking agents, 923 with calcium antagonists, 711 with both beta-blockers and calcium antagonists and 1,346 with neither drug. Compared with patients given beta-blockers, patients given calcium antagonists were more likely to suffer myocardial infarction or cardiac death (adjusted relative risk [RR] 1.4, 95% confidence interval [CI] 1.0 to 1.9), cardiac death (RR 1.6, 95% CI 1.0 to 2.7) and death from all causes (RR 1.7, 95% CI 1.1 to 2.6). Compared with patients given neither beta-blockers nor calcium antagonists, patients given calcium antagonists were not at increased risk of myocardial infarction or cardiac death (RR 1.0, 95% CI 0.8 to 1.3), cardiac death (RR 0.9, 95% CI 0.6 to 1.2) or death from all causes (RR 1.0, 95% CI 0.7 to 1.3). No excess in risk of myocardial infarction or cardiac death was observed among patients taking verapamil (RR 0.9, 95% CI 0.6 to 1.6), diltiazem (RR 1.1, 95% CI 0.8 to 1.4) or nifedipine (RR 1.3, 95% CI 0.7 to 2.2) compared,vith patients taking neither calcium antagonists nor beta-blockers. Conclusions. These results are consistent with randomized trial data showing benefit from beta blockers after myocardial infarction and no effect on the risk of recurrent myocardial infarction and death with the use of calcium antagonists. Comparisons between beta-blockers and calcium antagonists favor beta blockers because of the beneficial effects of beta-blockers and not because of adverse effects of calcium antagonists. (C) 1998 by the American College of Cardiology.

Characterization of calcium-activated chloride channels in patches excised from the dendritic knob of mammalian olfactory receptor neurons

We investigated the properties of calcium-activated chloride channels in inside-out membrane patches from the dendritic knobs of acutely dissociated rat olfactory receptor neurons. Patches typically contained large calcium-activated currents, with total conductances in the range 30-75 nS. The dose response curve for calcium exhibited an EC50 of about 26 mu M. In symmetrical NaCl solutions, the current-voltage relationship reversed at 0 mV and was linear between -80 and +70 mV. When the intracellular NaCl concentration was progressively reduced from 150 to 25 mM, the reversal potential changed in a manner consistent with a chloride-selective conductance. Indeed, modeling these data with the Goldman-Hodgkin-Katz equation revealed a P-Na/P-Cl of 0.034. The halide permeability sequence was P-Cl > P-F > P-I > P-Br indicating that permeation through the channel was dominated by ion binding sites with a high field strength. The channels were also permeable to the large organic anions, SCN-, acetate(-), and gluconate(-), with the permeability sequence P-Cl > P-SCN > gluconaie. Significant permeation to gluconate ions suggested that the channel pore had a minimum diameter of at least 5.8 Angstrom.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - II. Testing simulation capabilities of wheat growth

To simulate cropping systems, crop models must not only give reliable predictions of yield across a wide range of environmental conditions, they must also quantify water and nutrient use well, so that the status of the soil at maturity is a good representation of the starting conditions for the next cropping sequence. To assess the suitability for this task a range of crop models, currently used in Australia, were tested. The models differed in their design objectives, complexity and structure and were (i) tested on diverse, independent data sets from a wide range of environments and (ii) model components were further evaluated with one detailed data set from a semi-arid environment. All models were coded into the cropping systems shell APSIM, which provides a common soil water and nitrogen balance. Crop development was input, thus differences between simulations were caused entirely by difference in simulating crop growth. Under nitrogen non-limiting conditions between 73 and 85% of the observed kernel yield variation across environments was explained by the models. This ranged from 51 to 77% under varying nitrogen supply. Water and nitrogen effects on leaf area index were predicted poorly by all models resulting in erroneous predictions of dry matter accumulation and water use. When measured light interception was used as input, most models improved in their prediction of dry matter and yield. This test highlighted a range of compensating errors in all modelling approaches. Time course and final amount of water extraction was simulated well by two models, while others left up to 25% of potentially available soil water in the profile. Kernel nitrogen percentage was predicted poorly by all models due to its sensitivity to small dry matter changes. Yield and dry matter could be estimated adequately for a range of environmental conditions using the general concepts of radiation use efficiency and transpiration efficiency. However, leaf area and kernel nitrogen dynamics need to be improved to achieve better estimates of water and nitrogen use if such models are to be use to evaluate cropping systems. (C) 1998 Elsevier Science B.V.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.