Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. All rights reserved.

Axisymmetric shock wave interaction with a cone: a benchmark test

Results of the benchmark test are presented of comparing numerical schemes solving shock wave of M-s = 2.38 in nitrogen and argon interacting with a 43 degrees semi-apex angle cone and corresponding experiments. The benchmark test was announced in Shock Waves Vol. 12, No. 4, in which we tried to clarify the effects of viscosity and heat conductivity on shock reflection in conical flows. This paper summarizes results of ten numerical and two experimental applications. State of the art in studies regarding the shock/cone interaction is clarified.

MBCNS2 Example Book

MBCNS2 is a small collection of programs for the simulation of transient two-dimensional (or axisymmetric) flows. It is part of the larger collection of compressible flow simulation codes found at http://www.mech.uq.edu.au/cfcfd/. This manual is a collection of example simulations: scripts, results and commentary. It may be convenient for new users of the code to identify an example close to the situation that they wish to model and then adapt the scripts for that example.

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Characterization of turbulent jets from high-aspect-ratio rectangular nozzles

Turbulent free jets issuing from rectangular slots with various high aspect ratios (15-120) are characterized. The centerline mean and rms velocities are measured using hot-wire anemometry over a downstream distance of up to 160 slot heights at a slot-height-based Reynolds number of 10000. Experimental results suggest that a rectangular jet with sufficiently high aspect ratio (> 15) may be distinguished between three flow zones: an initial quasi-plane-jet zone, a transition zone, and a final quasi-axisymmetric-jet zone. In the quasi-plane-jet zone, the turbulent velocity field is statistically similar, but not identical, to those of a plane jet. (c) 2005 American Institute of Physics.

On the mechanism of microbial cell disruption in high-pressure homogenisation

The tensions produced in the wall of a rigid, thin-walled, liquid-filled sphere as it moves with an axisymmetric straining flow are examined. This problem has not been previously addressed. A generalised correlation for the maximum wall tension, expressed in dimensionless form as a Weber number (We), is developed in terms of the acceleration number (Ac) and Reynolds number (Re) of the straining flow. At low Reynolds number We is dominated by viscous forces, while inertial forces due to internal pressure gradients caused by sphere acceleration dominate at higher Re. The generalised correlation has been used to examine the case of a typical yeast cell (a thin-walled, liquid-filled sphere) passing through a typical high-pressure homogeniser (a straining-flow device). At 56 MPa homogenising pressure, a 6 mu m yeast cell experiences tensions in the inertially dominated regime (Re = 100). The correlation gives We = 0.206, corresponding to a maximum wall tension of 8 Nm(-1). This is equivalent to an applied compressive force of 150 mu N and compares favourably with the force required to break yeast cells under compressive micromanipulation (40-90 mu N). Inertial forces may therefore be an important and previously unrecognised. mechanism of microbial cell disruption during high-pressure homogenisation. Further work is required to examine the likelihood of cell deformation in the high-strain-rate short-residence-time environment of the homogeniser, and the effect that such deformation may have on the contribution of inertial forces to disruption. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

A small disturbance in the strong vortex flow

A small disturbance in the axisymmetric, bathtub-like flow with strong vorticity is considered and the asymptotic representation of the solution is found. It is shown that if the disturbance is smaller than a certain critical scale, the conventional strong vortex approximation cannot describe the field generated by the disturbance not only in the vicinity of the disturbance but also at the distances much larger than the critical scale. (C) 2001 American Institute of Physics.

Moderately strong vorticity in a bathtub-type flow

Theoretical and numerical analysis is performed for an inviscid axisymmetric vortical bathtub-type flow. The level of vorticity is kept high so that the image of the flow on the radial-axial plane (r-z plane) is not potential. The most significant findings are: (1) the region of validity of the strong vortex approximation is separated from the drain by a buffer region, (2) the power-law asymptote of the stream function, specified by Delta Psi similar to r(4/3) Deltaz, appears near the axis when vorticity in the flow is sufficiently strong and (3) the local Rossby number in the region of the 4/3 power-law the initial vorticity level in the flow and the global Rossby number.

Control of chaotic instability in a dual-spin spacecraft with dissipation using energy methods

Control of chaotic instability in a rotating multibody system in the form of a dual-spin spacecraft with an axial nutational damper is achieved using an algorithm derived using energy methods. The control method is implemented on two realistic spacecraft parameter configurations which have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude and consequently impair pointing accuracy. The control method is formulated from nutational stability results derived using an energy sink approximation for a dual-spin spacecraft with an asymmetric platform and axisymmetric rotor. The effectiveness of the control method is shown numerically and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and Bifurcation diagrams.

Investigation of the flow field of a highly heated jet of air

Measurements of mean and fluctuating velocity and temperature and their self- and cross-products to the third-order are presented for a heated axisymmetric air jet. Froude numbers in the range of 3500 13,190, Reynolds numbers in the range of 3470-8500 and non-dimensional streamwise distances. X*, from 0.27 to 1.98 are covered by the data. It was found that turbulence intensity decreases for the heated jet in the region between the inertia dominated and the buoyancy dominated regions which is contrary to findings with helium jets mixing with ambient air to produce density fluctuations. The effects of heating on the turbulent kinetic energy budget and the temperature variance budget show small differences for the inertia dominated region and the intermediate region which help to explain the transition process to the far field plume region. Constants are evaluated for the isotropic eddy diffusivity and generalised gradient hypothesis models as well as the scalar variance model. No significant effect of heating on the dissipation time-scale ratio was found. A novel wire array with an inclined cold wire was used. Measurements obtained with this probe are found to lead to asymmetries in some of the higher-order products. Further investigation suggested that the asymmetries are attributable to an as yet unreported interference effect produced by the leading prong of the inclined temperature wire, The effect may also have implications for inclined velocity wires which contain a temperature component when used in heated flows. (C) 2002 Elsevier Science Inc. All rights reserved.

Drop penetration into porous powder beds

The kinetics of drop penetration were studied by filming single drops of several different fluids (water, PEG200, PEG600, and HPC solutions) as they penetrated into loosely packed beds of glass ballotini, lactose, zinc oxide, and titanium dioxide powders. Measured times ranged from 0.45 to 126 s and depended on the powder particle size,viscosity, surface tensions, and contact angle. The experimental drop penetration times were compared to existing theoretical predictions by M. Denesuk et al. (J. Colloid Interface Sci. 158, 114, 1993) and S. Middleman (Modeling Axisymmetric Flows: Dynamics of Films, Jets, and Drops, Academic Press, San Diego, 1995) but did not agree. Loosely packed powder beds tend to have a heterogeneous bed structure containing large macrovoids which do not participate in liquid flow but are included implicitly in the existing approach to estimating powder pore size. A new two-phase model was proposed where the total volume of the macrovoids was assumed to be the difference between the bed porosity and the tap porosity. A new parameter, the effective porosity (epsilon)eff, was defined as the tap porosity multiplied by the fraction of pores that terminate at a macrovoid and are effectively blocked pores. The improved drop penetration model was much more successful at estimating the drop penetration time on all powders and the predicted times were generally within an order of magnitude of the experimental results. (C) 2002 Elsevier Science (USA).

Effect of cold wire inclination on temperature fluctuation measurements in a heated jet

Hot-wire anemometers at low operating currents are used as fast response resistance thermometers for the study of heated turbulent flows. Simultaneous measurement of temperature and velocity is generally performed with multi-wire arrays. In order to give good spatial resolution a new layout has been tested which uses an inclined temperature wire positioned parallel to the nearest inclined velocity wire. This leads to an asymmetric wire arrangement relative to the mean flow direction. As expected, a reduction in thermal interference from the velocity wires results when compared with an array containing a temperature wire placed normal to the flow. However, measurement of higher order moments of fluctuating quantities in an axisymmetric jet shows considerable distortion of radial distributions which is traced to alteration of the temperature field sensed by the temperature wire. When inclined velocity sensitive wires contain a temperature component, the latter may be affected by the same phenomenon.

An experimental and finite element poroelastic creep response analysis of an intervertebral hydrogel disc model in axial compression

A hydrogel intervertebral disc (lVD) model consisting of an inner nucleus core and an outer anulus ring was manufactured from 30 and 35% by weight Poly(vinyl alcohol) hydrogel (PVA-H) concentrations and subjected to axial compression in between saturated porous endplates at 200 N for 11 h, 30 min. Repeat experiments (n = 4) on different samples (N = 2) show good reproducibility of fluid loss and axial deformation. An axisymmetric nonlinear poroelastic finite element model with variable permeability was developed using commercial finite element software to compare axial deformation and predicted fluid loss with experimental data. The FE predictions indicate differential fluid loss similar to that of biological IVDs, with the nucleus losing more water than the anulus, and there is overall good agreement between experimental and finite element predicted fluid loss. The stress distribution pattern indicates important similarities with the biological lVD that includes stress transference from the nucleus to the anulus upon sustained loading and renders it suitable as a model that can be used in future studies to better understand the role of fluid and stress in biological IVDs. (C) 2005 Springer Science + Business Media, Inc.