Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

A Linear Stability Analysis Of Large-Prandtl-Number Thermocapillary Liquid Bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4 <= Pr <= 50). We focus on the relationships between the critical Reynolds number Re-c, the azimuthal wave number m, the aspect ratio F and the Prandtl number Pr. A detailed Re-c-Pr stability diagram is given for liquid bridges with various Gamma. In the region of Pr > 1, which has been less studied previously and where Re, has been usually believed to decrease with the increase of Pr, we found Re-c exhibits an early increase for liquid bridges with Gamma around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the "effective" part of liquid bridge. (c) 2008 Published by Elsevier Ltd on behalf of COSPAR.

Hydrothermal Instability Of Thermocapillary Convection In Large-Prandtl-Number Liquid Bridges Under Microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q-R. method. The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Effects of Mach number on turbulent separation behaviours induced by blunt fin

An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M-infinity = 7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin.

Parkes problems revisited - Application of response number

Response number R-n(n), proposed in [3, 4], is an important independent dimensionless number for the dynamic response of structures [2]. In this paper, the response number is applied to the dynamic plastic response of the well-known Parkes' problem, i.e., beams struck by concentrated mass.