Electrohydrodynamic instability in some nematic cyanobiphenyls in an a.c. electric field

Electrohydrodynamic instability in homeotropically oriented nematic samples of 4'-n-octyl-4-cyano-biphenyl and 4'n-alkyloxy-4-cyanobiphenyl, (n = 8.9) have been studied in an a.c. electric field. The domain patterns during the instability in these compounds in a very low frequency a.c. field are very similar to those in a d.c. field. The domain patterns observed at higher frequencies have been identified as 'maltese crosses' or 'crossed isogyres'. The electro-convective 'isotropic' flows near the electrode play an important role in the observed instability.     

Combined supercritical-fluid chromatography/mass spectrometry in the analysis of diuron in plasma using on-line phase-system switching

On-line sample pretreatment by means of the phase-system switching approach is an interesting technique for the analysis of aqueous samples, e.g., plasma, by means of supercritical-fluid chromatography. In order to analyse plasma samples the following analytical procedure is used. The plasma sample is injected on to a short precolumn, which is washed with water and subsequently dried with nitrogen. Next, the solutes are desorbed with the supercritical mobile phase, analysed with packed-column supercritical-fluid chromatography and detected with either a UV detector or a mass spectrometer, equipped with a moving-belt interface. The herbicide diuron is selected as a test compound to study the feasibility of this approach. Using a selective detector the procedure is sufficiently sensitive to detect diuron in plasma, but not appropriate to detect the diuron metabolites in a post-mortem plasma sample. These have been identified with liquid chromatography/mass spectrometry. The detection limit of diuron in plasma using the procedure described is about 30 ng/mL.     

Characterization of thermal diffusion in polymer solutions by thermal field-flow fractionation: Dependence on polymer and solvent parameters

The thermal diffusion coefficient D_T has been obtained for 17 polymer-solvent combinations, each of them spanning a range of polymer molecular weights, using thermal field-flow fractionation. The polymers examined include polystyrene, poly(alpha-methyl)styrene, polymethylmethacrylate, and polysioprene. The solvents include benzene, toluene, ethylbenzene, tetrahydrofuran, methylethylketone, ethylacetate, and cyclohexane. Although D_T was confirmed as essentially independent of polymer molecular weight, it was found to vary substantially with the chemical composition of polymer and solvent. The results were used to evaluate several thermal diffusion theories; the agreement with theory was generally found to be unsatisfactory. Attempts were then made to correlate the measured thermal diffusion coefficients with various physicochemical parameters of the polymers and solvent. A good correlation was found in which D_T increases with the thermal conductivity difference of the polymer and solvent and varies inversely with the activation energy of viscous flow of the solvent.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of Conditions of Electric-Spark Dispersion of Zinc on the Quality of the Product Obtained

Influence exerted by additives introduced into the working solution, by iron or copper formedin electric-spark dispersion of these metals and zinc, on the quality of the products obtained in electric erosionwas studied.     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulations of a filament in a flowing soap film

Experiments concerning the properties of soap films have recently been carried out and these systems have been proposed as experimental versions of theoretical two‐dimensional liquids. A silk filament introduced into a flowing soap film, was seen to demonstrate various stable modes, and these were, namely, a mode in which the filament oscillates and one in which the filament is stationary and aligns with the flow of the liquid. The system could be forced from the oscillatory mode into the non‐ oscillatory mode by varying the length of the filament. In this article we use numerical and computational techniques in order to simulate the strongly coupled behaviour of the filament and the fluid. Preliminary results are presented for the specific case in which the filament is seen to oscillate continuously for the duration of our simulation. We also find that the filament oscillations are strongly suppressed when we reduce the effective length of the filament. We believe that these results are reminiscent of the different oscillatory and non‐oscillatory modes observed in experiment. The numerical solutions show that, in contrast to experiment, vortices are created at the leading edge of the filament and are preferentially grown in the curvature of the filament and are eventually released from the trailing edge of the filament. In a similar manner to oscillating hydrofoils, it seems that the oscillating filaments are in a minimal energy state, extracting sufficient energy from the fluid to oscillate. In comparing numerical and experimental results it is possible that the soap film does have an effect on the fluid flow especially in the boundary layer where surface tension forces are large. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.