A partial velocity approach to subcycling structural dynamics

Subcycling, or the use of different timesteps at different nodes, can be an effective way of improving the computational efficiency of explicit transient dynamic structural solutions. The method that has been most widely adopted uses a nodal partition. extending the central difference method, in which small timestep updates are performed interpolating on the displacement at neighbouring large timestep nodes. This approach leads to narrow bands of unstable timesteps or statistical stability. It also can be in error due to lack of momentum conservation on the timestep interface. The author has previously proposed energy conserving algorithms that avoid the first problem of statistical stability. However, these sacrifice accuracy to achieve stability. An approach to conserve momentum on an element interface by adding partial velocities is considered here. Applied to extend the central difference method. this approach is simple. and has accuracy advantages. The method can be programmed by summing impulses of internal forces, evaluated using local element timesteps, in order to predict a velocity change at a node. However, it is still only statistically stable, so an adaptive timestep size is needed to monitor accuracy and to be adjusted if necessary. By replacing the central difference method with the explicit generalized alpha method. it is possible to gain stability by dissipating the high frequency response that leads to stability problems. However. coding the algorithm is less elegant, as the response depends on previous partial accelerations. Extension to implicit integration, is shown to be impractical due to the neglect of remote effects of internal forces acting across a timestep interface. (C) 2002 Elsevier Science B.V. All rights reserved.

Microstructural evolution of the Serido Belt, NE-Brazil: the effect of two tectonic events on development of c-axis preferred orientation in quartz

The polyphase evolution of the Serido Belt (NE-Brazil) includes D, crust formation at 2.3-2.1 Ga, D-2 thrust tectonics at 1.9 Ga and crustal reworking by D-3 strike-slip shear zones at 600 Ma. Microstructural investigations within mylonites associated with D-2 and D-3 events were used to constrain the tectono-thermal evolution of the belt. D-2 shear zones commenced at deeper crustal levels and high amphibolite facies conditions (600-650 degreesC) through grain boundary migration, subgrain rotation and operation of quartz Q-prism slip. Continued shearing and exhumation of the terrain forced the re-equilibration of high-T fabrics and the switching of slip systems from (c)-prism to positive and negative (a)-rhombs. During D-3, enhancement of ductility by dissipation of heat that came from syn-D-3 granites developed wide belts of amphibolite facies mylonites. Continued shearing, uplift and cooling of the region induced D-3 shear zones to act in ductile-brittle regimes, marked by fracturing and development of thinner belts of greenschist facies mylonites. During this event, switching from (a)-prism to a basal slip indicates a thermal path from 600 to 350 degreesC. Therefore, microstructures and quartz c-axis fabrics in polydeformed rocks from the Serido Belt preserve the record of two major events, which includes contrasting deformation mechanisms and thermal paths. (C) 2003 Elsevier Ltd. All rights reserved.

New fluorescence parameters for monitoring photosynthesis in plants - 1. The effect of illumination on the fluorescence parameters of the JIP-test

Chlorophyll fluorescence measurements have a wide range of applications from basic understanding of photosynthesis functioning to plant environmental stress responses and direct assessments of plant health. The measured signal is the fluorescence intensity (expressed in relative units) and the most meaningful data are derived from the time dependent increase in fluorescence intensity achieved upon application of continuous bright light to a previously dark adapted sample. The fluorescence response changes over time and is termed the Kautsky curve or chlorophyll fluorescence transient. Recently, Strasser and Strasser (1995) formulated a group of fluorescence parameters, called the JIP-test, that quantify the stepwise flow of energy through Photosystem II, using input data from the fluorescence transient. The purpose of this study was to establish relationships between the biochemical reactions occurring in PS II and specific JIP-test parameters. This was approached using isolated systems that facilitated the addition of modifying agents, a PS II electron transport inhibitor, an electron acceptor and an uncoupler, whose effects on PS II activity are well documented in the literature. The alteration to PS II activity caused by each of these compounds could then be monitored through the JIP-test parameters and compared and contrasted with the literature. The known alteration in PS II activity of Chenopodium album atrazine resistant and sensitive biotypes was also used to gauge the effectiveness and sensitivity of the JIP-test. The information gained from the in vitro study was successfully applied to an in situ study. This is the first in a series of four papers. It shows that the trapping parameters of the JIP-test were most affected by illumination and that the reduction in trapping had a run-on effect to inhibit electron transport. When irradiance exposure proceeded to photoinhibition, the electron transport probability parameter was greatly reduced and dissipation significantly increased. These results illustrate the advantage of monitoring a number of fluorescence parameters over the use of just one, which is often the case when the F-V/F-M ratio is used.

Scaling of decoherence effects in quantum computers

The scaling of decoherence rates with qubit number N is studied for a simple model of a quantum computer in the situation where N is large. The two state qubits are localized around well-separated positions via trapping potentials and vibrational centre of mass motion of the qubits occurs. Coherent one and two qubit gating processes are controlled by external classical fields and facilitated by a cavity mode ancilla. Decoherence due to qubit coupling to a bath of spontaneous modes, cavity decay modes and to the vibrational modes is treated. A non-Markovian treatment of the short time behaviour of the fidelity is presented, and expressions for the characteristic decoherence time scales obtained for the case where the qubit/cavity mode ancilla is in a pure state and the baths are in thermal states. Specific results are given for the case where the cavity mode is in the vacuum state and gating processes are absent and the qubits are in (a) the Hadamard state (b) the GHZ state.

Entanglement induced by spontaneous emission in spatially extended two-atom systems

The role of the collective antisymmetric state in entanglement creation by spontaneous emission in a system of two non-overlapping two-level atoms has been investigated. Populations of the collective atomic states and the Wootters entanglement measure (concurrence) for two sets of initial atomic conditions are calculated and illustrated graphically. Calculations include the dipole-dipole interaction and a spatial separation between the atoms that the antisymmetric state of the system is included throughout even for small interatomic separations. It is shown that spontaneous emission can lead to a transient entanglement between the atoms even if the atoms were prepared initially in an unentangled state. It is found that the ability of spontaneous emission to create transient entanglement relies on the absence of population in the collective symmetric state of the system. For the initial state of only one atom excited, entanglement builds up rapidly in time and reaches a maximum for parameter values corresponding roughly to zero population in the symmetric state. On the other hand, for the initial condition of both atoms excited, the atoms remain unentangled until the symmetric state is depopulated. A simple physical interpretation of these results is given in terms of the diagonal states of the density matrix of the system. We also study entanglement creation in a system of two non-identical atoms of different transition frequencies. It is found that the entanglement between the atoms can be enhanced compared to that for identical atoms, and can decay with two different time scales resulting from the coherent transfer of the population from the symmetric to the antisymmetric state. In addition, it was found that a decaying initial entanglement between the atoms can display a revival behaviour.

Stationary two-atom entanglement induced by nonclassical two-photon correlations

A system of two two-level atoms interacting with a squeezed vacuum field can exhibit stationary entanglement associated with nonclassical two-photon correlations characteristic of the squeezed vacuum field. The amount of entanglement present in the system is quantified by the well known measure of entanglement called concurrence. We find analytical formulae describing the concurrence for two identical and nonidentical atoms and show that it is possible to obtain a large degree of steady-state entanglement in the system. Necessary conditions for the entanglement are nonclassical two-photon correlations and nonzero collective decay. It is shown that nonidentical atoms are a better source of stationary entanglement than identical atoms. We discuss the optimal physical conditions for creating entanglement in the system; in particular, it is shown that there is an optimal and rather small value of the mean photon number required for creating entanglement.

Entangling two atoms via spontaneous emission

We discuss the creation of entanglement between two two-level atoms in the dissipative process of spontaneous emission. It is shown that spontaneous emission can lead to a transient entanglement between the atoms even if the atoms were prepared initially in an unentangled state. The amount of entanglement created in the system is quantified by using two different measures: concurrence and negativity. We find analytical formulae for the evolution of concurrence and negativity in the system. We also find the analytical relation between the two measures of entanglement. The system consists of two two-level atoms which are separated by an arbitrary distance r(12) and interact with each other via the dipole-dipole interaction, and the antisymmetric state of the system is included throughout, even for small interatomic separations, in contrast to the small-sample model. It is shown that for sufficiently large values of the dipole-dipole interaction initially the entanglement exhibits oscillatory behaviour with considerable entanglement in the peaks. For longer times the amount of entanglement is directly related to the population of the slowly decaying antisymmetric state.

Merging of low-mass systems and the origin of the Fundamental Plane

We present a new set of dissipationless N-body simulations to examine the feasibility of creating bright ellipticals (following the Kormendy relation, hereafter KR) by hierarchically merging present-day early-type dwarf galaxies, and to study how the encounter parameters affect the location of the end product in the (mu(e))-R-e plane. We investigate the merging of one-component galaxies of both equal and different masses, the merging of two-component galaxy models to explore the effect of dark haloes on the final galaxy characteristics, and the merging of ultracompact dwarf galaxies. We find that the increase of (mu(e)) with R-e is attributable to an increase in the initial orbital energy. The merger remnants shift down in the (mu(e))-R-e plane and fail to reach the KR. Thus, the KR is not reproducible by mergers of dwarf early-type systems, rendering untenable the theory that present-day dwarfs are responsible for even a small fraction of the present-day ellipticals, unless a considerable amount of dissipation is invoked. However, we do find that present-day dwarfs can be formed by the merger of ultracompact dwarfs.

Disaggregation of Microcystis aeruginosa colonies under turbulent mixing: laboratory experiments in a grid-stirred tank

Samples of the cyanobacterium Microcystis aeruginosa from a small pond were used in laboratory experiments with a grid-stirred tank to quantify the effect of turbulent mixing on colony size. Turbulent dissipation in the tank was varied from 10(-9) m(2) s(-3) to 10(-4) m(2) s(-3), covering the range of turbulence intensities experienced by M. aeruginosa colonies in the field and exceeding the maximum dissipation by two orders of magnitude. Large colonies broke up into smaller colonies during the experiments; the mass fraction of colonies with diameter less than 200 mum increased over time. Colony disaggregation was observed to increase with turbulent dissipation. The maximum stable colony diameter across all experiments was in the range 220-420 mum. The overall change in size distribution during the experiments was relatively small, and the colony size distribution remained very broad throughout the experiments. Since colony size affects migration velocity, susceptibility to grazing and surface area to volume ratios, more work is needed to determine how to best represent this broad size distribution when modelling M. aeruginosa populations.

Undular tidal bore dynamics in the Daly Estuary, Northern Australia

Measurements in the macro-tidal Daly Estuary show that the presence of an undular tidal bore contributed negligibly to the dissipation of tidal energy. No recirculation bubble was observed between a trough and the following wave crest in the lee waves following the undular bore. This differs to stationary undular bores in laboratory experiments at larger Froude numbers where a recirculation bubble exists. Secondary motions and the turbulence generated by the undular bore had no measurable influence on the sediment transport. This situation contrasts with the intense sediment resuspension observed in breaking tidal bores. The tidally averaged sediment budget in the Daly Estuary was controlled by the asymmetry of tidal currents. The undular bore may widen the river by breaking along the banks that it undercuts, leading to bank slippage. A patch of river-wide macro-turbulence of 3-min duration occurred about 20 min after the passage of the bore during accelerating tidal currents. (C) 2004 Elsevier Ltd. All rights reserved.

Accelerating precursory activity within a class of earthquake analogue automata

A statistical fractal automaton model is described which displays two modes of dynamical behaviour. The first mode, termed recurrent criticality, is characterised by quasi-periodic, characteristic events that are preceded by accelerating precursory activity. The second mode is more reminiscent of SOC automata in which large events are not preceded by an acceleration in activity. Extending upon previous studies of statistical fractal automata, a redistribution law is introduced which incorporates two model parameters: a dissipation factor and a stress transfer ratio. Results from a parameter space investigation indicate that a straight line through parameter space marks a transition from recurrent criticality to unpredictable dynamics. Recurrent criticality only occurs for models within one corner of the parameter space. The location of the transition displays a simple dependence upon the fractal correlation dimension of the cell strength distribution. Analysis of stress field evolution indicates that recurrent criticality occurs in models with significant long-range stress correlations. A constant rate of activity is associated with a decorrelated stress field.

Using continuous measurement to protect a universal set of quantum gates within a perturbed decoherence-free subspace

We consider a universal set of quantum gates encoded within a perturbed decoherence-free subspace of four physical qubits. Using second-order perturbation theory and a measuring device modelled by an infinite set of harmonic oscillators, simply coupled to the system, we show that continuous observation of the coupling agent induces inhibition of the decoherence due to spurious perturbations. We thus advance the idea of protecting or even creating a decoherence-free subspace for processing quantum information.

Matching conditional moments in PDF modelling of nonpremixed combustion

The rate of generation of fluctuations with respect to the scalar values conditioned on the mixture fraction, which significantly affects turbulent nonpremixed combustion processes, is examined. Simulation of the rate in a major mixing model is investigated and the derived equations can assist in selecting the model parameters so that the level of conditional fluctuations is better reproduced by the models. A more general formulation of the multiple mapping conditioning (MMC) model that distinguishes the reference and conditioning variables is suggested. This formulation can be viewed as a methodology of enforcing certain desired conditional properties onto conventional mixing models. Examples of constructing consistent MMC models with dissipation and velocity conditioning as well as of combining MMC with large eddy simulations (LES) are also provided. (c) 2005 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Prediction of chaotic instabilities in a dragline bucket swing

The occurrence of chaotic instabilities is investigated in the swing motion of a dragline bucket during operation cycles. A dragline is a large, powerful, rotating multibody system utilised in the mining industry for removal of overburden. A simplified representative model of the dragline is developed in the form of a fundamental non-linear rotating multibody system with energy dissipation. An analytical predictive criterion for the onset of chaotic instability is then obtained in terms of critical system parameters using Melnikov's method. The model is shown to exhibit chaotic instability due to a harmonic slew torque for a range of amplitudes and frequencies. These chaotic instabilities could introduce irregularities into the motion of the dragline system, rendering the system difficult to control by the operator and/or would have undesirable effects on dragline productivity and fatigue lifetime. The sufficient analytical criterion for the onset of chaotic instability is shown to be a useful predictor of the phenomenon under steady and unsteady slewing conditions via comparisons with numerical results. (c) 2005 Elsevier Ltd. All rights reserved.

On the modelling of wave breaking and set-up on coral reefs

An extended refraction-diffraction equation [Massel, S.R., 1993. Extended refraction-diffraction equation for surface waves. Coastal Eng. 19, 97-126] has been applied to predict wave transformation and breaking as well as wave-induced set-up on two-dimensional reef profiles of various shapes. A free empirical coefficient alpha in a formula for the average rate of energy dissipation [epsilon(b)] = (alpha rho g omega/8 pi)(root gh/C)(H-3/h) in the modified periodic bore model was found to be a function of the dimensionless parameter F-c0 = (g(1.25)H(0)(0.5)T(2.5))/h(r)(1.75), proposed by Gourlay [Gourlayl M.R., 1994. Wave transformation on a coral reef. Coastal Eng. 23, 17-42]. The applicability of the developed model has been demonstrated for reefs of various shapes subjected to various incident wave conditions. Assuming proposed relationships of the coefficient alpha and F-c0, the model provides results on wave height attenuation and set-up elevation which compare well with experimental data. (C) 2000 Elsevier Science B.V. All rights reserved.