Magnetic damping of levitated liquid droplets in AC and DC field

The intense AC magnetic field required to produce levitation in terrestrial conditions, along with the buoyancy and thermo-capillary forces, results in turbulent convective flow within the droplet. The use of a homogenous DC magnetic field allows the convective flow to be damped. However the turbulence properties are affected at the same time, leading to a possibility that the effective turbulent damping is considerably reduced. The MHD modified K-Omega turbulence model allows the investigation of the effect of magnetic field on the turbulence. The model incorporates free surface deformation, the temperature dependent surface tension, turbulent momentum transport, electromagnetic and gravity forces. The model is adapted to incorporate a periodic laser heating at the top of the droplet, which have been used to measure the thermal conductivity of the material by calculating the phase lag between the frequency of the laser heating and the temperature response at the bottom. The numerical simulations show that with the gradual increase of the DC field the fluid flow within the droplet is initially increasing in intensity. Only after a certain threshold magnitude of the field the flow intensity starts to decrease. In order to achieve the flow conditions close to the ‘laminar’ a D.C. magnetic field >4 Tesla is required to measure the thermal conductivity accurately. The reduction in the AC field driven flow in the main body of the drop leads to a noticeable thermo-capillary convection at the edge of the droplet. The uniform vertical DC magnetic field does not stop a translational oscillation of the droplet along the field, which is caused by the variation in total levitation force due to the time-dependent surface deformation.

Singular stresses at the tip of a crack terminating at frictional interface of fibre reinforced composites

The stress singularities at the tip of a crack that terminates at a frictional interface between two layers in anisotropic composites are investigated. The order of stress singularities is determined by solving the characteristic equations obtained from the boundary conditions and the frictional interface conditions for the cases concerned. The interface is assumed to be governed by Coulomb's law of friction. Numerical results are presented for the cases with a crack terminating at a frictional interface of a fibre reinforced composite, and it is shown that there is a big difference of stress singularities between cases with and without considering friction along the interface.

An experimental and numerical CFD study of turbulence in a tundish container

This paper describes work performed at IRSID/USINOR in France and the University of Greenwich, UK, to investigate flow structures and turbulence in a water-model container, simulating aspects typical of metal tundish operation. Extensive mean and fluctuating velocity measurements were performed at IRSID using LDA to determine the flow field and these form the basis for a numerical model validation. This apparently simple problem poses several difficulties for the CFD modelling. The flow is driven by the strong impinging jet at the inlet. Accurate description of the jet is most important and requires a localized fine grid, but also a turbulence model that predicts the correct spreading rates of jet and impinging wall boundary layers. The velocities in the bulk of the tundish tend to be (indeed need to be) much smaller than those of the jet, leading to damping of turbulence, or even laminar flow. The authors have developed several low-Reynolds number (low-Re) k–var epsilon model variants to compute this flow and compare against measurements. Best agreement is obtained when turbulence damping is introduced to account not only for walls, but also for low-Re regions in the bulk – the k–var epsilon model otherwise allows turbulence to accumulate in the container due to the restricted outlet. Several damping functions are tested and the results reported here. The k–ω model, which is more suited to transitional flow, also seems to perform well in this problem.

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.