Monotone scheme and boundary conditions for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number

A monotone scheme for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number is presented. The numerical stability is analysed with respect to the electromagnetic force. Standard central finite differences applied to finite volumes can only be numerically stable if the vector products involved in this force are computed with a scheme using a fully staggered grid. The electromagnetic quantities (electric currents and electric potential) must be shifted by half the grid size from the mechanical ones (velocity and pressure). An integral treatment of the boundary layers is used in conjunction with boundary conditions for electrically conducting walls. The simulations are performed with inhomogeneous electrical conductivities of the walls and reach high Hartmann numbers in three-dimensional simulations, even though a non-adaptive grid is used.

Two fluid approach to high speed impact interaction amongst solid structures

The issues surrounding collision of projectiles with structures has gained a high profile since the events of 11th September 2001. In such collision problems, the projectile penetrates the stucture so that tracking the interface between one material and another becomes very complex, especially if the projectile is essentially a vessel containing a fluid, e.g. fuel load. The subsequent combustion, heat transfer and melting and re-solidification process in the structure render this a very challenging computational modelling problem. The conventional approaches to the analysis of collision processes involves a Lagrangian-Lagrangian contact driven methodology. This approach suffers from a number of disadvantages in its implementation, most of which are associated with the challenges of the contact analysis component of the calculations. This paper describes a 'two fluid' approach to high speed impact between solid structures, where the objective is to overcome the problems of penetration and re-meshing. The work has been carried out using the finite volume, unstructured mesh multi-physics code PHYSICA+, where the three dimensional fluid flow, free surface, heat transfer, combustion, melting and re-solidification algorithms are approximated using cell-centred finite volume, unstructured mesh techniques on a collocated mesh. The basic procedure is illustrated for two cases of Newtonian and non-Newtonian flow to test various of its component capabilities in the analysis of problems of industrial interest.

Water transport in zinc oxy-chloride cements

The water uptake and water loss behaviour in three different formulations of zinc oxy-chloride cement have been studied in detail. Specimens of each material were subjected to a high humidity atmosphere (93% RH) over saturated aqueous sodium sulfate, and a low humidity desiccating atmosphere over concentrated sulfuric acid. In high humidity, the cement formulated from the nominal 75% ZnCl2 solutions gained mass, eventually becoming too sticky to weigh further. The specimens at 25% and 50% ZnCl2 by contrast lost mass by a diffusion process, though by 1 week the 50% cement had stated to gain mass and was also too sticky to weigh. In low humidity, all three cements lost mass, again by a diffusion process. Both water gain and water loss followed Fick's law for a considerable time. In the case of water loss under desiccating conditions, this corresponded to values of Mt/MĄ well above 0.5. However, plots did not go through the origin, showing that there was an induction period before true diffusion began. Diffusion coefficients varied from 1.56 x 10-5 (75% ZnCl2) to 2.75 x 10-5 cm2/s (50% ZnCl2), and appeared to be influenced not simply by composition. The drying of the 25% and 50% ZnCl2 cements in high humidity conditions occurred at a much lower rate, with a value of D of 2.5 x 10-8 cm2/s for the 25% ZnCl2 cement. This cement was found to equilibrate slowly, but total water loss did not differ significantly from that of the cements stored under desiccating conditions. Equilibration times for water loss in desiccating conditions were of the order of 2-4 hours, depending on ZnCl2 content; equilibrium water losses were respectively 28.8 [25% ZnCl2], 16.2 [50%] and 12.4 [75%] which followed the order of ZnCl2 content. It is concluded that the water transport processes are strongly influenced by the ZnCl2 content of the cement.

Solid phase velocity measurements utilising electrostatic sensors and cross correlation signal processing

Gas-solids two phase systems are widely employed within process plant in the form of pneumatic conveyors, dust extraction systems and solid fuel injection systems. The measurement of solids phase velocity therefore has wide potential application in flow monitoring and, in conjunction with density measurement instrumentation, solids mass flow rate measurement. Historically, a number of authors have detailed possible measurement techniques, and some have published limited test results. It is, however, apparent that none of these technologies have found wide application in industry. Solids phase velocity measurements were undertaken using real time cross correlation of signals from two electrostatic sensors spaced axially along a pipeline conveying pulverised coal (PF). Details of the measurement equipment, the pilot scale test rig and the test results are presented.

Laser Doppler velocimetry measurements of particle velocity profiles in gas-solid two-phase flows

The measurement of particle velocities in two-phase gas-solid systems has a wide application in flow monitoring in process plant, where two-phase gas-solids systems are frequently employed in the form of pneumatic conveyors and solid fuel injection systems. Such measurements have proved to be difficult to make reliably in industrial environments. This paper details particle velocity measurements made in a two phase gas-solid now utilising a laser Doppler velocimetry system. Tests were carried out using both wheat flour and pulverised coal as the solids phase, with air being used as the gaseous phase throughout. A pipeline of circular section, having a diameter of 53 mm was used for the test work, with air velocities ranging from 25 to 45 m/s and suspension densities ranging from 0.001 kg to 1 kg of solids per cubic meter of air. Details of both the test equipment used, and the results of the measurements are presented.

Velocity measurement of pneumatically conveyed solids using electrodynamic sensors

Theoretical and experimental studies of cross correlation techniques applied to non-restrictive velocity measurement of pneumatically conveyed solids using ring-shaped electrodynamic flow sensors are presented. In-depth studies of the electrodynamic sensing mechanism, and also of the spatial sensitivity and spatial filtering properties of the sensor are included, together with their relationships to measurement accuracy and the effects of solids' velocity profiles. The experimental evaluation of a 53 mm bore sensing head is described, including trials using a calibrated pneumatic conveyor circulating pulverized fuel and cement. Comparisons of test results with the mathematical models of the sensor are used to identify important aspects of the instrument design. Off-line test results obtained using gravity-fed solids flow show that the system repeatability is within +/-0.5% over the velocity range of 2-4 m s(-1) for volumetric concentrations of solids no greater than 0.2%. Results obtained in the pilot-plant trials demonstrate that the system is capable of achieving repeatability better than +/-2% and linearity within +/-2% over the velocity range 20-40 m s(-1) for volumetric concentrations of solids in the range 0.01-0.44%.

The CO2 microalgae biorefinery: high value products and biofuels using halophilic microalgae in the “D-Factory”

Fuel-only algal systems are not economically feasible because yields are too low and costs too high for producing microalgal biomass compared to using agricultural residues e.g. straw. Biorefineries which integrate biomass conversion processes and equipment to produce fuels, power and chemicals from biomass, offer a solution. The CO2 microalgae biorefinery (D-Factory) is a 10 million Euro FP7-funded project which will cultivate the microalga Dunaliella in highly saline non-potable waters in photobioreactors and open raceways and apply biorefinery concepts and European innovations in biomass processing technologies to develop a basket of compounds from Dunaliella biomass, including the high value nutraceutical, β-carotene, and glycerol. Glycerol now finds markets both as a green chemical intermediate and as a biofuel in CHP applications as a result of novel combustion technology. Driving down costs by recovering the entire biomass of Dunaliella cells from saline cultivation water poses one of the many challenges for the D-Factory because Dunaliella cells are both motile, and do not possess an external cell wall, making them highly susceptible to cell rupture. Controlling expression of desired metabolic pathways to deliver the desired portfolio of compounds flexibly and sustainably to meet market demand is another. The first prototype D-Factory in Europe will be operational in 48 months, and will serve as a robust manifestation of the business case for global investment in algae biorefineries and in large-scale production of microalgae.

Exploring multi-sensor satellite synergies to provide direction to high-resolution along-track altimtery currents

A new approach to compute along-track velocity components by combining altimetry-based across-track components and front directions from remote sensing maps of surface chlorophyll concentration is proposed. The analysis focuses on the South Madagascar region characterized by the strong East Madagascar Current and sharp gradients of surface tracers. The results are compared against in-situ observations from three moorings along the Jason-1 track 196. Accurate information on the total velocity direction is the key factor for obtaining accurate estimates of along-track velocities. Although with some limitations, surface tracer fronts can be successfully used to retrieve such information.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

A new look at the pulsating DB white dwarf GD 358: Line-of-sight velocity measurements and constraints on model atmospheres

We report on our findings of the bright, pulsating, helium atmosphere white dwarf GD 358, based on time-resolved optical spectrophotometry. We identify 5 real pulsation modes and at least 6 combination modes at frequencies consistent with those found in previous observations. The measured Doppler shifts from our spectra show variations with amplitudes of up to 5.5 km s-1 at the frequencies inferred from the flux variations. We conclude that these are variations in the line-of-sight velocities associated with the pulsational motion. We use the observed flux and velocity amplitudes and phases to test theoretical predictions within the convective driving framework, and compare these with similar observations of the hydrogen atmosphere white dwarf pulsators (DAVs). The wavelength dependence of the fractional pulsation amplitudes (chromatic amplitudes) allows us to conclude that all five real modes share the same spherical degree, most likely, l=1. This is consistent with previous identifications based solely on photometry. We find that a high signal-to-noise mean spectrum on its own is not enough to determine the atmospheric parameters and that there are small but significant discrepancies between the observations and model atmospheres. The source of these remains to be identified. While we infer Teff =24 kK and log g ~ 8.0 from the mean spectrum, the chromatic amplitudes, which are a measure of the derivative of the flux with respect to the temperature, unambiguously favour a higher effective temperature, 27 kK, which is more in line with independent determinations from ultra-violet spectra.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

High-frequency oscillations in a solar active region observed with the RAPID DUAL IMAGER

High-cadence, synchronized, multiwavelength optical observations of a solar active region (NOAA 10794) are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system: the rapid dual imager. Wavelet analysis is undertaken to search for intensity related oscillatory signatures, and periodicities ranging from 20 to 370 s are found with significance levels exceeding 95%. Observations in the H-α blue wing show more penumbral oscillatory phenomena when compared to simultaneous G-band observations. The H-α oscillations are interpreted as the signatures of plasma motions with a mean velocity of 20 km s-1. The strong oscillatory power over H-α blue-wing and G-band penumbral bright grains is an indication of the Evershed flow with frequencies higher than previously reported.

A high-resolution soft x-ray spectrometer on the MAST tokamak

A curved crystal spectrometer in Johann configuration has been implemented on MAST to obtain values for electron temperature, ion temperature and toroidal velocity. The spectrometer is used to examine medium Z impurities in the soft x-ray region by utilising a Silicon (111) crystal, bent using a 4 pin bending jig, and a CCD detector (Deltat=8 ms). Helium-like Argon emissions from 3.94 to 4.00 Angstrom have been examined using a crystal radius of 859.77 mm. The Bragg angle and crystal radius can be adjusted with relative ease. The spectrometer can be scanned toroidally and poloidally to include a radial view which facilitates absolute velocity measurements by assuming radial velocity =0. Doppler shifts of 2.3x10(-5) Angstrom (1.8 kms(-1)) can be measured. The line of sight is shared with a neutral particle analyzer, which enables in situ ion temperature comparisons. Ray tracing has been used for the development of new imaging spectrometers, using spherical/toroidal crystals, planned to be implemented on MAST. (C) 2004 American Institute of Physics.

In-situ FT-IR spectroscopic studies of fuel cell electro-catalysis: From single-crystal to nanoparticle surfaces

The work presented in this article shows the power of the variable temperature, in-situ FT-IR spectroscopy system developed in Newcastle with respect to the investigation of fuel cell electro-catalysis. On the Ru(0001) electrode surface, CO co-adsorbs with the oxygen-containing adlayers to form mixed [CO+(2x2)-O(H)] domains. The electro-oxidation of the Ru(0001) surface leads to the formation of active (1x1)-O(H) domains, and the oxidation of adsorbed CO then takes place at the perimeter of these domains. At 20 degrees C, the adsorbed CO is present as rather compact islands. In contrast, at 60 degrees C, the COads is present as a relatively looser and weaker adlayer. Higher temperature was also found to facilitate the surface diffusion and oxidation of COads. No dissociation or electro-oxidation of methanol was observed at potentials below approximately 950mV; however, the Ru(0001) surface at high anodic potentials was observed to be very active. On both Pt and PtRu nanoparticle surfaces, only one linear bond CO adsorbate was formed from methanol adsorption, and the PtRu surface significantly promoted both methanol dissociative adsorption to CO and its further oxidation to CO2. Increasing temperature from 20 to 60 degrees C significantly facilitates the methanol turnover to CO2.

SiO in G34.26: Outflows and shocks in a high mass star forming region

We have looked for SiO emission as evidence of shocks in the high mass star formation region G34.26+0.15. JCMT, VLA and FCRAO observations show that SiO emission is widespread across the region. The SiO emission highlights a massive, collimated out ow and other regions where stellar winds are interacting with molecular clumps. As in other star forming regions, there is also SiO at ambient velocities which is related to the out ow activity. No strong SiO abundance enhancement was measured in either the out ow or the low velocity gas, though abundances up to 10(-8) are possible if the SiO is locally enhanced in clumps and optically thick. SiO emission is not detected from the hot core itself, indicating either that SiO is not strongly enhanced in the hot core or that column densities in the region where grain mantle evaporation has taken place are low. In line of sight spiral arm clouds, we measure a SiO abundance of 0.4-2 x 10(-10), consistent with previous estimates for quiescent clouds.