Atomistic phenomena in dense fluid shock waves

The shock structure problem is one of the classical problems of fluid mechanics and at least for non-reacting dilute gases it has been considered essentially solved. Here we present a few recent findings, to show that this is not the case. There are still new physical effects to be discovered provided that the numerical technique is general enough to not rule them out a priori. While the results have been obtained for dense fluids, some of the effects might also be observable for shocks in dilute gases.     

Metallothionein as a biomarker for mercury in tissues of rat fed orally with cinnabar

Cinnabar, as one of the most widely used mineral drugs in traditional Chinese medicines, has been proven to have prominent curative effects in clinical use for more than 2000 years. But the safety and toxicity of the drug has been under constant debate in clinic usage. Metallothionein (MT) contains about 30% of cysteine in the molecule, and plays an important detoxification role against heavy metals. In this study, it was used as a biomarker to assess mercurial accumulation in rats fed orally with cinnabar. After feeding rats with cinnabar by gastric gavage at different dosages and at different times, the distribution of heavy metals (including mercury, copper and zinc) and MT was investigated among rat tissues, including liver, kidney, heart, brain, testis and blood. Metals and MT determinations were carried out using inductively coupled plasma mass spectrometry (ICP‐MS) and a modified mercury saturation assay technique respectively. The results indicated that mercury was easily accumulated in the tissues of rats exposed to cinnabar, especially in kidney. For example: at a feeding dosage of 5 g kg^?1 (bw) for 4 weeks, the mercury concentrations in kidney were 13, 8.7, 21.6 and 26 times those in liver, testis, brain and heart respectively; and at 2.5 g kg^?1 (bw) for 2 weeks, the mercury concentrations in kidney were 21, 2.1, 3 and 21 times those in liver, testis, brain and heart respectively. In addition, mercury in kidney and liver of all cinnabar groups was significantly higher than that of the control group (P < 0.01). A high positive correlation observed between MT concentrations and mercury levels in both liver and kidney (R² = 0.9299, P < 0.02 for liver; R² = 0.9923, P < 0.0008 for kidney) indicated that MT could be used as a biomarker for mercury in tissues. Copyright © 2004 John Wiley & Sons, Ltd.     

Operator‐splitting method for the analysis of cavitation in liquid‐lubricated herringbone grooved journal bearings

This paper presents an operator‐splitting method (OSM) for the solution of the universal Reynolds equation. Jakobsson–Floberg–Olsson (JFO) pressure conditions are used to study cavitation in liquid‐lubricated journal bearings. The shear flow component of the oil film is first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component is computed by the Galerkin finite element method. Present OSM solutions for slider bearings are in good agreement with available analytical and experimental results. OSM is then applied to herringbone grooved journal bearings. The film pressure, cavitation areas, load capacity and attitude angle are obtained with JFO pressure conditions. The calculated load capacities are in agreement with available experimental data. However, a detailed comparison of the present results with those predicted using Reynolds pressure conditions shows some differences. The numerical results showed that the load capacity and the critical mass of the journal (linear stability indicator) are higher and the attitude angle is lower than those predicted by Reynolds pressure conditions for cases of high eccentricities. Copyright © 2004 John Wiley & Sons, Ltd.     

A co‐ordinate system independent streamwise upwind method for fluid flow computation

A new co‐ordinate invariant streamwise upwind formulation for convection dominated flows is developed. The eddy diffusivity/viscosity is added directly to the equations in order to remove the oscillations in the solution. The equations then can be solved by any high‐order scheme and the solution retains the accuracy of the high‐order scheme. The accuracy and reduced lateral thickness growth rate are demonstrated with several numerical examples, including pure convective flows and lid‐driven cavity flow. The lateral spreading due to the numerical diffusion is controlled by the anisotropic tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

Grid‐free surface vorticity method applied to flow induced vibration of flexible cylinders

In order to study cross flow induced vibration of heat exchanger tube bundles, a new fluid–structure interaction model based on surface vorticity method is proposed. With this model, the vibration of a flexible cylinder is simulated at Re=2.67 × 10⁴, the computational results of the cylinder response, the fluid force, the vibration frequency, and the vorticity map are presented. The numerical results reproduce the amplitude‐limiting and non‐linear (lock‐in) characteristics of flow‐induced vibration. The maximum vibration amplitude as well as its corresponding lock‐in frequency is in good agreement with experimental results. The amplitude of vibration can be as high as 0.88D for the case investigated. As vibration amplitude increases, the amplitude of the lift force also increases. With enhancement of vibration amplitude, the vortex pattern in the near wake changes significantly. This fluid–structure interaction model is further applied to simulate flow‐induced vibration of two tandem cylinders and two side‐by‐side cylinders at similar Reynolds number. Promising and reasonable results and predictions are obtained. It is hopeful that with this relatively simple and computer time saving method, flow induced vibration of a large number of flexible tube bundles can be successfully simulated. Copyright © 2004 John Wiley & Sons, Ltd.     

Inviscid,laminar and turbulent opposed flows

This paper attempts to reproduce numerically previous experimental findings with opposed flows and extends their range to quantify the effects of upstream pipes and nozzles with inviscid, laminar and turbulent flows. The choice of conservation equations, boundary conditions, algorithms for their solution, the degree of grid dependence, numerical diffusion and the validity of numerical approximations are justified with supporting calculations where necessary. The results of all calculations on the stagnation plane show maximum strain rates close to the annular exit from the nozzles and pipes for lower separations and it can be expected that corresponding reacting flows will tend to extinguish in this region with the extinction moving towards the axis. With laminar flows, the maximum strain rate increased with Reynolds number and the maximum values were generally greater than with inviscid flows and smaller than with turbulent flows. With large separations, the strain rates varied less and this explains some results with reacting flows where the extinction appeared to begin on the axis. The turbulent‐flow calculations allowed comparison of three common variants of a two‐equation first‐moment closure. They provided reasonable and useful indications of strain rates but none correctly represented the rms of velocity fluctuations on the axis and close to the stagnation plane. As expected, those designed to deal with this problem produced results in better agreement with experiment but were still imperfect. Copyright © 2004 John Wiley & Sons, Ltd.     

Time resolved spectroscopic NMR imaging using hyperpolarized Xe

We have visualized the melting and dissolution processes of xenon (Xe) ice into different solvents using the methods of nuclear magnetic resonance (NMR) spectroscopy, imaging, and time resolved spectroscopic imaging by means of hyperpolarized ¹²⁹Xe. Starting from the initial condition of a hyperpolarized solid Xe layer frozen on top of an ethanol (ethanol/water) ice block we measured the Xe phase transitions as a function of time and temperature. In the pure ethanol sample, pieces of Xe ice first fall through the viscous ethanol to the bottom of the sample tube and then form a thin layer of liquid Xe/ethanol. The xenon atoms are trapped in this liquid layer up to room temperature and keep their magnetization over a time period of 11 min. In the ethanol/water mixture (80 vol%/20%), most of the polarized Xe liquid first stays on top of the ethanol/water ice block and then starts to penetrate into the pores and cracks of the ethanol/water ice block. In the final stage, nearly all the Xe polarization is in the gas phase above the liquid and trapped inside the pores. NMR spectra of homogeneous samples of pure ethanol containing thermally polarized Xe and the spectroscopic images of the melting process show that very high concentrations of hyperpolarized Xe (about half of the density of liquid Xe) can be stored or delivered in pure ethanol.     

The 7th Asian Symposium on Visualization

The 7th Asian Symposium on Visualization (7ASV) was successfully held in Singapore from 3^rd to 7^th November 2003. This event was originally scheduled from 26^th to 30^th May 2003, but had to be postponed because of the outbreak of SARS (Severe Acute Respiratory Syndrome) in some parts of the world. A total number of 122 participants attended the Symposium of which, 107 were foreign participants from 16 countries (excluding Singapore) worldwide. There were 104 papers covering broad range of topics presented at the Symposium which were delivered in 2×7 long sessions, and 9 keynote papers.