Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Comment on 'processing and microstructure of single grain uranium-doped Y-Ba-Cu-O superconductor'

We comment on the paper by N Hari Babu et al. (2002 Supercond. Sci. Technol. 15 104-10) and point out misinterpretations of the chemical composition of U-bearing deposits observed in Y123. The observed small deposits are those of new compounds which do not contain Cu, rather than refined Y211 plus U, as stated by the authors. We further note that extensive literature, not quoted, is in disagreement by nearly an order of magnitude concerning the values of Pt and U doping at which the optimum value of Jc is obtained. Other related information, presently in the literature, which may be helpful to those working with this high temperature superconducting chemical system, is presented.

Development of the microstructure of uranium-doped Nd-Ba-Cu-O

Melt grown Nd-Ba-Cu-O (NdBCO) has been reported to exhibit higher values of critical current density, Jc and irreversibility field, Hirr, than other (RE)BCO superconductors, such as YBCO. The microstructure of NdBCO typically contains 5-10 μm sized inclusions of the Nd4Ba2Cu2O10 phase (Nd-422) in a superconducting NdBa2Cu3O7-δ phase (Nd-123) matrix. The average size of these inclusions is characteristically larger than that of the Y2BaCuO5 (Y-211) inclusions in YBCO. As a result, there is scope to further refine the Nd-422 size to enhance Jc in NdBCO. Large grain samples of NdBCO superconductor doped with various amounts of depleted UO2 and containing excess Nd-422 have been fabricated by top seeded melt growth under reduced oxygen partial pressure. The effect of the addition of depleted UO2 on the NdBCO microstructure has been studied systematically in samples with and without added CeO2. It is observed that the addition of UO2 refines the NdBCO microstructure via the formation of uranium-containing phase particles in the superconducting matrix. These particles are of approximately spherical geometry with dimensions of around 1 μm. The average size of the nonsuperconducting phase particles in the uranium-doped microstructure is an order of magnitude less than their size in un-doped Nd-123 prepared with excess Nd-422. The critical current density of uranium-doped NdBCO is observed to increase significantly compared to the undoped material.

Measurement of the Mechanical Properties of Filtercakes

In order to improve drilling mud design to cater for specific well situations, a more comprehensive knowledge and understanding of filter cake failure is needed. This paper describes experimental techniques aimed at directly probing the mechanical properties of filter cakes, without having to take into account artefacts due to fluid flow in the substrate. The use of rheometers allows us to determine shear yield stress and dynamic shear modulii of cakes grown on filter paper. A new scraping technique measures the strength and moisture profiles of typical filter cakes with a 0.1 mm resolution. This technique also allows us to probe the adhesion between the filter cake and its rock substrate. In addition, œdometer drained consolidation and unloading of a filter cake give us compression parameters useful for Cam Clay modelling. These independent measurements give similar results as to the elastic modulus of different filter cakes, showing an order of magnitude difference between water based and oil based cakes. We find that these standard cakes behave predominantly as purely elastic materials, with a sharp transition into plastic flow, allowing for the determination of a well-defined yield stress. The effect ofsolids loading on a given type of mud is also studied.

The flexoelectric behaviour of a hypertwisted chiral nematic liquid crystal

The flexoelectric behaviour of a hypertwisted chiral nematic bimesogenic liquid crystal is presented. Through detailed electro-optic measurements, with particular emphasis on the switching properties, we demonstrate remarkably high optical axis tilt angles. The material studied possessed a room temperature nematic phase and aligned easily on cooling under the application of a moderate electric field. Switching times of the order of 500 μs and contrast ratios of 90:1 are readily achieved. The tilt angles, measured using the rotating analyser technique, were found to be practically temperature independent and linear with the applied field. Tilt angles of 22.5° were obtained with moderate applied fields of 9.4 V/μm whilst fields of 25 V/μm yielded tilt angles of 45°. We believe these are the highest tilt angles ever recorded for such fields. © 2001 OPA (Overseas Publishers Association) N.V. Published by license under the Gordon and Breach Science Publishers imprint, a member of the Taylor & Francis Group.

A three-dimensional tunnel model for calculation of train-induced ground vibration

The frequency range of interest for ground vibration from underground urban railways is approximately 20 to 100 Hz. For typical soils, the wavelengths of ground vibration in this frequency range are of the order of the spacing of train axles, the tunnel diameter and the distance from the tunnel to nearby building foundations. For accurate modelling, the interactions between these entities therefore have to be taken into account. This paper describes an analytical three-dimensional model for the dynamics of a deep underground railway tunnel of circular cross-section. The tunnel is conceptualised as an infinitely long, thin cylindrical shell surrounded by soil of infinite radial extent. The soil is modelled by means of the wave equations for an elastic continuum. The coupled problem is solved in the frequency domain by Fourier decomposition into ring modes circumferentially and a Fourier transform into the wavenumber domain longitudinally. Numerical results for the tunnel and soil responses due to a normal point load applied to the tunnel invert are presented. The tunnel model is suitable for use in combination with track models to calculate the ground vibration due to excitation by running trains and to evaluate different track configurations. © 2006 Elsevier Ltd. All rights reserved.

Contrasting code performances for computational aeroacoustics of jets

The numerical propagation of subcritical Tollmein-Schlichting (T-S), inviscid vortical and cut-on acoustic waves is explored. For the former case, the performances of the very different NEAT, NTS, HYDRA, FLUXp and OSMIS3D codes is studied. A modest/coarse hexahedral computational grid that starkly shows differences between the different codes and schemes used in them is employed. For the same order of discretization the five codes show similar results. The unstructured codes are found to propagate vortical and acoustic waves well on triangular cell meshes but not the T-S wave. The above code contrasting exercise is then carried out using implicit LES or Smagorinsky LES for and Ma = 0.9 plane jet on modest 0.5 million cell grids moving to circa 5 million cell grids. For this case, even on the coarse grid, for all codes results were generally encouraging. In general, the spread in computational results is less than the spread of the measurements. Interestingly, the finer grid turbulence intensity levels are slightly more under-predicted than those of the coarse grid. This difference is attributed to the numerical dispersion error having a favourable coarse grid influence. For a non-isothermal jet, HYDRA and NTS also give encouraging results. Peak turbulence values along the jet centreline are in better agreement with measurements than for the isothermal jets. Copyright © 2006 by University of Wales.

Electric conduction improvement of well-structured multi-walled carbon nanotubes

We characterized the electrical conductance of well-structured multi-walled carbon nanotubes (MWCNTs) which had post-treated by a rapid vacuum arc thermal annealing process and structure defects in these nanotubes are removed. We found that the after rapid vacuum arc annealing, the conductivity of well-structured MWCNTs can be improved by an order of magnitude. We also investigated the conductivity of MWCNTs bundle by the variation of temperatures. These results show that the conductance of annealed defect-free MWCNTs is sensitive to temperature imply the phonon scatting dominated the electron conductions. Compare to the well-structured MWCNTs, the defect scattering dominated the electron conduction in the as-grown control sample which has large amount of structure defects. A detail measurement of electron conduction from an individual well-structured MWCNT shows that the conductivity increases with temperatures which imply such MWCNTs exhibited semiconductor properties. We also produced back-gated field-effect transistors using these MWCNTs. It shows that the well-structured MWCNT can act as p-type semiconductor. © 2010 IEEE.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number W e = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20 % of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number We = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20% of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

On the stabilisation of turbulent lifted flames

This study explores the stabilisation mechanisms of turbulent lifted flames by examining the scalar dissipation rate (SDR) of both passive and reactive scalars and their cross dissipation (CDR) in the stabilisation region. DNS results of a laboratory scale hydrogen turbulent lifted flame has been used for this analysis. Various definitions of the flame leading edge (FLE) has been compared and differences are illustrated. Time and spatial averaged statistic of SDR and CDR were examined. It was found that the averaged SDR for the mixture fraction at FLE was well below the reference quenching value for stoichiometric mixture. The averaged SDR for the progress variable is in the same order of the unstrained premixed laminar flame value. It was observed that the averaged CDR changed from negative to weakly positive at FLE. The change in sign was explained by a change in the relative alignment of the gradients of mixture fraction and progress variable. It was thus evident that the CDR was a good marker for stabilisation region and an important quantity in stabilisation mechanism.

Calculation and measurement of the magnetostriction of multilayered composites

A new theoretical model that predicts the magnetostriction of multilayered composites has been developed. The model takes into account the shear stress between the composite layers and consequently predicts a non-uniform strain along their thickness. The model has been experimentally validated by producing composites formed from three materials with different magnetostrains and mechanical properties, and controlled layer thicknesses in the order of micrometers. Deformations of several ppm, up to 7.5% of the saturation magnetostrain were measured between the edge and the centre of such composites. © 2006 Elsevier B.V. All rights reserved.

Eigenmodes of lined flow ducts with rigid splices

This paper presents an analytic expression for the acoustic eigenmodes of a cylindrical lined duct with rigid axially running splices in the presence of flow. The cylindrical duct is considered to be uniformly lined except for two symmetrically positioned axially running rigid liner splices. An exact analytic expression for the acoustic pressure eigenmodes is given in terms of an azimuthal Fourier sum, with the Fourier coefficients given by a recurrence relation. Since this expression is derived using a Greens function method, the completeness of the expansion is guaranteed. A numerical procedure is described for solving this recurrence relation, which is found to converge exponentially with respect to number of Fourier terms used and is in practice quick to compute; this is then used to give several numerical examples for both uniform and sheared mean flow. An asymptotic expression is derived to directly calculate the pressure eigenmodes for thin splices. This asymptotic expression is shown to be quantitatively accurate for ducts with very thin splices of less than 1 % unlined area and qualitatively helpful for thicker splices of the order of 6 % unlined area. A thin splice is in some cases shown to increase the damping of certain acoustic modes. The influences of thin splices and thin boundary layers are compared and found to be of comparable magnitude for the parameters considered. Trapped modes at the splices are also identified and investigated. © 2011 Cambridge University Press.

Increasing the design strength of glass – fractography and stress testing

Brittleness is the unintended, but inevitable consequence of producing a transparent ceramic for architectural applications such as the soda-lime glass. Its tensile strength is particularly sensitive to surface imperfections, such as that from natural weathering and malicious damage. Although a significant amount of testing of new glass has been carried out, there has been surprisingly little testing on weathered glass. Due to the variable nature of the causes of surface damage, the lack of data on weathered glass leads to a considerable degree of uncertainty in the long-term strength of exposed glass. This paper presents the results of recent tests on weathered annealed glass which has been exposed to natural weathering for more than 20 years. The tests include experimental investigations using the co-axial ring setup as well as optical and atomic force microscopy of the glass surfaces. The experimental data from these tests is subsequently used to extend existing fracture mechanics-based models to predict the strength of weathered glass. It is shown that using an automated approach based directly on finite element analysis results can give an increase in effective design strength in the order of 70 to 100% when compared to maximum stress methods. It is also shown that by combining microscopy and strength test results, it is possible to quantitatively characterise the damage on glass surfaces.