Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

The thermodynamic characteristics of formation of vacancies in carbon subgroup element crystals

A method for calculating the parameters of formation of vacancies in crystals formed by spherically symmetrical atoms was developed. Both quantum effects at low temperatures and the possibility of the delocalization of atoms at high temperatures were studied. The parameters of formation of vacancies in carbon subgroup element crystals C-diam, Si, Ge, α-Sn, and Pb were calculated. The inclusion of the delocalization of atoms was shown to increase the enthalpy, entropy, and volume of vacancy formation. At low temperatures, the parameters of vacancy formation were found to depend strongly on the temperature, and the entropy of vacancy formation became negative. At high temperatures, close agreement with experimental data and theoretical estimates reported by other authors was obtained. The temperature dependence of vacancy parameters was studied for diamond heated isobarically from 100 to 4500 K. The applicability scope of the Arrhenius equation with a temperature-independent activation energy is discussed. The validity of the “compensation rule” (correlation between the entropy and enthalpy of vacancy formation) was demonstrated. It was also shown that the volume and entropy of vacancy formation were correlated over the whole temperature range studied.

     

K X-ray production cross sections for 40-180 keV protons

Considering the importance of the X-ray production cross sections for the determination of the element concentrations in a given material, we have measured them experimentally for Al Si, Sc, Ti, V, Fe, Co, Ni and Cu bombarded by protons with energies ranging from 40 to 180 keV. This revised version was published online in July 2006 with corrections to the Cover Date.     

Unsteady motion in single-line flow-injection systems

Unsteady motion in single-line flow-injection systems consisting of n tubular elements with valve or hydrodynamic injection is studied theoretically. A formula for the duration of the initial period of unsteady flow is derived. This initial period is much shorter than the mean residence time of the analyte in practical systems. It can therefore be neglected in mathematical modelling of such systems. Experimental data obtained with single-line systems with valve and hydrodynamic injection confirmed the validity of the theoretical equations.     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

Structure of solid monolayers and multilayers of n -hexane on graphite

We present all-atom molecular dynamics simulations ofn-hexane on the basal plane of graphite at monolayer and multilayer coverages. In keeping with experimental data, we find the presence of ordered adsorbed layers both at single monolayer coverage and when the adsorbed layer coexists with excess liquid adsorbate. Using a simulation method that does not impose any particular periodicity on the adsorbed layer, we quantitatively compare our results to the results of neutron diffraction experiments and find a structural transition from a uniaxially incommensurate lattice to a fully commensurate structure on increasing the coverage from a monolayer to a multilayer. The zig-zag backbone planes of all the alkane molecules lie parallel to the graphite surface at the multilayer coverage, while a few molecules are observed to attain the perpendicular orientation at monolayer coverage. Dedicated to Professor C N R Rao on his 70th birthday     

Mathematical modelling of the oxidation of nickel titanium alloys

A one-dimensional bulk reaction model for the oxidation of nickeltitanium is formulated, with preferential oxidation of titaniumbeing included. The modelling is directed at the better understandingof the dominant mechanisms involved in the oxidation processand their significance for the biocompatibility of the alloy.Two different regimes for the relative diffusivities of oxygenand the metals are investigated. By assuming fast bulk reactions,different asymptotic structures emerge in different parameterregimes and the resulting models take the form of moving boundaryproblems. Different profiles of nickel concentration are obtained:in particular a nickel-rich layer (observed in practice) ispresent below the oxide/metal interface for the case when oxygenand the metals diffuse at comparable rates.