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.

Computational modelling of bubbles, droplets and particles in metals reduction and refining

A multi-phase framework is typically required for the CFD modelling of metals reduction processes. Such processes typically involve the interaction of liquid metals, a gas (often air) top space, liquid droplets in the top space and injection of both solid particles and gaseous bubbles into the bath. The exchange of mass, momentum and energy between the phases is fundamental to these processes. Multi-phase algorithms are complex and can be unreliable in terms of either or both convergence behaviour or in the extent to which the physics is captured. In this contribution, we discuss these multi-phase flow issues and describe an example of each of the main “single phase” approaches to modelling this class of problems (i.e., Eulerian–Lagrangian and Eulerian–Eulerian). Their utility is illustrated in the context of two problems – one involving the injection of sparging gases into a steel continuous slab caster and the other based on the development of a novel process for aluminium electrolysis. In the steel caster, the coupling of the Lagrangian tracking of the gas phase with the continuum enables the simulation of the transient motion of the metal–flux interface. The model of the electrolysis process employs a novel method for the calculation of slip velocities of oxygen bubbles, resulting from the dissolution of alumina, which allows the efficiency of the process to be predicted.

Thermoelectric MHD effects on equiaxed crystal morphology

The effects of a constant uniform magnetic field on a growing equiaxed crystal are investigated using a 3-dimensional enthalpy based numerical model. Two cases are considered: The first case looks at unconstrained growth, where the current density is generated through the thermo-electric effect and the current circulates between the tips and roots of the dendrite, the second represents an imposed potential difference across the domain. A jump in the electrical conductivity between the liquid and solid causes the current density to be non uniform. In both cases the resulting Lorentz force drives fluid flow in the liquid phase, this in turn causes advection of the thermal and solute field altering the free energy close to the interface and changing the morphology of the dendrite. In the first case the flow field is complex comprising of many circulations, the morphological changes are modelled using a 2D model with a quasi 3D approximation. The second case is comparable to classic problems involving a constant velocity boundary.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Enhanced Nanoparticle Detection with Liquid Droplet Resonators

Whispering gallery mode particle sensing experiments are commonly performed with solid resonators, whereby the sensing volume is limited to the weak evanescent tail of the mode near the resonator surface. In this work we discuss in detail the sensitivity enhancements achievable in liquid droplet resonators wherein the stronger internal fields and convenient means of particle delivery can be exploited. Asymptotic formulae are derived for the relative resonance shift, line broadening and mode splitting of TE and TM modes in liquid droplet resonators. As a corollary the relative fraction of internal and external mode energy follows, which is shown to govern achievable sensitivity enhancements of solute concentration measurements in droplet sensors. Experimental measurements of nanoparticle concentration based on whispering gallery mode resonance broadening are also presented.

The application of calcined natural dolomitic rock as a solid base catalyst in triglyceride transesterification for biodiesel synthesis

Natural dolomitic rock has been investigated in the transesterification of C-4 and C-8 triglycerides and olive oil with a view to determining its viability as a solid base catalyst for use in biodiesel synthesis. XRD reveals that the dolomitic rock comprised 77% dolomite and 23% magnesian calcite. The generation of basic sites requires calcination at 900 degrees C, which increases the surface area and transforms the mineral into MgO nanocrystallites dispersed over CaO particles. Calcined dolomitic rock exhibits high activity towards the liquid phase transesterification of glyceryl tributyrate and trioctanoate, and even olive oil, with methanol for biodiesel production.

Electrooxidation of the iodides [C(4)mim]I, LiI, NaI, KI, RbI, and CsI in the room temperature ionic liquid [C(4)mim][NTf2]

The electrochemical oxidation of 1-butyl-3-methylimidazolium iodide, [C(4)mim]I, has been investigated by cyclic voltammetry at a platinum microelectrode at varying concentrations in the RTIL 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(4)mim][NTf2]. Two oxidation peaks were observed. The first peak is assigned to the oxidation of iodide to triiodide, in an overall two-electron process: 3I(-)- 2e(-) -> I-3(-). At higher potentials, the electrogenerated triiodide oxidizes to iodine, in an overall one-electron process: I-3(-) - e(-) -> 3/2I(2). An average diffusion coefficient, D, for I- of 1.55 x 10(-11) m(2) s(-1) was obtained. A digital simulation program was used to simulate the voltammetric response, and kinetic parameters were successfully extracted. The parameters deduced from the simulation include D for I-, I-3(-), and I-2 and K-eq,K-2, the equilibrium constant for the reaction of iodide and iodine to form triiodide. Values for these parameters are of the same order as those previously published for the oxidation of Br- in the same RTIL [Allen et al. J. Electroanal. Chem. 2005, 575, 311]. Next, the cyclic voltammetry of five different inorganic iodide salts was studied by dissolving small amounts of the solid in [C(4)mim][NTf2]. Similar oxidation peaks were observed, revealing diffusion coefficients of ca. 0.55, 1.14, 1.23, 1.44, and 1.33 x 10(-11) m(2) s(-1) and solubilities of 714, 246, 54, 83, and 36 mM for LiI, NaI, KI, RbI, and CsI, respectively. The slightly smaller diffusion coefficients for the XI salts (compared to [C(4)mim]I) may indicate that I- is ion-paired with Li+, Na+, K+, Rb+, and Cs+ in the RTIL medium.

The electrochemical oxidation and reduction of nitrate ions in the room temperature ionic liquid [C(2)mim][NTf2]; the latter behaves as a 'melt' rather than an 'organic solvent'

The electrochemical oxidation of 1-butyl-3-methylimidazolium nitrate [C(4)mim][NO3] was studied by cyclic voltammetry in the room temperature ionic liquid (RTIL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl) imide [C(2)mim][NTf2]. A sharp peak was observed on a Pt microelectrode (d = 10 mu m), and a diffusion coefficient at infinite dilution of ca. 2.0 x 10(-11) m(2) s(-1) was obtained. Next, the cyclic voltammetry of sodium nitrate (NaNO3) and potassium nitrate (KNO3) was studied, by dissolving small amounts of solid into the RTIL [ C2mim][ NTf2]. Similar oxidation peaks were observed, revealing diffusion coefficients of ca. 8.8 and 9.0 x 10(-12) m(2) s(-1) and solubilities of 11.9 and 10.8 mM for NaNO3 and KNO3, respectively. The smaller diffusion coefficients for NaNO3 and KNO3 (compared to [C(4)mim][NO3]) may indicate that NO3- is ion-paired with Na+ or K+. This work may have applications in the electroanalytical determination of nitrate in RTIL solutions. Furthermore, a reduction feature was observed for both NaNO3 and KNO3, with additional anodic peaks indicating the formation of oxides, peroxides, superoxides and nitrites. This behaviour is surprisingly similar to that obtained from melts of NaNO3 and KNO3 at high temperatures ( ca. 350 - 500 degrees C), and this observation could significantly simplify experimental conditions required to investigate these compounds. We then used X-ray photoelectron spectroscopy (XPS) to suggest that disodium( I) oxide (Na2O), which has found use as a storage compound for hydrogen, was deposited on a Pt electrode surface following the reduction of NaNO3.

Liquid structure of the ionic liquid 1,3-dimethylimidazolium bis{(trifluoromethyl)sulfonyl}amide

Neutron diffraction has been used to determine the liquid structure of 1,3-dimethylimidazolium bis{( trifluoromethyl) sulfonyl} amide ([dmim][NTf2]). Significantly smaller charge ordering is found in this liquid compared with analogous chloride and hexafluorophosphate salts due to the diffuse charge density and size of the [NTf2](-) anion. This is manifested in a much larger cation-cation and cation-anion separation and an overlap of the cation-cation and cation-anion shells. Comparison of the liquid structure with the crystal structure reported by Holbrey et al. ( Dalton Trans. 2004, 2267) indicates little correlation, for example, the [NTf2](-) anion adopts a trans orientation predominantly in the liquid whereas a cis orientation is found in the solid phase.

Liquid clathrate formation in ionic liquid-aromatic mixtures

1-Alkyl-3-methylimidazolium containing ionic liquids with hexafluorophosphate, bis(trifyl)imide, tetrafluoroborate, and chloride anions form liquid clathrates when mixed with aromatic hydrocarbons; in the system 1,3-dimethylimidazolium hexafluorophosphate-benzene, the aromatic solute could be trapped in the solid state forming a crystalline 2: 1 inclusion compound.

Evolution of Clay Morphology in Polypropylene/Montmorillonite Nanocomposites upon Equi-biaxial Stretching: A Solid-State NMR and TEM Approach

Solid-state NMR and TEM were used to quantitatively examine the evolution of clay morphology upon equibiaxial stretching of polypropylene/montmorillonite (PP-MMT) nanocomposites up to a stretch ratio (?= final length/initial length) of 3.5. 1 H spin-lattice relaxation times were measured by the saturation-recovery sequence. For the nanocomposites, initial portions of the magnetization recovery
curves (e~20 ms) were found to depend on v t, indicative of diffusion-limited relaxation and in agreement with calculations based on estimates of the spin-diffusion barrier radius surrounding the paramagnetic centers in the clay, the electron-nucleus coupling constant, and the spin-diffusion coefficient. Initial slopes of these magnetization recovery curves directly correlated with the fraction of clay/polymer interface. New clay surface was exposed as a near linear function of strain. Long-time portions of the magnetization recovery curves yielded information on the average interparticle separations, which decreased slowly before reaching a plateau at ?=~2.5 as particles aligned. TEM images supported these findings and were used to define and quantify degrees of exfoliation and homogeneity from the NMR data. Exfoliation, defined as (platelets/ stack)-1, increased from 0.38 (unstretched) to 0.80 at ? = 3.5 for PP-MMT nanocomposites stretched at
150 C and 16 s-1. A lower stretch temperature, 145 C, which is slightly below melting onset, led to an exfoliation degree of 0.87 at ?= 2.8, consistent with the ability of higher melt viscosities to allow for higher shear stress transfer. Exposure of new clay surface is attributed to aggregate breakup and orientation at low strains (? e ~2) and to platelets sliding apart at higher strains.

Measurements of the XUV mass absorption coefficient of an overdense liquid metal

The increase in the XUV mass absorption coefficient of liquid aluminium, produced by high-power-laser shock-compression, is measured using XUV laser radiography. At a photon energy of 63 eV a change in the mass absorption coefficient by up to a factor of similar to2.2 is determined at densities close to twice that of solid and electron temperatures of the order of 1 eV. Comparison with hydrodynamic simulations indicate that the absorption coefficient scales with density as rho (1.3 +/-0.2).

Interaction of Room Temperature Ionic Liquid Solutions with a Cholesterol Bilayer

Water solutions of representative (IC(4)mim][Cl] and [C(4)mim][Tf2N] room temperature ionic liquids (ILs) in contact with a neutral lipid bilayer made of cholesterol molecules has been investigated by molecular dynamics simulations based on an empirical force field model. The results show that both ILs display selective adsorption at the water-cholesterol interface, with partial inclusion of ions into the bilayer. In the case Of [C(4)mim][Cl], the adsorption of ions at the water-cholesterol interface is limited by a sizable bulk solubility of the IL, driven by the high water affinity of [Cl](-). The relatively low Solubility Of [C(4)mim][Tf2N], instead, gives rise to a nearly complete segregation of the IL component on the bilayer, altering its volume, compressibility, and electrostatic environment. The computational results display important similarities to the results of recent experimental measurements for ILs in contact with phospholipid model membranes (see Evans, K. O. Int. J. Mol. Sci. 2008, 9, 498-511 and references therein).