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.

Single mode Lamb waves in composite laminated plates generated by piezoelectric transducers

The technique of permanently attaching piezoelectric transducers to structural surfaces has demonstrated great potential for quantitative non-destructive evaluation and smart materials design. For thin structural members such as composite laminated plates, it has been well recognized that guided Lamb wave techniques can provide a very sensitive and effective means for large area interrogation. However, since in these applications multiple wave modes are generally generated and the individual modes are usually dispersive, the received signals are very complex and difficult to interpret. An attractive way to deal with this problem has recently been introduced by applying piezoceramic transducer arrays or interdigital transducer (IDT) technologies. In this paper, the acoustic wave field in composite laminated plates excited by piezoceramic transducer arrays or IDT is investigated. Based on dynamic piezoelectricity theory, a discrete layer theory and a multiple integral transform method, an analytical-numerical approach is developed to evaluate the input impedance characteristics of the transducer and the surface velocity response of the plate. The method enables the quantitative evaluation of the influence of the electrical characteristics of the excitation circuit, the geometric and piezoelectric properties of the transducer array, and the mechanical and geometrical features of the laminate. Numerical results are presented to validate the developed method and show the ability of single wave mode selection and isolation. The results show that the interaction between individual elements of the piezoelectric array has a significant influence on the performance of the IDT, and these effects can not be neglected even in the case of low frequency excitation. It is also demonstrated that adding backing materials to the transducer elements can be used to improve the excitability of specific wave modes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Hemodynamic and ventilatory effects of manual respiratory physiotherapy techniques of chest clapping, vibration, and shaking in an animal model

Chest clapping, vibration, and shaking were studied in 10 physiotherapists who applied these techniques on an anesthetized animal model. Hemodynamic variables (such as heart rate, blood pressure, pulmonary artery pressure, and right atrial pressure) were measured during the application of these techniques to verify claims of adverse events. In addition, expired tidal volume and peak expiratory flow rate were measured to ascertain effects of these techniques. Physiotherapists in this study applied chest clapping at a rate of 6.2 +/- 0.9 Hz, vibration at 10.5 +/- 2.3 Hz, and shaking at 6.2 +/- 2.3 Hz. With the use of these rates, esophageal pressure swings of 8.8 +/- 5.0, 0.7 +/- 0.3, and 1.4 +/- 0.7 mmHg resulted from clapping, vibration, and shaking respectively. Variability in rates and forces generated by these techniques was 80% of variance in shaking force (P = 0.003). Application of these techniques by physiotherapists was found to have no significant effects on hemodynamic and most ventilatory variables in this study. From this study, we conclude that chest clapping, vibration, and shaking 1) can be consistently performed by physiotherapists; 2) are significantly related to physiotherapists' characteristics, particularly clinical experience; and 3) caused no significant hemodynamic effects.

Numerical investigation of turbulent flow around a rotating stepped cylinder for corrosion study

Direct numerical simulation has been carried out for turbulent flow set up by a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction. This flow geometry creates a qualitatively similar flow pattern as observed near, a sudden, pipe expansion or a plane backward-facing step, characterized by flow separation and reattachment. A region of intense turbulence intensity and high wall-shear-stress fluctuations is formed in,the recirculating I region downstream of the step, where high mass-transfer capacity was also experimentally observed. Since, corrosion is frequently mass-transfer., controlled, our findings, put forward this apparatus as useful tool for future corrosion research.

Numerical study of hydrogenic effective mass theory for an impurity P donor in Si in the presence of an electric field and interfaces

In this paper we examine the effects of varying several experimental parameters in the Kane quantum computer architecture: A-gate voltage, the qubit depth below the silicon oxide barrier, and the back gate depth to explore how these variables affect the electron density of the donor electron. In particular, we calculate the resonance frequency of the donor nuclei as a function of these parameters. To do this we calculated the donor electron wave function variationally using an effective-mass Hamiltonian approach, using a basis of deformed hydrogenic orbitals. This approach was then extended to include the electric-field Hamiltonian and the silicon host geometry. We found that the phosphorous donor electron wave function was very sensitive to all the experimental variables studied in our work, and thus to optimize the operation of these devices it is necessary to control all parameters varied in this paper.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.

Usefulness of quantitative echocardiographic techniques to predict recovery of regional and global left ventricular function after acute myocardial infarction

The left ventricular response to dobutamine may be quantified using tissue Doppler measurement of myocardial velocity or displacement or 3-dimensional echocardiography to measure ventricular volume and ejection fraction. This study sought to explore the accuracy of these methods for predicting segmental and global responses to therapy. Standard dobutamine and 3-dimensional echocardiography were performed in 92 consecutive patients with abnormal left ventricular function at rest. Recovery of function was defined by comparison with follow-up echocardiography at rest 5 months later. Segments that showed improved regional function at follow-up showed a higher increment in peak tissue Doppler velocity with dobutamine therapy than in nonviable segments (1.2 +/- 0.4 vs 0.3 +/- 0.2 cm/s, p = 0.001). Similarly, patients who showed a > 5% improvement of ejection fraction at follow-up showed a greater displacement response to dobutamine (6.9 +/- 3.2 vs 2.1 +/- 2.3 mm, p = 0.001), as well as a higher rate of ejection fraction, response to dobutamine (9 +/- 3% vs 2 +/- 2%, p = 0.001). The optimal cutoff values for predicting subsequent recovery of function at rest were an increment of peak velocity > 1 cm/s, >5 mm of displacement, and a >5% improvement of ejection fraction with low-dose dobutamine. (C) 2003 by Excerpta Medica, Inc.