Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, PHYSICA, for later application to full dynamic fluid structure interaction.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

Maximising heat transfer efficiency in the cold crucible induction melting process

Induction heating is an efficient method used to melt electrically conductive materials, particularly if melting takes place in a ceramic crucible. This form of melting is particularly good for alloys, as electromagnetic forces set up by the induction coil lead to vigorous stirring of the melt ensuring homogeneity and uniformity in temperature. However, for certain reactive alloys, or where high purity is required, ceramic crucibles cannot be used, but a water-cooled segmented copper crucible is employed instead. Water cooling prevents meltdown or distortion of the metal wall, but much of the energy goes into the coolant. To reduce this loss, the electromagnetic force generated by the coil is used to push the melt away from the walls and so minimise contact with water-cooled surfaces. Even then, heat is lost through the crucible base where contact is inevitable. In a collaborative programme between Greenwich and Birmingham Universities, computer modelling has been used in conjunction with experiments to improve the superheat attainable in the melt for a,number of alloys, especially for y-TiAl intermetallics to cast aeroengine turbine blades. The model solves the discretised form of the turbulent Navier-Stokes, thermal energy conservation and Maxwell equations using a Spectral Collocation technique. The time-varying melt envelope is followed explicitly during the computation using an adaptive mesh. This paper briefly describes the mathematical model used to represent the interaction between the magnetic field, fluid flow, heat transfer and change of phase in the crucible and identifies the proportions of energy used in the melt, lost in the crucible base and in the crucible walls. The role of turbulence is highlighted as important in controlling heat losses and turbulence damping is introduced as a means of improving superheat. Model validation is against experimental results and shows good agreement with measured temperatures and energy losses in the cooling fluid throughout the melting cycle.

Identifying the damage mechanisms in paper using the acoustic emission monitoring technique

This paper investigates the use of the acoustic emission (AE) monitoring technique for use in identifying the damage mechanisms present in paper associated with its production process. The microscopic structure of paper consists of a random mesh of paper fibres connected by hydrogen bonds. This implies the existence of two damage mechanisms, the failure of a fibre-fibre bond and the failure of a fibre. This paper describes a hybrid mathematical model which couples the mechanics of the mass-spring model to the acoustic wave propagation model for use in generating the acoustic signal emitted by complex structures of paper fibres under strain. The derivation of the mass-spring model can be found in [1,2], with details of the acoustic wave equation found in [3,4]. The numerical implementation of the vibro-acoustic model is discussed in detail with particular emphasis on the damping present in the numerical model. The hybrid model uses an implicit solver which intrinsically introduces artificial damping to the solution. The artificial damping is shown to affect the frequency response of the mass-spring model, therefore certain restrictions on the simulation time step must be enforced so that the model produces physically accurate results. The hybrid mathematical model is used to simulate small fibre networks to provide information on the acoustic response of each damage mechanism. The simulated AEs are then analysed using a continuous wavelet transform (CWT), described in [5], which provides a two dimensional time-frequency representation of the signal. The AEs from the two damage mechanisms show different characteristics in the CWT so that it is possible to define a fibre-fibre bond failure by the criteria listed below. The dominant frequency components of the AE must be at approximately 250 kHz or 750 kHz. The strongest frequency component may be at either approximately 250 kHz or 750 kHz. The duration of the frequency component at approximately 250 kHz is longer than that of the frequency component at approximately 750 kHz. Similarly, the criteria for identifying a fibre failure are given below. The dominant frequency component of the AE must be greater than 800 kHz. The duration of the dominant frequency component must be less than 5.00E-06 seconds. The dominant frequency component must be present at the front of the AE. Essentially, the failure of a fibre-fibre bond produces a low frequency wave and the failure of a fibre produces a high frequency pulse. Using this theoretical criteria, it is now possible to train an intelligent classifier such as the Self-Organising Map (SOM) [6] using the experimental data. First certain features must be extracted from the CWTs of the AEs for use in training the SOM. For this work, each CWT is divided into 200 windows of 5E-06s in duration covering a 100 kHz frequency range. The power ratio for each windows is then calculated and used as a feature. Having extracted the features from the AEs, the SOM can now be trained, but care is required so that the both damage mechanisms are adequately represented in the training set. This is an issue with paper as the failure of the fibre-fibre bonds is the prevalent damage mechanism. Once a suitable training set is found, the SOM can be trained and its performance analysed. For the SOM described in this work, there is a good chance that it will correctly classify the experimental AEs.

Droplet oscillations in high gradient magnetic field

In high intensity and high gradient magnetic fields the volumetric force on diamagnetic material, such as water, leads to conditions very similar to microgravity in a terrestrial laboratory. In principle, this opens the possibility to determine material properties of liquid samples without wall contact, even for electrically non-conducting materials. In contrast, AC field levitation is used for conductors, but then terrestrial conditions lead to turbulent flow driven by Lorentz forces. DC field damping of the flow is feasible and indeed practiced to allow property measurements. However, the AC/DC field combination acts preferentially on certain oscillation modes and leads to a shift in the droplet oscillation spectrum.What is the cause? A nonlinear spectral numerical model is presented, to address these problems

Turbulence model performance in continuous casting simulations

This paper concerns the development and validation (using an oil/water system) of a finite volume computer model of the continuous casting process for steel flat products. The emphasis is on hydrodynamic aspects and in particular the dynamic behaviour of the metal/slag interface. Instability and wave action encourage the entrainment of inclusions into the melt affecting product quality. To track the interface between oil and water a new implicit algorithm was developed, called the Counter Diffusion Method. To prevent excessive damping, a time-filtered version of the k-e model, was found necessary, with appropriate density stratification terms representing interface turbulence damping.

Droplet oscillations in high gradient static magnetic field

In high intensity and high gradient magnetic fields the volumetric force on diamagnetic material, such as water, leads to conditions very similar to microgravity in a terrestrial laboratory. In principle, this opens the possibility to determine material properties of liquid samples without wall contact, even for electrically non-conducting materials. In contrast, AC field levitation is used for conductors, but then terrestrial conditions lead to turbulent flow driven by Lorentz forces. DC field damping of the flow is feasible and indeed practiced to allow property measurements. However, the AC/DC field combination acts preferentially on certain oscillation modes and leads to a shift in the droplet oscillation spectrum.What is the cause? A nonlinear spectral numerical model is presented, to address these problems.

Completely positive non-Markovian decoherence

We propose an effective Hamiltonian approach to investigate decoherence of a quantum system in a non-Markovian reservoir, naturally imposing the complete positivity on the reduced dynamics of the system. The formalism is based on the notion of an effective reservoir, i.e., certain collective degrees of freedom in the reservoir that are responsible for the decoherence. As examples for completely positive decoherence, we present three typical decoherence processes for a qubit such as dephasing, depolarizing, and amplitude damping. The effects of the non-Markovian decoherence are compared to the Markovian decoherence.

Eclipse observations of high-frequency oscillations in active region coronal loops

One of the mechanisms proposed for heating the corona above solar active regions is the damping of magnetohydrodynamic (MHD) waves. Continuing on previous work, we provide observational evidence for the existence of high-frequency MHD waves in coronal loops observed during the August 1999 total solar eclipse. A wavelet analysis is used to identify twenty 4 x 4 arcsec(2) areas showing intensity oscillations. All detections lie in the frequency range 0.15 - 0.25 Hz (7 - 4 s), last for at least 3 periods at a confidence level of more than 99% and arise just outside known coronal loops. This leads us to suggest that they occur in low emission-measure or different temperature loops associated with the active region.

Surface plasmon polariton propagation length: A direct comparison using photon scanning tunneling microscopy and attenuated total reflection

The propagation of surface plasmon polaritons (SPP's) is studied using a photon scanning tunneling microscope (PSTM) and conventional attenuated total reflection (ATR). The PSTM experiment uses localized (focused beam) launching or SPP's at a wavelength of 632.8 nm. Propagation of the SPP is observed as an exponentially decaying tail beyond the launch site acid the 1/e propagation length is measured directly for a series of Ag films of different thicknesses. The ATR measurements are used to characterize the thin film optical and thickness parameters, revealing, notably, the presence of a contaminating adlayer of Ag2S of typical dielectric function, 8.7 + i2.7, and thickness 1-2 nm. Values of the SPP propagation length, based on the ATR- derived film parameters used in the four-media implicit SPP dispersion relation, show very good agreement with those based on the PSTM images for the case of undercoupled or optimally coupled SPP modes. The observed propagation lengths are quantitatively analyzed taking explicit account of additional intrinsic damping due to the growth of the Ag2S layer and of reradiation of the SPP back into the prism outside the launch site. Finally, the PSTM images show excellent SPP beam confinement in the original propagation direction.

Comparison of LIMPET contra-rotating wells turbine with theoretical and model test predictions

The performance of the contra-rotating Wells turbine installed in the LIMPET wave power station is compared to the predicted performance from theoretical analysis and model tests. A reasonable agreement was found between the predicted and measured turbine damping characteristic, however the turbine efficiency was found to be poorly predicted. It is postulated that this is due to the unsteady nature and mal-distribution of flow through the LIMPET turbine, which were not considered in the predictions. It is suggested that the reduced performance of the contra-rotating Wells turbine makes biplane or monoplane Wells turbines with guide vanes better solutions for OWC's.

Experimental demonstration of decoherence-free one-way information transfer

We report the experimental demonstration of a one-way quantum protocol reliably operating in the presence of decoherence. Information is protected by designing an appropriate decoherence-free subspace for a cluster state resource. We demonstrate our scheme in an all-optical setup, encoding the information into the polarization states of four photons. A measurement-based one-way information-transfer protocol is performed with the photons exposed to severe symmetric phase-damping noise. Remarkable protection of information is accomplished, delivering nearly ideal outcomes.

One-way quantum computing in a decoherence-free subspace

We introduce a novel scheme for one-way quantum computing (QC) based on the use of information encoded qubits in an effective cluster state resource. With the correct encoding structure, we show that it is possible to protect the entangled resource from phase damping decoherence, where the effective cluster state can be described as residing in a decoherence-free subspace (DFS) of its supporting quantum system. One-way QC then requires either single or two-qubit adaptive measurements. As an example where this proposal can be realized, we describe an optical lattice set-up where the scheme provides robust quantum information processing. We also outline how one can adapt the model to provide protection from other types of decoherence.

Discovery of Spatial Periodicities in a Coronal Loop Using Automated Edge-Tracking Algorithms

A new method for automated coronal loop tracking, in both spatial and temporal domains, is presented. Applying this technique to TRACE data, obtained using the 171 angstrom filter on 1998 July 14, we detect a coronal loop undergoing a 270 s kink-mode oscillation, as previously found by Aschwanden et al. However, we also detect flare-induced, and previously unnoticed, spatial periodicities on a scale of 3500 km, which occur along the coronal loop edge. Furthermore, we establish a reduction in oscillatory power for these spatial periodicities of 45% over a 222 s interval. We relate the reduction in detected oscillatory power to the physical damping of these loop-top oscillations.

Anomalous helicon wave absorption and parametric excitation of electrostatic fluctuations in a helicon-produced plasma

The nonlinear nature of the rf absorption in a helicon-produced plasma was recently evidenced by the observation that the helicon wave damping as well as the level of short-scale electrostatic fluctuations excited in the helicon plasma increases with rf power. Correlation methods using electrostatic probes as well as microwave back-scattering at the upper-hybrid resonance allow identifying the fluctuations as ion-sound and Trivelpiece– Gould waves satisfying the frequency and wavenumber matching conditions for the parametric decay instability of the helicon pump wave. Furthermore, the growth rates and thresholds deduced from their temporal growth are in good agreement with theoretical predictions for the parametric decay instability that takes into account realistic damping rates for the decay waves as well as a non-vanishing parallel wavenumber of the helicon pump. The close relationship between the rf absorption and the excitation of the fluctuations was investigated in more detail by performing time- and space-resolved measurements of the helicon wave field and the electrostatic fluctuations.