(465f) Decay of Mechanically Driven Axial Counter-Current in a High Speed Rotating Cylinder Using Dsmc Simulations

The decay of mechanically driven axial counter-current along the axial direction in a high speed rotating cylinder is studied for wall pressure P_w in the range 20 to 100 m-bar using two dimensional Direct Simulation Monte Carlo (DSMC) simulations [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. The shape & magnitude of the radial-profile of the axial mass flux is investigated quantitatively at various axial locations, and the axial-decay is characterized by a universal exponential function with varying exponent & pre-exponential factor based on the wall pressure and hence the hold up. The analysis shows that as the wall pressure is increased from 20 to 100 m-bar, the shift in the inversion point (corresponds to zero axial mas flux) along the axial length is significant ((Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159); (Kumaran & Pradhan, J. Fluid Mech., vol. 753, 2014, pp. 307-359)). The analysis further indicates that the decay of axial counter-current influences both the flow profile efficiency (E_F) and the circulation efficiency (E_C) to a great extent, and plays an important role in deciding the separation performance of the gas centrifuge machine. The DSMC simulation results are compared with the analytical results for the decay length based on Dirac equation of high speed approximation ( Z_D = (1/ 2η) (1/(4.82 A⁶)) ((P_wall M_W)/(R_g T)) (V_θ R²_wall ) [ 1 + (((γ-1) M_W V_θ²)/(4 γ R_g T))² ] ^1/2 , and found good agreement (error within 15%). Here, Z_D is the decay length, η is the gas viscosity, A is the stratification parameter A= (M_W V_θ²/(2 R_g T))^1/2, P_wall is the wall pressure, M_W is the molecular weight, R_g is the universal gas constant, T is the uniform gas temperature, V_θ is the peripheral velocity, R_wall is the radius of the cylinder, γ is the specific heat ratio (C_P/C_V), and the parameter B = (((γ-1) M_W V_θ²)/(4 γ R_g T)) represents the ratio of adiabatic force to angular momentum force.

Keywords: High-speed rotating flow, Mechanically driven axial counter-current, Decay length, DSMC Simulations.

References:

1. PRADHAN, S. & KUMARAN, V. 2011 The generalized Onsager model for the secondary flow in a high-speed rotating cylinder. J. Fluid Mech. 686, 109.

2. KUMARAN, V & PRADHAN, S. 2014 The generalized Onsager model for a binary gas mixture. J. Fluid Mech. 753, 307.

3. OLANDER, D. R. 1981 The theory of uranium enrichment by the gas centrifuge. Prog. Nucl. Energy., 1981, 8, 1 - 33.

4. WOOD, H. G. & SANDERS, G. 1983 Rotating compressible flows with internal sources and sinks. J. Fluid Mech. 127, 299.

5. BIRD, G. A. 1994 Molecular gas dynamics and the direct simulation of gas flows. Clarendon Press, Oxford.

6. BIRD, G. A. 1963 Approach to translational equilibrium in a rigid sphere gas. Physics of fluids 6, 1518.

7. BIRD, G. A. 1987 Direct simulation of high-vorticity gas flows. Phys. Fluids. 30, 364.

8. BIRD, G. A. 1981 Monte-Carlo simulation in an engineering context. In S. S. Fischer, editor, Rarefied Gas Dynamics, proceedings of the 12th international symposium - part-1., 239.

9. BORGNAKKE, C. & LARSEN, P. S. 1975 Statistical collision model for Monte Carlo simulation of polyatomic gas mixture, J. Comp. Phys. 18, 405.

10. WALZ, A. & VOLOSCIUK, K. & SCHUTZ, H. 1983 Numerical investigations of the flow field near a model of a scoop using the vortex transport equations. In Proc. Fifth Workshop on Gases in Strong Rotation, Charlottesrille (ed. H.G.Wood), 425.

11. VINCENTI. 1967 Introduction to Physical Gas Dynamics. Wiley, New York.